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Abstract 



Spontaneous chiral symmetry breaking is one of the most important fea- 
tures in low-energy QCD. The chiral symmetry is expected to be restored at 
very high temperature and/or density Accompanied by the chiral phase tran- 
sition, properties of hadrons will be changed especially near the critical point. 
The study of the phenomena associated with the chiral phase transition will 
give us some clues on the connection between the chiral symmetry and the 
low-energy hadron dynamics. 

We develop the theory based on the hidden local symmetry (HLS) at finite 
temperature, which is an effective field theory of QCD and includes pions and 
vector mesons as the dynamical degrees of freedom, and study the chiral phase 
transition in hot matter. We show that the chiral symmetry is restored as 
the vector manifestation (VM), in which the massless degenerate pion (and 
its flavor partners) and the longitudinal p meson (and its flavor partners) as 
the chiral partner. We also present several predictions based on the VM. 
We estimate the critical temperature T c and show the following phenomena 
near T c : the vector charge susceptibility becomes equal to the axial-vector 
charge susceptibility; the vector dominance of the electromagnetic form factor 
of the pion is largely violated; the pion velocity is close to the speed of light. 
Furthermore, we show that the remnant of the VM can be clearly seen in 
the system of heavy mesons. We expect that the VM and its predictions are 
testable by current and future experiments and the lattice analysis. 
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Chapter 1 
Introduction 



Quantum Chromo Dynamics (QCD) is the theory which governs the dynamics of quarks and 
gluons, and eventually the hadron dynamics. The QCD Lagrangian has an approximate global 
chiral symmetry in the light quark sector. The QCD vacuum, however, holds no longer the 
chiral symmetry which is spontaneously broken down caused by its strongly interacting dy- 
namics. Spontaneous chiral symmetry breaking is one of the most important features in 
low-energy QCD. The phenomena of the light pseudoscalar mesons (the pion and its flavor 
partner), which are regarded as the approximate Nambu-Goldstone bosons associated with 
the symmetry breaking, are well described by the symmetry property in the low-energy re- 
gion. 

For the purpose to clarify the connection between the low-energy phenomena and the chiral 
symmetry, one of the most important techniques to study them is using effective field theories 
(EFTs) in which a systematic low-energy expansion is done, as in the chiral perturbation 
theory (ChPT). The EFTs are constructed based on the chiral symmetry, in which the non- 
perturbative effects are "dressed" around the (bare) quarks and gluons and the fundamental 
degrees of freedom are replaced with the relevant ones in the low-energy region, i.e., hadrons. 
In general, the EFTs have many parameters like hadron masses and coupling constants. Some 
of them can be determined in low-energy limit by the chiral symmetry, however others are not 
completely restricted by only the symmetry. When one performs the matching of the EFT 
to QCD to determine the parameters of the EFT, one can clarify how the hadron physics 
succeeds the dynamics of underlying QCD. Furthermore, including quantum effects into the 
parameters, we can understand the hadron physics based on QCD from many different angles 
and clarify the non-perturbative aspects of QCD. 

However, the mechanism how hadrons acquire their masses as a result of the dynamics 
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is still an open problem of QCD. In order to find out some clues on this issue, it is efficient 
to consider a world where the QCD vacuum holds the chiral symmetry. Although the chiral 
symmetry is spontaneously broken in the real world, it is expected to be restored under some 
extreme conditions, e.g., very high temperature and/or density. Properties of hadrons will be 
changed at finite temperature/density, especially near the critical point of the chiral symmetry 
restoration El El EH 13 

Both theoretically and experimentally, it is important to investigate the physics near the 
phase transition point. In fact, the CERN Super Proton Synchrotron (SPS) observed an 
enhancement of dielectron (e + e~) mass spectra below the p/u resonance jH]. This can be 
explained by the dropping mass of the p meson (see, e.g., Refs. j^J El El) associated with 
the chiral symmetry restoration following the Brown- Rho scaling proposed in Ref. 03] . The 
Relativistic Heavy Ion Collider (RHIC) has started to measure several physical processes in 
hot matter which include the dilepton energy spectra. In near future, the CERN Large Hadron 
Collider (LHC) and GSI are planning to measure them. These experiments will further clarify 
the properties of vector mesons in hot matter. 

As indicated by the SPS data, the vector meson is the one of important degrees of freedom 
in studying the phenomena near the QCD phase transition. The theory based on the hidden 
local symmetry jTJ ^] is an EFT including both pions and vector mesons, where the low- 
energy expansion is systematically done (for details, see Ref. [E]). The HLS theory has a wider 
energy range of the validity than the ordinary ChPT. In this thesis, we review the construction 
of the HLS theory in hot matter and its application to the chiral phase transition following 

Refs. [ii [na eii nil El- 

There is no strong restriction concerning vector meson masses in the standard scenario of 
chiral symmetry restoration, where the pion joins with the scalar meson in the same chiral 
representation [see for example, Refs. ED]]- However there is a scenario which cer- 
tainly requires the dropping mass of the vector meson and supports the Brown-Rho scaling: 
In Ref. j2Ij, it was proposed that there can be another possibility for the pattern of chiral 
symmetry restoration, the vector manifestation (VM). The VM was proposed as a novel man- 
ifestation of the chiral symmetry in the Wigner realization, in which the chiral symmetry is 
restored by the massless degenerate pion (and its flavor partners) and the longitudinal p meson 
(and its flavor partners) as the chiral partner. In terms of the chiral representations of the 
low-lying mesons, there is a representation mixing in the vacuum |2T| IT3] . Approaching the 
critical point, we find that there are two possibilities for the pattern of the chiral symmetry 
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restoration (see chapter EJ): One possible pattern is the standard scenario and another is the 
VM. Both of them are on an equal footing with each other in terms of the chiral representa- 
tions. It is worthwhile to study the physics associated with the VM as well as that with the 
standard scenario of the chiral symmetry restoration. 

It has been shown that the VM is formulated at a large number of flavor [21], critical 
temperature and critical density [22] by using the HLS theory, where a second order or 
weakly first order phase transition was assumed. In the VM at finite temperature and/or 
density, the intrinsic temperature and/or density dependences of the bare parameters of the 
HLS Lagrangian played important roles to realize the chiral symmetry restoration consistently 
(see chapters 0] and 03): Since the HLS theory is an EFT of underlying QCD, the parameters 
of the HLS Lagrangian do depend on the temperature and/or density. In the framework of 
the HLS the equality between the axial-vector and vector current correlators at the critical 
point can be satisfied only if the intrinsic thermal and/or density effects are included. The 
intrinsic temperature and/or density dependences are nothing but the information integrated 
out which is converted from the underlying QCD through the Wilsonian matching jTjU^]. In 
other words, the intrinsic effects are the signature that the hadron has an internal structure 
constructed from the quarks and gluons. 

The VM explains an anomalous enhancement of dielectron mass spectra observed at SPS. 
This is a theoretical support of the dropping mass of p meson following the Brown- Rho scaling, 
which can explain the SPS data. 

In this thesis, we assume that only the pions and vector mesons are the relevant (light) 
degrees of freedom until near the critical temperature, say T c — e, and that the chiral symmetry 
is partially restored already at T c — e. Then we study whether the dropping p can be formulated 
in some EFT. We should note that, although we assume the dropping p meson mass m p ~ O(e') 
as an input, it is non-trivial that the dropping p can actually be formulated in the framework 
of quantum field theory. 

Now, what happens to other mesons? In Fig. 11.11 we sketch our view of the temperature 
dependences of the light mesons masses. We assume that the A\ meson is still heavy at T c — e. 
The scalar mesons and their fate near T c may be more controversial. Several pictures on the 
construction of the scalar meson are proposed [see e.g., Refs. I24[ 1251 125] ]- In our present 
picture, we do not introduce the scalar mesons as the dynamical degrees of freedom into the 
EFT because of the following naive speculations based on two kinds of picture on the quark 
content of scalar mesons. 

Two-quark state : The light scalar meson is a bound state of a quark and an anti-quark. The 
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Figure 1.1: Our picture. The solid lines denote the temperature dependences of light mesons, 
7i,a,p,Ai and /o(1370). Throughout this thesis, we concentrate on the region near T c — e 
(denoted by the green ellipsis) where we assume that only the pions and vector mesons are 
relevant degrees of freedom. 

sigma meson mass m a ~ 600 MeV in the real world can be explained by the existence of 
the disconnected quark diagram joined by the gluon exchange |27j . At finite temperature, 
we expect that the effective QCD coupling constant becomes weaker than that in the 
vacuum and that m a may go up. In that case, the m a may go to zero in the limit T — > T c 
although it might be still heavy until around T c — e. There may be a possibility that the 
/o(1370) (or the mixture of a and / (1370)) goes down at T — > T c . 

Four-quark state : The light scalar meson is a bound state of a diquark and an anti-diquark. 
The sigma meson may be melted at finite temperature less than T c since the binding 
energy of qq is smaller than that of qq. When there exists the light scalar meson near T c 
(above T c — e), such scalar meson might be related to the /o(1370) at T = 0. 

In order to understand which universality class the VM belongs to, we may need to include the 
scalar and A\ mesons other than pions and vector mesons for the analysis of critical exponents. 
In the chiral symmetric phase, on the other hand, an existence of some mesonic (soft) modes 
is proposed [see e.g, Refs. [21 12S EH EU W- 

Throughout this thesis, we concentrate on the region out of the window near T c — e where 
we assume that only the pions and vector mesons are relevant degrees of freedom. In that 
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region, we will see that the dropping p can be formulated within the EFT including both 
pions and vector mesons, related to the chiral symmetry restoration. Furthermore, the present 
framework gives several physical predictions which are closely related to the dropping p. These 
predictions can be tested in future experiments and numerical simulations on the lattice. 

Now let us assume that the chiral symmetry of QCD is restored as the VM. Then can we 
see some remnant of the VM in the real world? In order to answer this question, in chapter El 
we study the heavy meson system which consists of one heavy quark and one light anti-quark 
following Ref. [HHj- We construct an effective Lagrangian based on the VM and heavy quark 
symmetry at the restoration point, and then we introduce the effects of the spontaneous chiral 
symmetry breaking to go back to the real world. We will show that the EFT can describe the 
recent experimental data. 

This paper is organized as follows: 

In chapter |2J we give a general concept of EFTs and the matching between an EFT and 
the underlying theory in the Wilsonian sense. Then we apply the Wilsonian matching to the 
QCD in hot and/or dense matter. We give an account of the general idea of the intrinsic 
temperature and/or density dependences of the parameters, which are inevitably introduced 
cis du result of the Wilsonian matching. 

In chapter 01 we briefly review the HLS and the ChPT with HLS in hot matter following 
Refs. [2J [1^1: In section |3~T1 we review the theory based on the HLS and in section we 
show details of the calculation of the hadronic thermal corrections as well as the quantum 
corrections in the background field gauge. In section we construct the axial- vector and 
vector current correlators. Then in section 13.41 we study the temperature dependences of the 
vector meson mass, the pion decay constant and the pion velocity. We also study the validity 
of the vector dominance (VD) in hot matter. In section ETKj we show the bare HLS Lagrangian 
without Lorentz invariance. 

In chapter HJ we briefly review the Wilsonian matching at finite temperature following 
Refs. O 113 EE EH EE] : We extend the Wilsonian matching at zero temperature to the 
version at non-zero temperature. We also study in detail the effect of Lorentz non-invariance 
of the bare pion decay constants. 

In chapter we show how the VM is realized at the critical point in the framework of 
the HLS theory. We also study the predictions of the VM in hot matter: In section 15.11 we 
present the conditions for the bare parameters through the Wilsonian matching at the critical 
temperature. Further, we take into account the Lorentz non-invariance in the bare HLS theory 



CHAPTER 1. INTRODUCTION 



10 



and present the extended conditions. Then we show that the conditions are still protected as 
the fixed point of the RGEs. In section I5~2l we review the formulation of the VM in hot matter 
following Ref. ^Bj- Then in section we show how the pion decay constants vanish at the 
critical temperature following Ref. ^Hj- Next, in section E3J we summarize the predictions of 
the VM in hot matter studied in Refs. ^3J E3 El El E] : in subsection 15.4.11 we estimate the 
value of the critical temperature. In subsection 15.4.21 we focus on the vector and axial- vector 
susceptibilities very near the critical temperature T c . Then we address the issue of what the 
relevant degrees of freedom can be at the critical point. In subsection 15.4. H| we shed some 
light on the validity of the vector dominance (VD) of electromagnetic form factor of the pion 
near the critical temperature. We find that the VM predicts a large violation of the VD at the 
critical temperature. This indicates that the assumption of the VD may need to be weakened, 
at least in some amounts, for consistently including the effect of the dropping mass of the 
vector meson. In subsection 15.4.41 we study the pion velocity at the critical temperature based 
on the VM. There we show the non-renormalization property on the pion velocity at the critical 
temperature, which is protected by the VM. Then using this property and through matching 
to the OPE, we estimate the value of the pion velocity. In section 15 .5[ we discuss the critical 
behavior of m p (T) near the critical temperature. 

In chapter El we study the system of heavy mesons like D meson, following Ref. [33] . We 
construct an EFT based on the VM and the heavy quark symmetry motivated by the recent 
discovery of new D mesons (J p = + and 1 + ) in Babar, CLEO and Belle (HU ESI IHHj- We 
show that the mass splitting between new D mesons and the existing D mesons ( J p = 0~ and 
1~) is directly proportional to the light quark condensate. Further we study the characteristic 
decay modes and give the predictions on the decay widths and the branching ratios. 

We give a brief summary and discussions in chapter 



Chapter 2 

Effective Field Theory and 
Renormalization Group 

Use of effective field theories (EFTs) is one of important implements to study strongly inter- 
acting system like low-energy QCD. An EFT works only in low energy below a characteristic 
scale, which is dependent on dynamics of the EFT, and is expanded in powers of momentum. 
At a given order, the theory has a finite number of counter terms and thus we can perform 
an order-by- order renormalization in the sense of momentum expansion although the EFT is 
conventionally non-renormalized. 

In this chapter, we give a general idea of EFTs and the matching between an EFT and the 
underlying theory in the Wilsonian sense. 

2.1 General Concept of Effective Field Theory 

We consider the QCD Lagrangian with Nf massless quarks: 

4°cd = 1^(9, - WsG^q - ^[G^GH , (2.1) 

with Gfj, and g s being the gluon field and the QCD gauge coupling constant. The field strength 
G^u is defined by 

GW = d,G u - d u G, - ig s [G„ G„] . (2.2) 
We define the left- and right-handed quarks by 

<?l = 7^(1 -7s)<?, to = t;(1 + 7s)<?- (2.3) 

The interacting term between qi and qn never appear in the Lagrangian. It implies that 
the QCD Lagrangian has an SU(Nf)L x SU(Nf)n x U(l)v X U(1)a symmetry. In the real 
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Figure 2.1: Schematic view of effective field theories (EFTs). Here q and g denote quarks and 
gluons respectively. 

world, however, the chiral symmetry is spontaneously broken down by the strongly interacting 
dynamics, which makes pions to be much lighter than other hadrons. The pion are identified 
as massless Nabmu-Goldstone (NG) bosons associated with the spontaneous chiral symmetry 
breaking. 

In the construction of an EFT valid in some energy region, it is important to specify 
the relevant degrees of freedom in the energy scale as well as to respect symmetries of the 
fundamental theory. The relevant degrees of freedom may be changed according to the energy 
scale that one is interested in. We show a schematic view of the relation between EFTs and 
their underlying QCD in Fig. 12.11 In low energies, the relevant degrees of freedom are hadrons. 
Especially pion is the lightest hadron, which is the most relevant degree of freedom in low- 
energy limit, and its dynamics is well described by EFTs based on the chiral symmetry of 
QCD like the chiral perturbation theory (ChPT) jHZ[ EH1 EH]- Here we consider an EFT, say 
EFT-A, which includes only pions as the relevant degrees of freedom. In the energy region 
where the EFT-A works well, other hadrons, e.g., vector mesons, nucleons and excited states, 
are sitting as the background. However at some energy scale, the EFT-A breaks down, and 
then one has to go to a more suitable theory which might be an EFT: As the energy scale 
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increases, other mesons like p mesons are elevated from the background and become relevant 
degrees of freedom as well as pions. The background matter now consists of other hadrons 
except for pions and p mesons. Thus one can switch over the EFT-A to new EFT, EFT-B, 
which includes both pions and p mesons as the relevant degrees of freedom. 
Let us consider the QCD Lagrangian with external source fields: 

£qcd = £g£ D + q(l^ + 757^)9 - Q{S ~ iV)q , (2.4) 

where the external vector (V M ), axial- vector (Ap), scalar (<S) and pseudoscalar (V) fields are 
the hermitian Nf x Nf matrices in flavor space. If there is an explicit chiral symmetry breaking, 
it is introduced as the vacuum expectation value (VEV) of S. In the case of Nf = 3, 

(S) = M = diag(m u ,m d ,m s ). (2.5) 

The following functional of a set of source fields denoted by J, 



Z[J] = j VqVqVGexp % J d 4 x£ Q c D [g, q, G; J] 



(2.6) 



generates Green's functions which have important information of the system. The basic concept 
of the EFT is based on the following assumption: The effective Lagrangian, which has the most 
general form constructed from the chiral symmetry, can give the same generating functional 
as in Eq. (JUJ): 

(2.7) 



Z[J] = J VUexp % J d 4 xC cS [U;J] 



where U denotes the relevant hadronic fields such as the pion fields and £ e g is an Lagrangian 
expressed in terms of these hadrons. In the ChPT, Nf x Nf special-unitary matrix U is 
parametrized by the pseudoscalar fields tt = Tc a T a as 

JJ = e 2in/F % ( 2 _g) 

where F w denotes the pion decay constant. The leading Lagrangian in the chiral limit is given 

by pa SO! 



C^ PT = ^tr[D^D^U], (2.9) 



F 2 
- 

T 

where D^ is the covariant derivative acting on U. 

Corresponding to each Green's function derived from Eq. ([2.6)1 . we have the same Green's 
function obtained from the EFT through Eq. ([2.7jh We should note that according to the 
energy scale, the relevant degrees of freedom are different although both fundamental (q and 
g) and hadronic (U) degrees of freedom are integrated over all configuration space. 
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Figure 2.2: Separation over two regions (hard and soft). A crossover region where both theories 
are available is assumed. In this region, we can perform the matching between two theories. 

The idea mentioned above is the most general construction of EFTs, where only the chiral 
symmetry gives constraints for the forms of the effective Lagrangian. Thus various EFTs 
consistent with the chiral symmetry form an "EFT space" and the real QCD holds a some 
region in the EFT space. 

2.2 Matching in the Wilsonian Sense 

In general, an EFT has many parameters like masses and coupling constants among hadrons. 
We can determine them by performing the matching between two Green's functions, Eqs. (I2.fij) 
and (J2.7|) at a scale. 

In some matching schemes, the renormalized parameters of the EFT are determined from 
QCD. On the other hand, the matching in the Wilsonian sense is performed based on the 
following general idea: 

Step 1 We imagine a some high-energy scale A which separates the perturbative (hard) region 
from the non-perturbative (soft) region [see Fig. l2.2j . Although the quarks and gluons are 
relevant degrees of freedom in the hard region, they are not relevant anymore in the soft 
region. Then we assume that one can replace the quarks and gluons with hadrons since 
"bare" quarks and gluons in the QCD Lagrangian become "dressed" ones by integrating 
out the high-energy quarks and gluons (i.e., the quarks and gluons above A), which 
eventually form hadrons. 

Step 2 Under the above replacement of the degrees of freedom at A, we define the EFT at 
A as the bare EFT. Thus the generating functional derived from the bare EFT leads to 



energy scale 



crossover 
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the same Green's function as that derived from Eq. ()2.6|) at A. The bare parameters of 
the EFT are determined after integrating out the quarks and gluons above A, and then 
carry information of the underlying QCD. 

Once the bare theory is determined through the matching following the above steps, we 
can obtain the physical quantities in low-energy region within the EFT including quantum 
corrections, for instance, by using renormalization group equations (RGEs). 

It is worthwhile mentioning the difference from the way to determine the renormalized 
parameters of the EFT through the matching to QCD. In this scheme, the theory at the 
matching scale A is full quantum (i.e., renormalized) theory which is achieved by starting over 
in low energies. While in the Wilsonian matching, one starts from the bare theory defined 
at A, and in order-by-order renormalizes quantum effects into the parameters to go down 
to low-energy scales. Both of them are almost equivalent in the case that A is much larger 
than a characteristic scale M of the EFT, e.g, energy, temperature and density. In this case, 
corrections 0(M/A) are neglected and we can replace the bare parameter with the renormalized 
one at A and vice versa. 

2.3 Application to Hot and/or Dense Matter 

In this section, we apply the Wilsonian matching to the QCD in hot and/or dense matter. We 
first give an account of the general idea of the intrinsic temperature and/or density dependences 
of the bare parameters, which are inevitably introduced as a result performed the Wilsonian 
matching. We also discuss the effect of Lorentz symmetry violation at bare level. 

In the following, we consider that one performs the Wilsonian matching of the EFT-B shown 
in Fig. 12.11 to perturbative QCD at finite temperature and/or density: In order to obtain 
the bare EFT-B, we integrate out high frequency modes in hot and/or dense matter. The 
hard modes integrated out are converted into the bare parameters and thus the temperature 
and/or density carried by them become temperature and/or density dependences of the bare 
parameters. We call the thermal and/or dense effects as intrinsic temperature and/or density 
dependences. They carry the information without (or behind) the EFT-B, i.e., the information 
of the underlying theory QCD. Ordinary hadronic thermal and/or dense effects are carried by 
the relevant degrees of freedom below A, ir and p, which are generated by dynamics of the 
soft region. Thus in the Wilsonian matching, we can for convenience distinguish two kinds of 
temperature and/or density dependences: 

Intrinsic temperature and/or density dependences : Carried by the relevant degrees of 
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freedom above A (hard region). Such degrees of freedom are integrated out in the EFT-B 
and contribute to bare parameters. 

Hadronic temperature and/or density dependences : Carried by the relevant degrees 
of freedom below A (soft region). Such degrees of freedom interact to heat /particle bath 
and generate thermal and/or dense effects which appear in the Bose/Fermi distribution 
functions. 

In the present case, the hard modes integrated out are all hadronic degrees of freedom except 
for 7T and p relevant below A. It should be noted that the intrinsic effects never appear in the 
bare parameters if one does not take into account these "residual" degrees of freedom. However, 
temperature and/or density of the system are provided by the hadronic matter formed by all 
hadronic degrees of freedom. The matter including only ir andp is not identical with the QCD 
matter. Thus we need to take into account "all the others" never incorporated in system 
within the EFT, which enable us to include the intrinsic effects. 

Further, according to Step 1 mentioned in the previous section, we replace the hard modes 
with the quarks and gluons above A. The replacement indicates the following: The intrinsic 
temperature and/or density dependences are nothing but the signature that the hadron has 
an internal structure constructed from the quarks and gluons. This is similar to the situation 
where the coupling constants among hadrons are replaced with the momentum- dependent 
form factor in higher energy region since the hadronic picture becomes irrelevant. Thus the 
intrinsic thermal and/or dense effects play more important roles in the higher temperature 
region, especially near the phase transition point. 

Intrinsic thermal/dense effects also cause Lorentz non-invariance in the bare theory: At 
finite temperature and/or density, theory has no longer the Lorentz invariance since we specify 
the rest frame of medium. The Lorentz non-invariance of hot and/or dense QCD is generated 
through the dynamics of quarks and gluons. Then the QCD Lagrangian takes the Lorentz 
invariant form. On the other hand, parameters of EFT have implicitly the information of the 
fundamental theory as the intrinsic effects. Thus the intrinsic temperature and/or density 
dependences carry the Lorentz non-invariance of hot/dense QCD, and eventually they cause 
the Lorentz symmetry breaking in the bare EFT. 



Chapter 3 



Hidden Local Symmetry Theory at 
Finite Temperature 

The hidden local symmetry (HLS) theory [TT| IT2] is a natural extension of the non-linear 
chiral Lagrangian and it can describe a system including both pseudoscalars (pions) and vector 
mesons. In the present framework, the vector meson is introduced as the gauge boson of the 
HLS. We can introduce the vector mesons in other frameworks, e.g., matter field [H], massive 
Yang-Mills (121 IIBI El E3 EE] and anti-symmetric tensor field [HBII1ZI- It was shown that these 
models are equivalent with the HLS theory at tree level jUJHEHHIllSniEIlESEB]- However it is 
the most important advantage using the HLS theory that the systematic low-energy expansion 
can be performed. 

In this chapter, we first briefly review the HLS and the chiral perturbation theory (ChPT) 
with HLS. (For details, see Ref. [E].) Next we apply the HLS theory to QCD at finite temper- 
ature following Refs. [TH fT3| EE El EH] • We present the predictions in the low temperature 
region: We show the temperature dependences of the vector meson mass and the pion pa- 
rameters, and the validity of the vector dominance of the electromagnetic form factor of the 
pion. 

3.1 Hidden Local Symmetry 

The HLS theory #1 is based on the G g i bai x -f^iocai symmetry, where G = SU (Nf) L x SU(Nf) r 

is the chiral symmetry and H = SU(Nf)y is the HLS. The basic quantities are the HLS gauge 

^In the modern interpretation |54) . implementing HLS in the chiral Lagrangian can be associated with the 
"ultraviolet completion" to the fundamental theory of strong interactions, i.e., QCD. The matching to QCD 
at a matching scale is therefore a crucial ingredient of the approach. 
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Nonlinear sigma model 
SU(3) L X SU(3) R — SU(3) V 
F 2 

l= ^tr[D^UD^U] 
U = ^ 2 = exp[27ti/F 7l ] 



Other models of vector meson 

* matter field (CCWZ) 

* massive Yang-Mills 

* anti-symmetric tensor field 



♦ 



including vector mesons 



Hidden Local Symmetry (HLS) 
[SU(3) L X SU(3) R ] global Xj^SU^ 

^L,R= exp[ ia/F ] exp[ + i7i/F n ] 



Figure 3.1: Non-linear chiral Lagrangian and the HLS theory, 
boson and two matrix valued variables and £r(x) which transform as 

£l,r{ x ) Z'l,r{ x ) = h i x )iLA x )9L,R , 



(3-1] 



where h(x) G H\ oca \ and gL,R G [SU(A r /)L,R] g i bai- These variables are parameterized as 

£ LR {x) = e *(x)/F* e **(x)/F* 9 (3.2) 



where n = 7i a T a denotes the pseudoscalar Nambu-Goldstone bosons associated with the spon- 
taneous symmetry breaking of G g i bai chiral symmetry, and a = a a T a denotes the Nambu- 
Goldstone bosons associated with the spontaneous breaking of -ffi OC ai * 2 - This a is absorbed 
into the HLS gauge boson through the Higgs mechanism. F n and F a are the decay constants 
of the associated particles. The phenomenoiogicaliy important parameter a is defined as 

a = —^T . (3.3) 

The covariant derivatives of £,lr ar e given by 



9^r - iV^ R + itnlln, 



(3.4) 



^ 2 The a is not to be confused with the scalar that figures in two-flavor linear sigma model. 
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where is the gauge field of H\ oca i, and £ M and TZ^ are the external gauge fields introduced 
by gauging Global symmetry. 

The basic quantities in constructing the Lagrangian are the following two 1-forms: 

% = j. {p^ R ■ 4 + ■ ii) , 

&i* = i (d^r ■ 4 - ■ £j) • (3.5) 

They transform as 

|| — >• • «x || " ^ ( x ) • (3-6) 

When HLS is gauge-fixed to the unitary gauge, o = 0, £l and £r are related with each other 
by 

£ = & = £ = e^- . (3.7) 

This unitary gauge is not preserved under the G g i bai transformation, which in general has the 
following form 

Ggiobai : £ —> £' = £ • g R = 9l ■ £ 

= exp [icr'iir, g R , g^/F^] exp [in'/ F v ] 

= expfm' /F 7r ]exp[-icr / (7r,g'R,g'L)/-P 1 a -] • (3.8) 

The unwanted factor exp [icr'(7r, g R , gi)/F a \ can be eliminated if we simultaneously perform the 
-f^iocai gauge transformation with 

#iocai : h = exp[ia'(iv,g R ,g L )/F a ] = h(iv,g R ,g L ) . (3.9) 

Then the system has a global symmetry G = S77(3)l x S77(3);a under the following combined 
transformation: 

G : £ -> fe(7r, £r, st) • £ • g R = gL ■ £, ■ h(n, gR, ql) ■ (3.10) 
Under this transformation the HLS gauge boson field in the unitary gauge transforms as 

G : —> h(ir,g R ,g L ) ■ ■ h} (n , g R , g L ) - idh(ir , g R , g L ) ■ h j (ir,g R ,g L ) , (3.11) 

which is nothing but the transformation property of Weinberg's "p meson" [55]. The two 
1-forms dy 1 and aij_ transform as 

"ill -> /l ( 7r ' 9r, 9l) ■ • ^(n, g R , g L ) . (3.12) 
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Then, we can regard these 1-forms as the fields belonging to the chiral representations (1, 8) 
(8, 1) and (1, 8) - (8, 1) under SU(3) L x SU(3) R . 

The HLS Lagrangian with lowest derivative terms at the chiral limit is given by ^TJ Ej 



£( 2 ) = -^ 2 tr[d ±At Q^] + F^trt&i^ajf] 



29 



(3.13) 



where g is the HLS gauge coupling, V pu is the field strength of V M and 



(3.14) 



We expand the Lagrangian ()3.13|) in terms of the 7r field taking the unitary gauge of the HLS 
(a = 0) to obtain 



(3.15) 



C(2) = tr [d^n] + ag 2 F 2 tr [p^\ + 2i \-agj tr \pP [d^ , n}} 
- 2 (agF 2 ) tr \p^} +2i(l- |) tr [V [<V , tt]] + ■ ■ • 
where the vector meson field p M is introduced by 

and vector external gauge field V M is defined as * 3 



V„ 



(^ + £. 



(3.16) 



(3.17) 



2 v 

From the expansion in Eq. f)3. 15j) we find the following expressions for the mass of vector meson 
M p , the pirn coupling g pnn , the p-j mixing strength g p and the direct ^titx coupling g 77r7r : 

M n 2 



9 P 



ag 2 F n 2 , 

\ag, 
agF 2 , 

e l - -a). 
v 2 ; 



(3.18) 
(3.19) 
(3.20) 
(3.21) 



By taking the parameter choice a = 2 at tree level, the HLS reproduces the following phe- 
nomenological facts [TT] : 
# 3 Note that the photon field for Nf = 3 is embedded into V M as 

1 2/3 \ 
V M = eA^Q, Q = -1/3 



V 



-1/3 J 



with e being the electromagnetic coupling constant. 
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1- 9 P mr = 9 (universality of the p coupling) 

2. M p 2 = 2g 2 nn F n 2 (the KSRF relation, version II) jEj; 

3. g 17T7T = (vector dominance of the electromagnetic form factor of the pion) 

Moreover, the KSRF relation (version I) j37j is predicted as a low-energy theorem of the 
HLS ISHJ: 

fi'p = 2i^r Qp-K-K-, (3.22) 

which is independent of the parameter a. This relation was first proved at tree level and 
at one-loop level jHO] and then at any loop order [BT] . 

Chiral perturbation theory with HLS 

In the HLS theory it is possible to perform the derivative expansion systematically owing to 
the gauge invariance fHl E21 ES] • In this ChPT with HLS the vector meson mass is considered 
as small compared with the chiral symmetry breaking scale A x , by assigning 0(p) to the HLS 
gauge coupling EH] : 

9 ~ 0{p). (3.23) 

According to the entire list shown in Ref. [BH], there are 35 counter terms at 0(p A ) for general 

Nf. However, only three terms are relevant when we consider two-point functions at chiral 
limit: 

£ (4 ) = zitr[V M „l>1 + ^trLV>H + z 3 tv[V^V^\, (3.24) 

where 

Ap V = l - [6^*4 - ZlC^I] , (3.25) 

V M „ = \ [£ R K^ R + , (3.26) 

with TZ^u and being the field strengths of TZ^ and 

There is a limit that the vector meson mass is small, in large Nf QCD as shown in Ref. |13j . 
Then it can be justified to perform the derivative expansion by m^/A 2 , under such an extreme 
condition. We obtain physical quantities in the real world with Nf = 2 or 3 by extrapolating 
the results in large Nf. Numerically the expansion parameter is not very small, w? p j A 2 X ~ 0.5. 
However, the Wilsonian matching at T = with Nf = 3, through which the parameters of 
the HLS Lagrangian are determined by the underlying QCD at the matching scale A (see 
chapter HJ, was shown to give several predictions in remarkable agreement with experiments 
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|SHll3| . where the validity of the ChPT with HLS is essential. This strongly suggests that the 
extrapolation mentioned above is valid even numerically. 

As discussed in Ref. ^H], A A x can be justified in the N c limit of QCD * 4 : In this limit, 
F%(A) scales as N c , which implies that A x becomes large in the large N c limit. On the other 
hand, the meson masses do not scale with N c , so that we can introduce the A which has no 
large N c scaling property. Thus the quadratic divergent at n-th loop order is suppressed by 
[A 2 /A x ] n ~ [1/A c ] n . As a result, we can perform the loop expansion with quadratic divergences 
in the large N c limit, and extrapolate the results to the real QCD with N c = 3. 

In the following, we will apply the ChPT with HLS combined with the Wilsonian matching 
to finite temperature. There the expansion parameter is T/F n (A) instead of T/F n (0) used in 
the standard ChPT. Since F n (A) > F n (0), the present formalism can be applied in the higher 
temperature region than the standard ChPT. 

3.2 Two- Point Functions in Background Field Gauge 

In the present approach hadronic thermal effects are included by calculating pseudoscalar and 
vector meson loop contributions. In this section we show details of the calculation of the 
hadronic thermal corrections as well as the quantum corrections to the two-point functions in 
the background field gauge. 

3.2.1 Background field gauge 

In this subsection we briefly review the background field gauge following Ref. jHSHEUEj- We 
introduce the background fields £ L R as 

£l,r = £l,r£l,ri (3.27) 
where £l,r are the quantum fields. It is convenient to write 

£ P = exp[if a T a /F n ] , £ s = exp[ia a T a /F a ) , (3.28) 

with 7r and a being the quantum fields corresponding to the pseudoscalar NG boson ir and the 
would-be NG boson a. The background gauge field V ^ for the HLS gauge boson is introduced 

as 

V„ = + gfa (3.29) 

# 4 Wc consider that the scale at which the HLS theory breaks down is the same order as A x since the chiral 
symmetry is non-linearly realized. 
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where p M is the quantum field. It is convenient to introduce the following fields corresponding 
to dj^ and dy^ + 

A = ■ t R - o,l L ■ f L ] + \\^ R n,t R - l L c,t L ] , 



(3.30) 



The field strengths of An and V» are defined as 



A 



fJ,V 



v 



V A 

v 111 



A V 



A A 



(3.31) 



Note that both A^ and V \ v do not include any derivatives of the background fields £ R and 



V 



A 



fill 



(3.32) 



Then A^ u and V M1/ correspond to A^ u and V^j, in Eqs. f!3.25|) and ()3.26|) . respectively. 

In the background field gauge the quantum fields as well as the background fields £ RL 
transform homogeneously under the background gauge transformation, while the background 
gauge field V ^ transforms inhomogeneously: 



a 



h{x) ■ ■ h\x) + ih{x) ■ d^h\x) 
h(x) ■ £ l ,r • h)(x) , 
h(x) ■ 7r • h*(x) , 
h{x) ■ & ■ h*(x) , 
h(x) ■ p M • h\x) . 



(3.33) 



Then the expansion of the Lagrangian in terms of the quantum field manifestly keeps the HLS 
of the background field V ^ and thus the gauge invariance of the result is manifest. We fix 
the background field gauge in 't Hooft-Feynman gauge as 



£ GF = -tr + gF a af\ , (3.34) 

where D^ denotes the covariant derivative with respect to the background fields: 

Vp v = & t p v -tfsF,p v \. (3.35) 
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The FP ghost associated with C GF is given by jHS] 



C FP = 2*tr C(D^D, + (gF a ) 2 )C +■■■ 



(3.36) 



where ellipses stand for the terms including at least three quantum fields. 

Now the 0{p 2 ) Lagrangian, £( 2 ) + £gf + £fp, is expanded in terms of the quantum fields, 
7T, a, p, C and C. The terms which do not include the quantum fields are nothing but the 
original 0(p 2 ) Lagrangian with the fields replaced by the corresponding background fields. 
The terms which are of first order in the quantum fields lead to the equations of motions for 
the background fields: [SU UH] 



The quantum correction as well as the hadronic thermal corrections are calculated from the 
terms of quadratic order in the quantum fields. The explicit forms of the terms as well as the 
Feynman rules obtained from them are shown in Ref. [TB*] . 

3.2.2 Two-point functions at T = 

In this subsection, we calculate the quantum corrections from pseudoscalar and vector mesons 
to the two-point functions of A^, V M and V ^ at zero temperature. In the present analysis, as 
was done in Ref. ^3] (see also subsection I4.2J1 . we neglect possible Lorenz symmetry violating 
effects caused by the intrinsic temperature dependences of the bare parameters, and use the 
bare HLS Lagrangian with Lorenz invariance even at non-zero temperature. Then the results 
obtained in this subsection are identified with the quantum corrections from pseudoscalar and 
vector mesons at non-zero temperature by taking the intrinsic temperature dependences of 
the parameters into account. * 5 In the following, we write the two-point functions of A^-Au, 
V^-Vj,, V^-Vu and V^-V u as LT^, n^, Tly and n^, respectively. 

In the present analysis, it is important to include the quadratic divergences to obtain the 
RGEs in the Wilsonian sense. Here, following Ref. [631 IMl ITB*] . we adopt the dimensional 
regularization and identify the quadratic divergences with the presence of poles of ultraviolet 
origin at n = 2 jHEj. This can be done by the following replacement in the Feynman integrals: 



^ 5 See next subsection for dividing the quantum corrections from the hadronic thermal corrections. 



D»A» 



-z(o-l) ^-V^A"] +0( P 4 ) , 

V") = Q(p*) , 

g 2 Fl (V" - ) + Otf) . 



(3.37) 
(3.38) 
(3.39) 




(3.40) 
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Figure 3.2: Diagrams for contributions to TL^ at one-loop level. The circle (o) denotes the 
momentum-independent vertex and the dot (•) denotes the momentum-dependent vertex. 

On the other hand, the logarithmic divergence is identified with the pole at n = 4: 

i + l^lnA 2 , (3.41) 

where 

- = J?— -1e + ln(4vr) , (3.42) 
e 4 — n 

with 7s being the Euler constant. 

Let us start from the one-loop corrections to the two-point function of A p -A u denoted by 
n^. We show the diagrams for contributions to 11^ at one-loop level in Fig. 13.21 We obtain 
these contributions as * 6 

nf^(p) = -N f aM p 2 g^B ( ™ c \p;M p ,0) , (3.43) 
Ilf^ip) = N f ~B( Yac ^ v (p;M p ,0) , (3.44) 

nf^(p) = N f (a- 1)^^(0), (3.45) 
where the functions £?Q Vac ' ) , £?( vac W anc l Aq are defined by the following integrals |13j : 

(p ' M " M2) - J i(2*y\Mi-B][Nq-(k-„y] ■ (3 ' 48) 

Next we calculate the one-loop corrections to the V^-Vp two-point function denoted by 
n{| 1!/ . We show the diagrams for contributions to at one-loop level in Fig. 13.31 We obtain 



# 6 Hcre we put the superscript (vac) in the functions in the right-hand-side to distinguish them from the 
functions expressing the hadronic thermal corrections [see next subsection]. 
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7C 



71 




PP 

(a) 



(b) 



71 

(c) 



(d) 



Figure 3.3: Diagrams for contributions to TLu at one-loop level 




(d) 



(e) 



o 

(f) 




71 



7t 

(g) 



Figure 3.4: Diagrams for contributions to ITy at one-loop level 



these contributions as 



n[ a> » = -N f M p 2 g^B^ c \ P] M p ,M p ) 

nj^Go) = lN f B(™^{p;M p ,M p ), 

11 o 

n j c )^(p) = ^ 2 ~ a ^ JV / ^vac)^(p. 0>0 ) > 

1 8 



(3.49) 
(3.50) 

(3.51) 
(3.52) 



The one-loop corrections to Hy 1 , the diagrams for contributions to which are shown in 



Fig. 13 A\ are evaluated as 



n 



n 
n 



(b)fiv 



(c)/xf 



(P) 
(P) 
(P) 



_ NfB (^ {p . Mp; Mp) + 4iV / (^y - p^)Br\p; M p , M p ), (3.53) 



-N f M p 2 g^B^\ P] M p ,M p ), 
-N f B^ v {p;M p ,M p ), 



(3.54) 
(3.55) 
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Figure 3.5: Diagrams for contributions to n^,, at one-loop level. 

n!f"» = nN f g^A^(M p ), (3.56) 

4 e) ^(p) = -2A^4 vac) (M p ), (3.57) 

4 />y (p) = lN f B^^(p-M p ,M p ), (3.58) 

o 

n^(p) = ^>B<™>'»'(p;0,0) ) (3.59) 



where n denotes the dimension of the spacetime. 

We also show the one-loop diagrams for contributions to in Fig. 13.51 These are given 

by 



U^f(p) = N f M 2 p g^B^ c \p;M p ,M p ), (3.60) 

n ?!r>) = \NfB^»(p;M p ,M p ), (3.61) 

11 o 

nfhp) = a ^~ a) N f B^ ac ^(p;0,0), (3.62) 
" 8 

41"" (P) = ^/^4 Vac) (0), (3.63) 

n^f(p) = ^/^^(A^)- (3-64) 



At tree level the two-point functions of A p , V p and V p are 

nf cc) ^(p) = < bare <r + tz^tfr - pv) 
[ (t 
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n^(p) = ^^-^-bV-PW, 

"bare 

n ^(p) = -F CT 2 jbare ^ + Z 3 ,baxe(pV l/ -^) • (3.65) 

Thus the one-loop contributions to give the quantum corrections to F 2 and z 2 . Similarly, 
each of the one- loop contributions to II j^, U v u and includes the quantum corrections to 
two parameters shown above. For distinguishing the quantum corrections to two parameters 
included in the two-point function, it is convenient to decompose each two-point function as 

W v {p) = II s \p)g^ + Tv LT {p){g^p 2 - p»p v ). (3.66) 

Using the above decomposition, we identify n^ 1 " loop ' ) (p 2 ) with the quantum correction to F 2 , 
IL I f r( - 1Aoop \p 2 ) with that to z 2 , and so on. It should be noticed that the following relation is 
satisfied [531 E3|: 

n^(p 2 ) = nf(p 2 ) = -n^( P 2 ) . (3.67) 

Then the quantum correction to F 2 can be extracted from any of ITj| U/ , Uy and Ryu- The 
divergent contributions coming from the diagrams shown in Figs. l3.2ll3.5l are listed in Refs. [H3J 
ITB*] . We list them in Appendix for convenience. 

3.2.3 Two-point functions at T ^ 

Let us consider the loop corrections to the two-point functions at non-zero temperature. For 
this purpose it is convenient to introduce the following functions: 

A)(M 2 ;T) 
B {p Q ,P;M 1 ,M 2 ;T) 
B^(p Q ,P;M l ,M 2 ;T) 

where p = \p\ and the Oth component of the loop momentum is taken as k° = i2mrT, while 
that of the external momentum is taken as p° = i2n'nT [n,n': integer]. Using the standard 
formula (see, e.g., Ref. [HZ]), these functions are divided into two parts as 

A {M;T) = A { ™ c) {M)+A {M;T) , (3.71) 
B (p ,p;M u M 2 ;T) = B^ip; M u M 2 ) + B (p ,p; M u M 2 ;T) , (3.72) 
B^{p ,p;M u M 2 ;T) = B^ u {p; M h M 2 ) + B^^p; M u M 2 ;T) , (3.73) 



T E 

n=—c 
oo 

T E 

n= — o 
oo 

T E 



d 3 k 



(2tt) 3 M 2 - k 2 ' 
d 3 k 1 



(2tt) 3 [M 2 - k 2 ][M 2 — (k — p) 2 } 

d 3 k {2k - pf {2k - p) u 
(2tt) 3 [M 2 - k 2 ][M 2 - {k - p) 2 } 



(3.68) 
(3.69) 
(3.70) 
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where < ac) , Pj vac) and express the quantum corrections given in Eqs. (j3.46j) - (J3.48l) . 

and Aq, Bo and B^ v the hadronic thermal corrections. We summarize the explicit forms of the 
functions Aq, B and B^ u in various limits relevant to the present analysis in Appendix El Note 
that the Oth component of the momentum p in the right-hand-sides of the above expressions 
is analytically continued to the Minkowsky variable: po is understood as po + ie (e — > +0) for 
the retarded function and po — ie for the advanced function. Then, Aq™^, Bq™^ and f?( vac )' u ' 
have no explicit temperature dependence, while they have intrinsic temperature dependence 
which is introduced through the Wilsonian matching as we will see later. 

Now, the contributions to the two-point functions at non-zero temperature are obtained by 
replacing, in the expressions listed in the previous subsection, the functions A { ™ ] , B^ and 
S (vac)^ defined in Eqs. wi th A , B Q and B^ u defined in Eqs. (ET771j> . 

At non-zero temperature there are four independent polarization tensors, which we choose 
as defined in Appendix |X] We decompose the two-point functions IF_^, liy 1 n^jj and Yl^ u as 

= uvn* + (g^ - u v)ir + pf n L + pfu T , (3.74) 

where = (1,0) and P£ v and P£ v are defined in Eq. ([A.ljl . Similarly, we decompose the 
function B^ v as 

B^ = uVB* + (g^ - u^u u )B s + PfB L + P£ U B T . (3.75) 

Furthermore, similarly to the division of the functions into one part for expressing the quantum 
correction and another for the hadronic thermal correction done in Eqs. (I3.71|) - (J3.73)) . we divide 
the two-point functions into two parts as 

n^(po,p; T) = U^ c ^( Po ,p) + W v (p Q ,p; T) . (3.76) 

Accordingly each of four components of the two-point functions defined in Eqs. f!3.74j) is divided 
into the one part for the quantum correction and another for the hadronic thermal correction. 
We note that all the divergences are included in zero temperature part IF va ^ . 

From the contributions in Eqs. ()3.43J) - (j3.45J) . the temperature dependent parts are obtained 

as 

-N f aM 2 p B (p ,p; M p , 0; T) + N f ^B\p ,p; M p , 0; T) 
+ N f (a-1)A (0;T), (3.77) 
-N f aM*B (p ,p; M p , 0; T) + N f ^B s (p ,p; M p , 0; T) 
+ N f {a-l)Ao(0;T), (3.78) 
N f ^B L (p ,p;M p ,0;T), (3.79) 

N f ^B T (p o ,p;M p ,0;T). (3.80) 



ni(po,p;T) = 
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As we will see later in subsection 13.4.11 we define the vector meson mass as the pole of 
longitudinal or transverse component of the vector meson propagator. From the contributions 
in Eqs. (|3.53]) - (|3.59jl . we obtain the hadronic thermal corrections to U. 1 ^ as 

T? v (p ,p;T) = 2N f A (M p ;T)-M 2 p N f B (p ,p;M p ,M p ;T) 

+ jN f B\p ,p;0,0;T) + ~N f B\p o ,p;M p ,M p ;T) , (3.81) 

n s v (p ,p;T) = 2N f MM P \T)-M 2 p N f B Q {p^p-M p ,M p -T) 
a 2 „^s, 9 ^ 



N f B s (p ,p;0,0;T) + -N f B s (p ,p;M p ,M p ,T) , (3.82) 



U v (p ,p;T) = -4p 2 N f B ( P(h p;M p ,M p ,T) 

+ jN f B L (p ,p; 0, 0; T) + ^N f B~ L (p ,p; M p , M p ; T), (3.83) 
%{p^p-T) = -Ap 2 N f B (p ,p;M p ,M p ;T) 

+ yJV/B T (po, p; 0, 0; T) + ^N f B T (p , p; M p , M p ; T). (3.84) 

Another two-point function associated with the vector current correlator is n^. From the 
corrections in Eqs. I|3.6()j) - (|3.fi4j) . we get the temperature dependent parts as 



lf v]] (p ,p;T) = -N f A (0;T) + -N f A (M p ,T) + N f M 2 B (p ,p;M p ,M p ,T) 



^N f A o (0;T) + ± 

+ ]-N f B\p ,p;M p ,M p ,T) + ° (2 ~ a) N f B\p ,p; 0, 0; T), (3.85) 
o o 

n^||(po,p;T) = ^ / Ao(0;T) + ^ / A (M p ;T) + ^ / M p 2 5o(po,p;M p ,M p ;T) 

+ ^/F(p , P! M p , M p ; T) + a(2 ~ a) iV / F(po, P! 0, 0; T), (3.86) 

n^(po,p;r) = ±N f B L (p Q ,p; M p , M p ; T) + ° (2 ~ NfB~ L (po,p; 0, 0; T% (3.87) 

nL(p ,p;T) = liV^ r (p , p; M p , M p ; T) + ° (2 ~ a) iV^ T (p , p; 0, 0; T). (3.88) 

o o 

In subsection 13.4.21 we will define the parameter a from the direct / -fim coupling using the 
two-point function H^ u . From the corrections in Eqs. (j3.49|) - (j3.52|) . we get the temperature 
dependent parts as 

Tl\(p ,p;T) = -(a-l)N f M0;T)-N f M 2 B (p ,p;M p ,M p ;T) 



,2 



+ ^N f B\p ,p;M p ,M p ,T) + ^-^N f B\p ,p;0,0;T) , (3.89) 
nJ(p ,p;T) = -(a-l)N f A (0-,T)-N f M 2 B {p ,p-,M p ,M p] T) 

+ ^N f B s (p ,p;M p ,M p ;T)+ ( - 2 ~ 8 a > NfB~ s (p Q ,p; 0, 0; T) , (3.90) 
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n[(po,p;T) = lN f B L {p ,p-M p ,M p -T) + ^-^N f B L {p ,p-OAT), (3.91) 
IIJWjT) = ±N f B T (p ,p;M p ,M p ;T) + ^^N f B T (p ,p; 0,0; T). (3.92) 

At the end of this subsection, we note that, using the relations shown in Eq. (jB.14|) . we 
obtain 

% = -Kw = sj 

= n^ = -n; ll = nj 

= -N f - [A (M P ; T) + A (0; T)} - N f M^B Q (p Q ,p; M p , M p] T) . (3.93) 



Since the quantum corrections to the corresponding components are identical as we have shown 
in Eq. (|3.ti7|) . the above relation implies that the components Hy, n^y, ITy, Yl v , H v ^ and II| 
agree: 

u v = -u vll = nj =u s v = -u s vll = njj ■ ( 3 - 94 ) 

This relation is important to prove the conservation of the vector current correlator as we will 
show in subsection 13.3.21 * 7 

3.3 Current Correlators 

In this section, following Ref . ^3] , we review how to construct the axial- vector and vector cur- 
rent correlators from the two-point functions calculated in the previous section. The correlators 
are defined by 

G%(po = vJn,P\T)6 ab = dr J d 3 xe^ + ^ ^Jg a (r,x)J» b (0,0)) ^ , 
G v v (p = iu n} p;T)5 ab = dr j d 3 xe~^ + ^ (j£(r, 5)^(0, 0)}^ , (3.95) 



where J^ a and J£ are, respectively, the axial-vector and vector currents, u n = 2nnT is the 
Matsubara frequency, (a, b) = l,...,Nj — 1 denotes the flavor index and ( )p the thermal 
average. The correlators for Minkowski momentum are obtained by the analytic continuation 
of p . 

# 7 Actually, for the conservation of the vector current correlator, H v = ~riy|| = Ily and H v = — Hy\\ = njj 
are enough, and Uy and Hy can be generally different. 
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3.3.1 Axial-vector current correlator 

Let us consider the axial- vector current correlator G A V . For constructing it, we first parameter- 
ize the background field £r and £l- It is convenient to take the unitary gauge of the background 
HLS. Then, £ R and £l are parameterized 

£ L = e-*, £ R = e^, = a T a , (3.96) 

where <j) is the interpolating background field corresponding to the pion field. In terms of <p, 
the background A^ is expanded as 

A li = A lt + id li $ + --- , (3.97) 

where ellipses denote the terms including two or more fields and the axial-vector external gauge 
field A^ is defined as 

A, = \ (TZ, - £„) . (3.98) 

There are two contributions to the axial- vector current correlator G A U : one from the 0- 
exchange diagram and another from the one-particle-irreducible (1PI) diagram of A^-Ay in- 
teraction. By adding these two contributions, G A V is expressed as 

G7 = W ^ + nf • (3.99) 
Substituting the decomposition of IT^ in Eq. ()5.17j) into the above , we obtain 

G% = Pl v G\ + PPG%, (3.100) 

where 

ql = p 2 njnj 

-[pgn* ± -^ni] +ii± ' 

g t a = -ni + n^. (3.101) 

From this form we can easily see that the current conservation p^G^ = is satisfied, and the 
pion couples only to the longitudinal part G L A . 

3.3.2 Vector current correlator 

In the background field gauge the vector meson inverse propagator is related to the two-point 
function of the background field V M , Tly as 

Hp- 1 ) 1 " = IT7 = u'VIlV + (g^ - u^Uy + PfTly + piflfy . (3.102) 
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It is convenient to decompose the full propagator in a same way: 



- iD^ = u»u u D v + (g^ - u»u u )D s v + P£ V D V + P$ V D 



T 

V ■ 



(3.103) 



After some algebra, these four components are expressed as 



D L 



D 



n? 



P 2 u v (u v - u v ) - p 2 u v (u v - n* v ) ' 



P 2 (u^ - n* v ) 



P 2 w v (u^ - u s v ) - p 2 n v (n v - nv) : 

2-nL 



p z U 



P gnV(n^ - n^) - p 2 n^(n£ - n' v ) ' 



D s v + 



T\T _ Tfs 



(3.104) 



where we define p — \p\. By using the above vector meson propagator —iD^ u and the two-point 
functions ITj^ and the vector current correlator Gy, which consists of the contributions 
from the vector meson exchange and 1PI diagrams, is obtained as 

G v v = n^A^nJii + n f • ( 3 - 105 ) 

Substituting Eq. (13.1031) together with Eqs. (I3~T1) into Eq. (I3.1()5|) . we obtain G y v as 



G v u = u»u u G v 



IIV 



u»u y )G s v + P^G V + P^G 



T 
\ ■ 



(3.106) 



where each component is expressed as 



G\ 



Gy 



^U^UyUy^ 



v 



-{-p 2 u vll (u v - n£) +p 2 (u vll Wy - u vll u v )} 



+ nf. 



SIS V 

v ~ Tl 17 



^{n^-n^,,} 



p 



n? 



v 



, {-^(n^n - u vll u v ) + p 2 u s vll (u v - u s v )} 



v 



L n2 



+ nf, 



n 



+ 



v 
(n 



-pnviiWn^-n^n*,,,} 

2 ^{-Pl^ylly^ - Ily^Uy) + P 2 W VII (U V - Tly)} 



S 

ni 



+ n 



(3.107) 
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One might worry that the above form does not satisfy the current conservation p^G 
However, since the following equalities are satisfied [see Eq. (|3.94[) ]. 



V 



n 



-n 



VI! 



1L v ~ v\\ 



II||: 

nf. 



(3.108) 



each component of the above G v is rewritten as 



Gy 
Gy 



Gy 



Uy 



0. 



U v Uy(Uy+2U vll ) + 



p 2 Uy-p 2 Ut 
P 2 



: n ^ii) 2 



n^(n^ + 2n^) + (n^) 2 



+ n . 



n 



(3.109) 



Now we can easily see that the current conservation p^Gy = is satisfied since p^P 1 ^ = 
p^P^ u = 0. In the present analysis we further have Tl v = Yl v as can be seen from Eq. ()3.94|) . 



Finally we obtain the following G 



/IV _ 



(jy 



(3.110) 



where 



Gy 



Gy 



Uy(Uy + 2Hl 

Ky - Uy 

nf,(n£ + 2ir 



\ 'II 



+ nf, 
+ n[. 



(3.111) 



Here we have dropped the terms (riyn) 2 and (n^,,) 2 since they are of higher order. 



3.4 Thermal Properties of n-p Parameters 

In this section, we briefly summarize the temperature dependences of several parameters of 
pions and vector mesons in hot medium following Refs. EE] • 



3.4.1 Vector meson mass 

In this subsection, we first define the vector meson pole masses of both the longitudinal and 
transverse modes at non-zero temperature from the vector current correlator in the background 
field gauge and show the explicit forms of the hadronic thermal corrections from the vector 
and pseudoscalar meson loop. 
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Let us define pole masses of longitudinal and transverse modes of the vector meson from 
the poles of longitudinal and transverse components of the vector current correlator in the rest 
frame, which are given by Eqs. (j3.110|) and (j3.111j) . Then, the pole masses are obtained as the 
solutions of the following on-shell conditions: 



= Re [U s v (po = m L p (T), 0; T) - U v (p = m L p (T), 0; T 
= Re [ILMpo = <(T), 0; T) - n£(p = mJ(T), 0; T 



(3.112) 



where m p (T) and m p (T) denote the pole masses of the longitudinal and transverse modes, 
respectively. As we have calculated in section l3~2"| Tly(p ,p]T), Tl v (p ,p;T) and Yly(p ,p;T) 
in the HLS at one-loop level are expressed as 



U v ( Po ,p-T) 
Il v (p ,P;T) 

Uy(p 0l p;T) 



F 2 (M p ) + U s v (p 2 ) + Uy{p ,p;T) 

-.2 



P 



U v T (p 2 ) + U v (p ,P;T) 



9 2 {M p 
P 2 



+ U^ T (p 2 ) + Uy(p ,p;T) 



(3.113) 



9 2 (M P ) 

where the explicit forms of the finite renormalization effects TL v (p 2 ) and U v T (p 2 ) are given in 



Eqs. (IC31J) and ()C.32j) . and those of the hadronic thermal effects H v (po,p;T), U v (p ,p;T) 
and Tly(p ,p;T) are given in Eqs. (JH331), (|3~£3jl and flUBD- Substituting Eq. (j3.113jl into 
Eq. (I3.112jl . we obtain 



M 2 + g 2 (M p 



Ren^(^) + Ren y (p ,0;T) 
-Reff^(p2) -ReTly{po,0;T) 



Jp =m£(T) 



[m T p {T 



M 2 + g 2 (M p ) 



Ren^(^) + Ren;(p ,0;T) 



Ren^(p2)-R e n y (p ,0;T) 



Jpo=mJ(T) 



(3.114) 



We can replace m p (T) and mJ(T) with M p in the hadronic thermal effects as well as in 
the finite renormalization effect, since the difference is of higher order. Then, noting that 
Refl£(p 2 = M 2 ) = Refl v T {p 2 = M 2 ) = as shown in Eq. (fHT9jl . we obtain 



[m L p {T 
[m T p (T 



p ■ 

r2 „2/ 



M 2 - g 2 (M p ) Re ]ll v (M p , 0; T) - U v (M p ; 0; T) 
M 2 - g\M p ) Re \%{M p ; 0; T) - %{M p] 0; T) 



(3.115) 
(3.116) 
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s ^ j> 

Let us calculate the explicit form of the pole mass using the expression of U v , YL V and YL V 
obtained in section 13.21 Here we note that Eqs. (jB.24|) and (|B.25|) imply that B L — B s agrees 

e £j a —T S 

with B — B in the rest frame. Then in the rest frame Tl v — YL V agrees with YL V — YL V . Thus 
the longitudinal pole mass is the same as the transverse one: 

m L p (T)=m T p (T)=m p (T) . (3.117) 

By using the low momentum limits of the functions shown in Eqs. (jB.23)) and (jB.24)) . YL V — lf v 
in the rest frame is expressed as 



Tly(p 0j p = 0;T) -Ily{p ,p = 0;T) 



N f 



■ M l -■>, , 3 ~ A 



r ^G 2 (p ; T) - Jl{M p] T) + C~f- rf)i^(po; M p , T) + ^F*(p o] M p ; T) (3.118) 



where the functions /„, J^, and are defined in Appendix |Dj Substituting the above 
expression into Eq. (|3.115jl and using the relation 

-^Fi(M p ;M p ;T) = \p 2 (M p , M p , T) - JL jf( Mp , T) , (3.119) 

we find that the vector meson pole mass is expressed as 

a 2 * , . , _ 4 ~o , , , 33 



m 2 (T) = M p 2 + N f g 2 



--G 2 (M p , T) + -J 2 (M p , T) + -M p 2 F 3 2 (M p ; M p ; T) 



(3.120) 



The contribution in this expression agrees with the result in Ref. [Hj which is derived from the 
hadronic thermal correction calculated in the Landau gauge in Ref. [HE] by taking the on-shell 
condition (|3.115|) . #8 

We study the behavior of the pole mass in the low temperature region, T <C M p . In this 
region the functions and are suppressed by e~ Mp ^ T , and thus give negligible contribu- 
tions. Since G 2 ~ ~T5jpzi the vector meson pole mass increases as T 4 in the low temperature 
region, dominated by the 7r-loop effects: 

m 2 (T)^M 2 + ^^T A for T <C M p . (3.121) 
p y ' p 360F W 2 p y ' 

Note that the lack of T 2 -term in the above expression is consistent with the result by the 
current algebra analysis 



# 8 The functions used in this paper are related to the ones used in Ref. [Tl| as Jf = ^iJ'i- G2 = ^1 G2, and 
so on. 
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3.4.2 Pion parameters 

First we study the on-shell structure of the pion. For this, we look at the pole of the longitudinal 
component G L A in Eq. (j3.1()lj) . Since both IT^ and IT^ have imaginary parts, we choose to 
determine the pion energy E from the real part by solving the dispersion formula 

= [p 2 ReU t ± (p ,p;T)-p 2 ReW ± (p Q ,p ] T)] po=E , (3.122) 

where p = \p\. As remarked in section l3~2] in HLS at one-loop level, rF ± (po,p;^) and 
nj_G°cbP; T) are of the form 

ni(p ,p;T) = F2(0) + ni(p 2 ) + nl(p ,p;T) , 

ni(p ,p;T) = F2(0) + ni(p 2 ) + ni(p ,p;T) , (3.123) 

where IT^ (p 2 ) is the finite renormalization contribution, and rP_|_(p ,P; T) and n^_(p ,p;T) are 
the hadronic thermal contributions. Substituting Eq. ()3.123J) into Eq. (J3.122J) . we obtain 

= (E 2 — p 2 ) \f 2 (0) + Re fl s ± (p 2 = E 2 — p 2 ) 

+ E 2 Ren' ± (E,p;r) -p 2 ReW L {E,p-T) . (3.124) 

The pion velocity v n (p) = E/p is then obtained by solving 

_ F 2 (0) + Reni(p,p;T) 

ffCp) " i?(o) + Heni(p > p;r) " (3 - 125) 

Here we replaced -E by p in the hadronic thermal terms Il i L (£', p) and fij^i?, p) as well as in the 
finite renormalization contribution IFj (p 2 = E 2 — p 2 ), since the difference is of higher order. 
[Note that U s ± (p 2 = 0) = 0.] 

Next we determine the wave function renormalization of the pion field, which relates the 
background field to the pion field n in the momentum space as 

= 7f/F(p;T). (3.126) 

We follow the analysis in Ref. [70] to obtain 

F 2 (p;T) = Ren t ± (E,p;T) = F 2 (0) + Ren t ± (p,p;T) . (3.127) 

Using this wave function renormalization and the velocity in Eq. (j3.125jl . we can rewrite the 
longitudinal part of the axial- vector current correlator as 

r ,L, p 2 n^(po,p;r)ni(po,p;T)/F 2 (p ; r) L 

G> °' P) = -H-vUP)P> + n,( Pa ,p;T)] + n xfa.P^)- (3-128) 
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where the pion self energy U. n (p ,p; T) is given by 
U n (p ,p;T) 



ReW ± (E,p;T) 

pi { III (p , p; T) - Re ni ( E, p; T) } - p 2 {III (p , p; T) - Re ni ( E , p; T) } 



(3.129) 



Let us now define the pion decay constant. A natural procedure is to define the pion decay 
constant from the pole residue of the axial- vector current correlator. From Eq. (|3. 128(1 . the 
pion decay constant is given by 

ni(£,p; T)U S ± (E, p;T) 



F2(p;T) 

[F|(0) + fli(p,p; T)} [F 2 (0) + ni(p,p; T)} 



(3.130) 



F 2 (p;T) 

We now address how / 2 (p;T) is related to the temporal and spatial components of the pion 
decay constant introduced in Ref. [2D]- Following their notation, let /* denote the decay 
constant associated with the temporal component of the axial-vector current and the one 
with the spatial component. In the present analysis, they can be read off from the coupling of 
the 7f field to the axial- vector external field A^. 

Ui(E,p;T) F 2 (0) + ni (p,p;T) 



f s MT) 



F(p;T) F(p;T) 
U s ± (E,p;T) F 2 (0) + tll(p,p;T) 



(3.131) 
(3.132) 



F(p;T) F(p;T) 
Comparing Eqs. (I3~T3TJ) and (EHS^ with Eqs. 1ETT25I) . (I3~T27() and (jUSDD , we have [2011101 

F(p;T) = Re/*(p;T), (3.133) 
tf(p;T) = /*(p;T)/;(p;T) . (3.134) 

By using the above decay constants, the pion velocity at one-loop level is expressed as [201 EG] 

#9 



F 2 



Re/^p;T)-Re/*(p;T) 



(3.135) 



Substituting Eqs. (|B.1() and (|B.2|) into Eqs. ((3.77J) and ()3.78|) . we obtain n_|_ and U ± for 
the on-shell pion as 

TT ± (p =p + te J p;T)=N f (a-l)I 2 (T) 



# 9 This form of the pion velocity looks slightly different from Eq. (|3.125|) . However, Eq. (|3.135|l is actually 
equivalent to Eq. (|3.125f) at one-loop order, and more convenient to study the temperature dependence of the 
pion velocity in the low temperature region. 
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N f 

H -a 

2 



5 (p + te,p; M pi O;T) - 2M 2 B (p + te 7 p; M p ,0;T) 



K(p =p+ie,P;T)=N f (a-l)I 2 (T) 



N f 

H 

2 



-B s (p + *e,p; M p , 0; T) - 2M?B (p + *e,p; M p , 0; T) 



(3.136) 



(3.137) 



where we put e — > +0 to make the analytic continuation for the frequency po = i^nnT to the 



Minkowski variable. We show the explicit forms of the functions B , B and So m Eqs. f)B.3|) . 
(JEU) and (EH). 

In general the pion velocity in medium does not agree with the value at T = due to the 
interaction with the heat bath. Below T c , since the temporal component does not agree with 
the spatial one due to the thermal vector meson effect, II ± ^ Il^_, the pion velocity v n (p;T) 
is not the speed of light. As we will see below, in the framework of HLS the pion velocity 
receives a change from the p-loop effect for < T < T c . When we take the low temperature 
limit (T <C M p ) and the soft pion limit (p <C M p and p < T), the real parts of IT^ and Yf ± are 
approximated as 



Rell^Qoo = p + ie,p;T) 
ReIl^(po = p + ie,p;T) 



N f I 2 (T) + N 



a ~ 



Nf a 



Nf I 2 (T)-N f —h(T) + N f a 
p 



M n 



M n 



e - A/p/T T 3/2 (3.138) 

e -M p /T T 3/2( 3 139) 



By using Eqs. (I3.138J) and ()3.139|) and neglecting the terms proportional to the suppression 
factor e~ Mp / T , the pion velocity is expressed as 

2a 



N 



F 2 M 2 

TV p 



h{T) 



N f 7 T 4 

1 — CL7X 



15 



F 2 M 2 

n p 



< 1. 



(3.140) 



This shows that the pion velocity is smaller than the speed of light already at one-loop level 
in the HLS due to the p-loop effect. This is different from the result obtained in the ordinary 
ChPT including only the pion at one- loop [see for example, Ref. j2D] and references therein]. 
Generally, the longitudinal p contribution to Tl ± expressed by B l in Eq. ()3.136p differs from 
that to by B s in Eq. ()3.137|) . which implies that, below the critical temperature, there 
always exists a deviation of the pion velocity from the speed of light due to the longitudinal 
p-loop effect: 

vl(p; T) < 1 for < T < T c . (3.141) 

Next, we study the temporal and spatial pion decay constants in hot matter defined by 
Eqs. ()3. 131)1 and ()3.132j) . The imaginary parts of 11^ and lfj_ in the low temperature region 
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are given by 



t S<CT N f t 

lmU ± (p =p + te,p;M p ,0;T) ~ -f- a ImB (p = p + ie,p; M p , 0; T) 



T<M P jV> 

8tt 



MV Mp/T , (3.142) 



Imn ± (p =p + 2e,p;M p ,0;T) ~ a ImB (po = p + ie, p; M p , 0; T) 

~ --^aMe- M ' /T (3.143) 

— t s 

Thus the contributions from the imaginary parts ImlT^ are small because of the suppression 
factor e" M "/ T . From Eqs. lCTT33|l and (gUSj with J 3 (T) = T 2 /12 and J 4 (T) = vr 2 T 4 /30, we 
obtain the following results for the real parts of /* and 



[Re/*]F ~ F 2 -iV 



(X 7T „, / M, 



12 30M 2 V 8tt 



T 4 + aW^e- M "/ T r 3 / 2 



[Re/*]F ~ F 2 -N 



12 + 30Mf a V8vr3 e 



(3.144) 



We note that F = Re/£ as shown in Eq. ()3.133|) . Then, neglecting the terms suppressed by 
e~ Mp ' T , we obtain the difference between (/*) 2 and as 

A/" 

(£) 2 -££ ^ if a7r2 M| >0 - ^ 145 ) 

This implies that the contribution from the p-loop (the second and third terms in the brackets) 
generates a difference between the temporal and spatial pion decay constants in the low tem- 
perature region. Similarly, at general temperature below T c , there exists a difference between 
f*F and f%F due to the p-loop effects: [Re/*]F > [Re/*]F. Since we can always take F to be 
positive, we find that Re/* is larger than Re/*: 

Re/*(p;T)>Re/*(p;T) for < T < T c . (3.146) 

This result is consistent with Eq. (|3.141j) because v 2 — 1 = (F/F 2 )[Re/* — Re/*] by definition. 
The difference between (f*) 2 and f*f* shown in Eq. (|3.145j) is tiny, and the hadronic thermal 
corrections to them are dominated by the first term (T 2 /12) in the bracket in Eq. (|3.144|) . 
Then, in the very low temperature region, the above expressions are further reduced to 

fl = (£)(£) ~ Ul? ~ F 2 - ^T 2 , (3.147) 
which is consistent with the result obtained in Ref. |71j . 
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Figure 3.6: Leading contributions to the ^nir interaction. 
3.4.3 Vector meson dominance 

As shown in section 13. 1| the HLS theory can reproduce the vector dominance (VD) of the 
electromagnetic form factor of pion, which is phenomenologically successful in mesonic system. 
In this subsection, we study the validity of the VD in hot matter. 

We first study the direct '-firir interaction at zero temperature. We expand the Lagrangian 
()3.13j) in terms of the ir field with taking the unitary gauge of the HLS (a = 0) to obtain 

£ ( 2) = tr [d^ir] + ag 2 F* tr [p^] + 2i {^ag^j tr [ff [d^ir , tt]] 

- 2 [agFl) tr [^V] + 2z (l - |) tr [V [d,ir , tt]] + • • ■ , (3.148) 
where the vector meson field p M is introduced by 

V, = gp,, (3.149) 
and vector external gauge field V M is defined as 

V„ = ^ + £„) . (3.150) 

At the leading order of the derivative expansion in the HLS, there are two contributions shown 
in Fig. l3.6[ i.e., the direct '-firir interaction [Fig. 13.61 (a)] and the p-mediated interaction [Fig. 13.61 
(b)]. The form of the direct '-firir interaction is easily read from Eq. ()3.148j) as 

^% M) = <9-kr{l-^), (3.151) 

where e is the electromagnetic coupling constant and q and k denote outgoing momenta of the 
pions. This shows that, for the parameter choice a = 2, the direct ^yim coupling vanishes (the 
VD of the electromagnetic form factor of the pion). 

At the next order there exist quantum corrections. In the background field gauge adopted 
in the present analysis the background fields and V M include the photon field A^ and the 
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background pion field tt as 



All 



V, 



eQA, 



i e 



A* [Q , 7f ] 



[9 M 7T , 7r] + • • 



(3.152) 



2i*(0) 

where -F„-(0) is the on-shell pion decay constant (residue of the pion pole) in order to identify 
the field tt with the on-shell pion field, and the charge matrix Q is given by 



Q 



( 



\ 



2/3 



-1/3 



\ 



■1/3 J 



(3.153) 



The direct 77T7r interaction including the next order correction is determined from the two- 
point functions of V^-V u and A^-A u and three-point function of V^-Aa-Ap. We can easily show 
that contribution from the three-point function vanishes in the low energy limit as follows: Let 
T^^j denote the V ^-A a -Ap three-point function. Then the direct ^irix coupling derived from 

Since the legs a and (3 of T^^- are 
carried by q or k, qJ^pT^-^ is generally proportional to two of q 2 , k 2 and q ■ k. Since the loop 
integral does not generate any massless poles, this implies that QakpT^j vanishes in the low 
energy limit q 2 = k 2 = q ■ k = 0. Thus, the direct ^tttt interaction in the low energy limit can 
be read from the two-point functions of V p-V v and A^-A u as 



VAA 

this three-point function is proportional to q a kpT^j. 



"/WIT 



(3.154) 



where p p 



+ k) u is the photon momentum. Substituting the decomposition of the two-point 

- 0, we obtain 



function given in Eq. ()3.66|) and taking the low-energy limit q 2 = k 2 

r(vac)5^ 2 = 0) - 



p- 



77T7T 



e(q — ky 



2 



(3.155) 



where we used n^ ac (q 2 = 0) = F 2 (0). Comparing the above expression with the one in 
Eq. ([3.151)1 . we define the parameter a(0) at one-loop level as 

n[, vac) V = o) 



a(0) 



^ 2 (0) 



(3.156) 



We note that, in Ref. [21], a(0) is defined by the ratio F 2 (M p )/F 2 (0) by neglecting the finite 
renormalization effect which depends on the details of the renormalization condition. While 
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the above a(0) in Eq. ()3.156|) denned from the direct '-firir interaction is equivalent to the one 
used in Ref. ^3]. In the present renormalization condition ()C.16|) . n|| Vac ' ) ' s (p 2 = 0) is given by 

nj vac) V = 0) = F 2 a (M p ) + ~^yM 2 p {2 - V3 tan" 1 V3) . (3.157) 
By adding this, the parameter a(0) is expressed as 

a(0) = ^rm + w~ 2 ^k {2 ~ >/5taa ~ 1 ^ ' (3 ' 158) 

Using M p = 771.1 MeV, F w (p = 0) = 86.4 MeV estimated in the chiral limit [Ml HOI EH] #1 ° 
and F 2 (M p ) /F 2 (0) = 2.03 obtained through the Wilsonian matching for A QC d = 400 MeV and 
the matching scale A = 1.1 GeV in Ref. [T3] . we estimate the value of a(0) at zero temperature 

as 

a(0)~2.31. (3.159) 

This implies that the VD is well satisfied at T = even though the value of the parameter a 
at the scale M p is close to one |72*j . 

Now, let us study the direct ^tttt coupling in hot matter. In general, the electric mode and 
the magnetic mode of the photon couple to the pions differently in hot matter, so that there 
are two parameters as extensions of the parameter a. Similarly to the one obtained at T = in 
Eq. (|3.154|l . in the low energy limit the direct '-firir interaction derived from A^-Av and Vp-V u 
two-point functions is expressed as 



rWp;g,*) 

l 



F(q;T)F(kT) 



q u n^(g , q\ T) - k v IP^ (k , k; T) 
-l( q -k) v U^(po,P;T) 



(3.160) 



where q and k denote outgoing momenta of the pions and p = (q + k) is the photon momentum. 
F is the wave function renormalization of the background 7f field given in Eq. ()3.133|) . Note 
that each pion is on its mass shell, so that go = v n{q)q and ko = v n (k)k. To define extensions 
of the parameter a, we consider the soft limit of the photon: po — > and p — > 0. * X1 Then the 
pion momenta become 

g = -ko , q = -k . (3.161) 
# 10 In Ref. PS], it was assumed that the scale dependent F^fj) agrees with the scale-independent parameter 
Fir in the ChPT at fi = to obtain F, T (fj, = 0) = 86.4 ± 9.7 MeV. Here, we simply use this value to get a rough 
estimate of the value of the parameter a(0) as done in Ref. 

# n Note that we take the po — > limit first and then take the p — > limit for definitencss. This is natural 
since the form factor in the vacuum is defined for space-like momentum of the photon. 
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Note that while only two components and appear in q^H 1 ^ or k u U^, (q — k) u H^ 1/ 
includes all four components n|, ITjj, Il| and ILj. Since the tree part of and IlJ is —2z2p 2 
which vanishes at p 2 = 0, it is natural to use only n| and ITjj to define the extensions of the 
parameter a. By including these two parts only, the temporal and the spatial components of 
are given by 



1^(0; g,-g) 
TL T (0;g,-g) 



2<?o 



77T7T \ 



F\q-T) 

F\q-T) 



Thus we define a*(T) and a s (T) as 



a\q-T) 



T) 



nl(go,g;T)--n| l (o,0;r) 

ni(g ,g;T)-^(0,0;T) 



nj^OjT) 

K(Qo,q;T) 
nj!(0,0;T) 



(3.162) 



ni(g ,g;T) ' 

Here we should stress again that the pion momentum q p is on mass-shell: go — v n (q)q. 
In the HLS at one-loop level the above a'(g; T) and a s (q; T) are expressed as 



(3.163) 



a\q-T) 
a s (q;T) 



o(0) 
a(0) 



n || (0,0;T)-a(0)n ± (g,g;T) 

a(0)F2(0;T) 
n;(0,0;T)-a(0)nl(g,g;T) 



(3.164) 
(3.165) 



o(o)if?(o ; r) 

where a(0) is defined in Eq. (|3.156() . From Eq. ()3.93() together with Eq. f)B.9|) we obtain II 
and TTJ in Eqs. (I3.164jl and | 



as 



111.(0,0; T) 



nJ(0,0;T) 



a 2 I 2 (T) - Jf(M p , T) + 2J_ 1 (M p ; T) 



(3.166) 



Let us study the temperature dependence of the parameters a*(g; T) and a s (g; T) in the low 
temperature region. At low temperature T <C M p the functions Jf(M p ; T) and J° 1 (M P ; T) are 
suppressed by e~ Mp l T . Then the dominant contribution is given by h(T) = T 2 /12. Combining 
this with Eq. (13.138(1 and (13.139(1 . we obtain 



o(0) 



1 + 



TV/ 
12 



4a(0)J F^-T) 



(3.167) 



where a is the parameter renormalized at the scale \x = M p , while a(0) is defined in Eq. (|3. 156(1 . 
We expect that the temperature dependences of the parameters are small in the low tempera- 
ture region, so that we use the values of parameters at T = to estimate the hadronic thermal 
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corrections to a M (T). By using F„(0) = 86.4 MeV, a(0) ~ 2.31 given in Eq. (I3.159j) and 
a(M p ) = 1.38 obtained through the Wilsonian matching for (Aqcd , A) = (0.4, 1.1) GeV and 
Nf = 3 [T3*] . a* and a s in Eq. ()3.167j) are evaluated as 

f T \ ^ 

1 + 0.066 



a 



o(0) 



T V 
50 MeV J 



(3.168) 



This implies that the parameters a* and a s increase with temperature in the low temperature 
region. However, since the correction is small, we conclude that the vector dominance is well 
satisfied in the low temperature region. 



3.5 Lorentz Non- invar iance at Bare Level 



As stressed in section 12 .3^ the Lorentz non-invariance appears in the bare HLS theory as a 
result of including the intrinsic temperature dependence. Once the temperature dependence 
of the bare parameters is determined through the matching with QCD mentioned above, the 
parameters appearing in the hadronic corrections pick up the intrinsic thermal effects through 
the RGEs. Then the HLS Lagrangian in hot and/or dense matter is generically Lorentz non- 
invariant. Its explicit form was presented in Ref. [22]- in this section, we show the HLS 
Lagrangian at leading order including the effects of Lorentz non-invariance. 
Two matrix valued variables £l(x) and £r(x) are now parameterized as * 12 



e ur{x)/F* e ^in(x)/F* 



(3.169) 



where and F* denote the temporal components of the decay constant of ir and a, respec- 
tively. The leading order Lagrangian with Lorentz non-invariance can be written as 



C 



raXu. + (9^ - u n u u) tr [a^a 



(Fl) 2 u,u u + F\F' (g 



tr 



\\ a \\ 



+ 



9l 



2g 2 T 



ti[V^V aP ] . (3.170) 



Here F£ denote the spatial pion decay constant and similarly F* for the a. The rest frame of 
the medium is specified by = (1, 0) and is the field strength of V^. gi and correspond 
in medium to the HLS gauge coupling g. The parametric 7r and a velocities are defined by 



ps/pt y2 = ps/pt 



(3.171) 



# 12 The wave function renormalization constant of the pion field is given by the temporal component of the 
pion decay constant |70| . Thus we normalize 7r and a by F* and F* respectively. 
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The axial-vector and vector current correlators at bare level are constructed in terms of 
bare parameters and are divided into the longitudinal and transverse components: 



^ Ay 



r A,V 



A,Vi 



(3.172) 



where P£ V T are the longitudinal and transverse projection operators, respectively. The axial- 
vector current correlator in the HLS around the matching scale A is well described by the 
following forms with the bare parameters |2"2*1 ITS] : 

^2 j?t ps 



Ci(HLS)(P0,P) 



7T F l 

" 7r,barc 7r, bare 



^P 2 2,bare) 



G A(HLS)(P0,P) 



"-^,bare-^7?,bare — ^ (po^bare — ^^bare) i (3.173) 

where z% bare and £^ bare are the coefficients of the higher order terms, and K- ibare is the bare 
pion velocity related to F* bare and F' barc as 



^ 2 h 

7r,barc 



7r,bare 

F* 

7r,bare 



(3.174) 



Similarly, two components of the vector current correlator have the following forms: 



P r cr,\ 



^2 zpt ps 

■,bare <r,bare 



.1 ^9L,baxe z 3,baxe) 



- 2pV bare + 0(p« 



- [Pi - y^ C P 2 - ^bare] 
Gy(HLS)(P0,P) 

o\bare cr.bare 



K - VI bare p 2 - M 2 bare ] 



x 



Pi (1 - 2g 2 Li 



bare"^3,bare 



) " VS Mlc p 2 (1 - 2gl 



bare^3,barc 



T 



~ 2 (plzt harc - p 2 < bare ) + 0(p 4 ) , (3.175) 

where 2bare and £^ 2i bare denote the coefficients of the higher order terms. In the above 
expressions, the bare vector meson mass in the rest frame, M Pjbare , is 



A/f 2 = n 2 F^ F s 

p,bare £/L,bare <r,bare <7,bare 



We define the bare parameters a bare and a barc as 



''bare 



F t x 2 
<T,barc 

F* 

7r,bare 



and the bare a and transverse p velocities as 

F 



v: 



bare 



bare 



cr,bare 



bare 



V 2 

^T.bare 



F s x 2 
<r,bare 



,bare 



,bare 
S'T.bare 



(3.176) 



(3.177) 



(3.178) 



Chapter 4 



Wilsonian Matching at Finite 
Temperature 

In order to fix full temperature dependences of physical quantities, parameters of the EFT 
should be determined through the matching to QCD. In this chapter, we perform the matching 
in the Wilsonian sense discussed in section 12.21 The bare parameters of HLS theory at zero 
temperature were originally determined in Ref . [S3] , where they matched the bare HLS theory 
to the operator product expansion (OPE). The Wilsonian matching well describes the real 
world (for details, see Ref. fSj). Applying this scheme to QCD in hot /dense matter, we obtain 
the bare parameters in terms of the OPE variables like expectation value of an operator. The 
intrinsic thermal effects of bare parameters are caused by the temperature dependences of such 
condensations at a matching scale, which are evaluated in the thermal vacuum. 

In this chapter, we briefly review the Wilsonian matching proposed at zero temperature. 
Next we extend the Wilsonian matching to the version at non-zero temperature in order 
to incorporate the intrinsic thermal effect into the bare parameters of the HLS Lagrangian 
following Refs. [TH ITHj ITT)| IT7| ITH] . There we also discuss the effect of Lorentz symmetry 
violation at bare level. 

4.1 Wilsonian Matching Conditions at T = 

The Wilsonian matching proposed in Ref. jM] is done by matching the axial- vector and vector 
current correlators derived from the HLS with those by the operator product expansion (OPE) 
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in QCD at the matching scale A. * x The axial-vector and vector current correlators in the 
OPE up until 0{1/Q 6 ) at T = are expressed as [73] 



N, 



8tt 2 V 3 

7T 



3(iV c 2 -l) a,V2 



+ 



N c IT J ~ /I 2 

—G pv G^ u ) 7T 3 96(iV c 2 - 1) (1 1 \ a s (qq) 2 



N c Q 4 N c N£ \2 3NJ Q 



(4-1) 



4 QCD) (^) 



N f 



8tt 2 V 3 

7T 



3(jV 2 -i) q.V Q 8 



n c 7T y ~ /i 2 

—G^G^) 7T 3 96(iV 6 2 - 1) /l l\a s (gg> 2 



iV c Q 4 iV c iV 2 V2 3iV c ; Q 



(4.2) 



where /x is the renormalization scale of QCD and we wrote the iV c -dependences explicitly (see, 
e.g., Ref. [13]). In the HLS the same correlators are well described by the tree contributions 
with including 0(p 4 ) terms when the momentum is around the matching scale, Q 2 ~ A 2 : 



gT S \Q 2 ) = ^~2z 2 (A), (4.3) 



4 HLS) (Q 2 ) = m 2 (A) + Q 2 t 1 ~ 2 ^>M] ~ 2*i(A) , (4.4) 

where we defined the bare p mass M p (A) as 

M p 2 (A) = ^ 2 (A)F CT 2 (A) . (4.5) 

We require that current correlators in the HLS in Eqs. ()4.3|) and (|4.4jl can be matched 
with those in QCD in Eqs. (j4.1j) and (|4.2|l . Of course, this matching cannot be made for any 
value of Q 2 , since the Q 2 -dependences of the current correlators in the HLS are completely 
different from those in the OPE: In the HLS the derivative expansion (in positive power of 
Q) is used, and the expressions for the current correlators are valid in the low energy region. 
The OPE, on the other hand, is an asymptotic expansion (in negative power of Q) , and it is 
valid in the high energy region. Since we calculate the current correlators in the HLS including 
the first non-leading order [(9(p 4 )], we expect that we can match the correlators with those in 
the OPE up until the first derivative * 2 . Then we obtain the following Wilsonian matching 

^For the validity of the expansion in the HLS the matching scale A must be smaller than the chiral symmetry 
breaking scale A x as we stressed in chapter [5] 

# 2 If there exists an overlapping area around a scale A, we can require the matching condition at that A. In 
fact, the Wilsonian matching at T — in three flavor QCD was shown to give several predictions in remarkable 
agreement with experiments |fi4l 113) . This strongly suggests that there exists such an overlapping region. 
As discussed in Ref. ^3]> A -c Ahls can be justified in the large N c limit, where Ahls denotes the scale at 
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conditions [HH CS1 #3 

g(A) 
A 2 



3 



3(iV 2 - 1) a, 

71 



2tt 2 (^G^G^) 



288tt(AT c 2 



A 4 

/- \2 



F 2 (A) A 4 [l-2^ 2 (A)z 3 (A)] 



A 2 (M/(A)+A 2 ) 2 



3 



3(iV 2 



8iV c 
288tt(A^ 2 



0:.s 



7T 



^ 2 (A) F 2 (A)[1-2^ 2 (A)4a)] 



+ 3iV c J 


A 6 


2tt 2 <^ 


GfiuG^^ 




A 4 


_ 1 > 


a s (qq) 2 



3N r 



A 6 



A 2 



M P 2 (A) + A 2 
4vr(iV 2 - 1) a s (qq) 2 



2[* 2 (A)-zi(A)] 



A 6 



(4.6) 



(4.7) 



(4.8) 



The above three equations ()4.8j) . ()4.6|) and ()4.7|) are the Wilsonian matching conditions pro- 
posed in Ref. [01]. They determine several bare parameters of the HLS without much ambi- 
guity. Especially, the first condition (|4.fij) determines the ratio F n (A)/A directly from QCD. 



4.2 Wilsonian Matching Conditions at T ^ 

Next we consider the extension of the matching conditions at T = to the analysis in hot 
matter. We present the matching conditions to determine the bare pion decay constants 
including the effect of Lorentz symmetry breaking at the bare level which is caused by the 
intrinsic thermal effect. 



Case neglecting Lorentz non-invariance 

Before going to study the matching conditions taking into account the possible Lorentz non- 
invariance, we consider a naive extension of the matching to the one in hot matter: Strictly 
speaking, inclusion of intrinsic effects generates Lorentz non-invariance in bare theory. Then 

which the HLS theory breaks down (see also section m|l . We obtain the matching conditions in N c = 3 by 
extrapolating the conditions in large N c . As we mentioned above, the success of the Wilsonian matching at 
T = with taking the matching scale as A = 1.1 GeV shows that this extrapolation is valid. 
# 3 One might think that there appear corrections from p and/or it loops in the left-hand-sides of Eqs. (|4.6|l 

and (|4.7(l . However, such corrections are of higher order in the present counting scheme, and thus we neglect 
them here at Q 2 ~ A 2 . In the low-energy scale we incorporate the loop effects into the correlators. 
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we should include the Lorentz non-scalar operators such as q , y ll D u q into the form of the current 
correlators derived from the OPE [75], which leads to a difference between the temporal and 
spatial bare pion decay constants. However, we neglect the contributions from these operators 
since they give a small correction compared with the main term 1 + — as discussed in Ref. |14j . 
This implies that the Lorentz symmetry breaking effect in the bare pion decay constant is 
small, F% hgxe ~ F^baie P-^l - 111 i&ct, we will see below that the difference between them is 
caused by an existence of a higher spin operator in the OPE side of the Wilsonian matching 
condition. Thus to a good approximation we determine the pion decay constant at non-zero 
temperature through the matching condition at zero temperature, putting possible temperature 
dependences into the gluonic and quark condensates [TH I15j : 



F 2 (A;T)_ 1 r a s , 2vr 2 (fG^G^) T 3 U08a s (qq) 



A 2 8tt 2 



1 + — + — ^ 7, '— + TT 



7T 3 A 4 27 A 6 



(4.9) 



Through this condition the temperature dependences of the quark and gluonic condensates 
determine the intrinsic temperature dependences of the bare parameter F n (A; T), which is then 
converted into those of the on-shell parameter F n (fi = 0; T) through the Wilsonian RGEs. 

Case taking into account Lorentz non-invariance 

Now we study the Wilsonian matching conditions for the bare pion decay constants without 
Lorentz invariance. From the bare Lagrangian with replacement of the parameters by the bare 
ones in Eq. (j3.170|) . the current correlator at the matching scale is constructed as Eq. (j3.178|) : 

n L I -\ _ Tr.bare 7r,bare r> L 

Lt A(HLS)W0, q) — r 2 T/2 =21 ZZ 2,bare 5 
WO K 7r,barcy J 

F l F s n 2 z L — n 2 z T 

^vr / — \ 7r,bare 7r,bare o"0 2, bare H 2, bare / A -, 

Ga(hls)(?o, q) = ~ 2 2 , (4.10) 



where Kr,bare = F£ bare / Fl bare is the bare pion velocity. To perform the matching, we regard 
G^ T as functions of — q 2 and q 2 instead of qo and q, and expand G^ T in a Taylor series around 
q = \q\ = in q 2 /(—q 2 ) as follows: 

G L A {-q\ f) = Gf\-q 2 ) + G L /\-q 2 )f + • • • , 

G^(-g 2 , g 2 ) = G T /\-q 2 ) + G?\-f)t + ■ • • ■ (4.11) 

In the following, we determine the bare pion velocity Kr.bare from and via the 

matching. 
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Expanding G^ LS ^ L in Eq. (|4.1()jl in terms of q 2 / (—q 2 ), we obtain 

rpt rps 

ry (HLS)L(0)/ 2\ _ 7r,bare Tr.barc n L /.« -, \ 

J - ^ 22; 2,bare ) l 4 - i2 J 

C 1 * PS /i y2 \ 

^y(HLS)L(l)/ 2\ 7r, bare 7r, bare V 7r,bare/ q\ 

G A [-Q ) = (-q 2 ) 2 ' 

On the other hand, the correlator G^ in the QCD sector to be given in OPE is more involved. 
Our strategy goes as follows. Since the effect of Lorentz non-invariance in medium has been 
more extensively studied in dense matter, we first examine the form of the relevant correlator 
in dense matter following Refs. [7^1 EZj. The current correlator G^ v constructed from the 
current defined by 

4 q) =qj»q, or = q l5ltl q, (4.14) 



is given by 



- c oln|Q 2 | + E^ n ' n 



^^[-^W 2 + g mx q v q m + <t<t x g vva + g^g^Q 2 } 



r=2 k=l 



n2k 

•- Vt,t 



x <7 M3 • ■ ■ n^ 2k A 2k+1 

y y Q4fc+r-2 Mi-A^fc 

r a^a u i 2 2k 



9 

T=2 fc = l 

•>2 _ „2 



where Q 2 = —q 2 . r = d — s denotes the twist, and s = 2k is the number of spin indices of 
the operator of dimension d. Here A n,n represents the contribution from the Lorentz invariant 
operators such as A A,i = 5 (-G 2 ) . A 2k+T ' T and C 2k+ l' T are the residual Wilson coefficient 
times matrix element of dimension d and twist r; e.g., A®; 2 ,, ,, ,, is given by 

4£*» w = * < 5T (Q^D, 2 D m D, 4 q)) p , (4.16) 

where we have introduced the symbol <ST which makes the operators symmetric and traceless 
with respect to the Lorentz indices. The general tensor structure of the matrix element of 
A 2k t T ' T is given in Ref. |75j . For k — 2, it takes the following form: 



A, 



P / 

PaPpPxPa - -^-(PaPpgXo + PaPxgpa + PaP a gxp + PpP\9a 



4 

+ ppp<,g a \ + p\Pag a p) + -^{g a p9xu + g &xg $0 + g a *gpx) A 4 (4.17) 



For r = 2 with arbitrary fc, we have |76| : 



4 2fc+2,2 _ ^9 49 , 4G 

^2* — Ly 2,2fc^ i 2fc "T L/ 2,2fc /1 2fc 

r< 2k+2,2 _ ^,9 ig 1 /^G 4G 18 \ 

°2fc — °L,2fc /1 2fc + ^L,2k A 2k > l 4 - i5 J 
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where C\ 2k 



1 + 0(a s ), C 



q,G 
L,2k 



0{a s 



and C 2,2k ~ 



0(a s ) (with the superscripts q and G 



standing respectively for quark and gluon) are the Wilson coefficients in the OPE jTH]. The 
quantities A q n and A„ are defined by 



dx x n [q(x, jj) + q(x, jj)} 



2 / dxx n - 1 G{x^) 



(4.19) 



where q(x,[i) and G(x,/i) are quark and gluon distribution functions respectively. We observe 
that ()4.15|) consists of three classes of terms: One is independent of the background, i.e., 
density in this case, the second consists of scalar operators with various condensates (0) p and 
the third is made up of non-scalar operators whose matrix elements in dense matter could not 
be simply expressed in terms of various condensates (0) p . 

It is clear that Eq. ()4.15|) is a general expression that can be applied equally well to heat- 
bath systems. Thus we can simply transcribe (j4.15j) to the temperature case by replacing 
the condensates (0) p by (O)t and the quantities A 2 ^+ T £ k and C^+Tg by the corresponding 
quantities in heat bath. (We show the matching condition on the bare pion velocity at finite 
density in Appendix |FJ) The higher the twist of operators becomes, the more these operators 
are suppressed since the dimensions of such operators become higher and the power of 1/Q 2 
appear. Thus in the following, we restrict ourselves to contributions from the twist 2 (r = 2) 
operators. Then the temperature-dependent correlator can be written as 



i 



2tt 2 



^)ln 

7T 



f< 2 



6Q 4 



—G 2 

71 



2ira 



Q 6 

~9Q> 



u,d, 

K u lp X a u + d^X a d) Y]q^\ a q 



T 



gWiq«qM 



q^q^g^ 2 



\—A 4 ' 2 + —n^n^ A 6 ' 2 
lg 4 /i mM2 + Q8 q q A M 



(4.20) 



where G^ is constructed from the axial-vector current associated with the iso-triplet channel 
defined by 

1 



J< 



5/j 



(4.21) 



and we keep terms only up to the order of 1/Q 8 for A 2 ^ 2 ^. The A a denote the SU(3) color 
matrices normalized as tr[A a A 6 ] = 2S ab . Here we have dropped the terms with C 2k ^ 2 ' 2 in the 
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non-scalar operators since they are of higher order in both l/(Q 2 ) n and a s compared to the 
terms in the first line of Eq. (|4.20j) . The temperature dependence of A A ' 2 and , 



implicit in Eq. (|4.20)l . will be specified below. 

In order to effectuate the Wilsonian matching, we should in principle evaluate the conden- 
sates and temperature-dependent matrix elements of the non-scalar operators in Eq. ([4.20j) at 
the given scale A and temperature T in terms of QCD variables only. This can presumably be 
done on lattice. However no complete information is as yet available from model-independent 
QCD calculations. We are therefore compelled to resort to indirect methods and we adopt 
here an approach borrowed from QCD sum-rule calculations. 

Let us first evaluate the quantities that figure in Eq. (j4.20j) at low temperature. In low 
temperature regime, only the pions are expected to be thermally excited. In the dilute pion- 
gas approximation, (O)t is evaluated as 

(0)t * (O)o + g J ^^(^ a (p}\0\^(p))n B (e/n (4-22) 



where e = a/p 2 + vn\ and n B is the Bose- Einstein distribution. As an example, we consider the 
operator of (r, s) = (2,4) that contributes to both G L ^ and G^\ Noting that G^(q ,q) = 
GAoo/f, we evaluate G A oo(qo,q)- 

XQY^f£n B (e/T), (4.23) 

where A^l^J is given by #4 



^afSXa 



P 2 , 

PaPpPxPa ~ ^{PaPpgXa + PaPx9f3a + P a Po9Xfi + Pf3Px9aa 
o 

p 4 1 

+ P/3PagaX + VWa9a.fi) + ^r(0a/90Aff + 9a\9fa + &*t09a)J A\ , (4.24) 

where A\ carries the temperature dependence. Taking m 2 = 0, we see that the terms with p 2 
and p 4 in Eq. (|4.24|) are zero. 

From Eqs. ijOil) . and (jOl, we obtain 

(OPE)L(O)/ 2\ _ _1 nFp (OPE)(0) 



{-(1 ) - —9 ^A^v 



3 
-1 

T 



_2tt 2 V 7T/ ) 6Q 4 \ vr It 

# 4 For the general tensor structure of the matrix elements with a polarized spin-one target, say, along the 
beam direction in scattering process, see Ref. |78) . 
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2iia s 

w 



T 



u,d, 



7T 



+ 30Q 4 ^ 2 



2 T 4 A7r ( u +d) 167r 4 T 6 „■(„ 



63 g 



4" 



(4.25) 



^-,{ope)l(i) ^.gj^gg f n ow i n g form 



32 A T 6 .-n-^+ji) 



&A ~ 105* Q^ 4 



(4.26) 



We now proceed to estimate the pion velocity by matching to QCD. 

First we consider the matching between (^^ ILS ) i ( ) anc j G^ PE ) L( W. From Eqs. ()4.12j) and 
(|4~25|) . we obtain 



r/ 



2\ G A 



(HLS)L(O) 



7r,bare 7r,bare 



Q 2 



2\ G A 



(OPE)L(O) 



7T 



2tt 2 (^G 2 ) T+ ^U08a s (qq) 2 T 



Q 4 



27 g e 



7T 2 T 4 ^(u+d) , 167T 4 T 6 A K( u +d) 



15 Q 



4^2 



+ 



21 Q 



4" 



(4.27) 



Matching them at Q 2 = A 2 , we obtain 



7r,bare 7r,barc 

A 2 



8tt 2 



1 + ^ 

7T 



2tt 2 (^G 2 ) 



7T 2 T 4 



.4 



7r(M+d) 



15 A 4 2 



3 A 4 
16vr 4 T 6 

^Ta6 



T '-7T 3 



1408 a s (qq) 



2 i 



.4 



7r(u+a 



(4.28) 



Next we consider the matching between q^ ls ^ l W anc [ q(ope)l(i) ^ p rom ^q S _ (gTEH) an d ()4.26|) . 
we have 

T?t Tps /I T/2 \ qo 

(4 29) 



32 /T a7T (u+g 
— A A 



A 2 ~ 105 A 6 4 

Noting that the right-hand-side of this expression is positive, we verify that 

Vtt, bare ^ 1 



(4.30) 



which is consistent with the causality. 

The bare pion velocity can be obtained by dividing Eq. (J4.29)) with Eq. ()4.28j) . What we 
obtain is the deviation from the speed of light: 

1 32 ,T & 



5 bare _ 1 V; )bare - — — tt 4 — A A 



7r(u+Q 



(4.31) 
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This should be valid at low temperature. We note that the Lorentz non-invariance does not 
appear when we consider the operator with s — 2, and that the operator with s = 4 generates 
the Lorentz non-invariance. This is consistent with the fact that G L A including up to the 
operator with s = 2 is expressed as the function of only Q 2 jZH] • Equation (j4.31j) implies that 
the intrinsic temperature dependence starts from the 0(T 6 ) contribution. On the other hand, 
the hadronic thermal correction to the pion velocity starts from the 0(T A ). [There are 0(T 2 ) 
corrections to [/*] 2 and [f%f%], but they are canceled with each other in the pion velocity. See 
subsection 13.4.21 ] Thus the hadronic thermal effect is dominant in low temperature region. 



Chapter 5 

Vector Manifestation in Hot Matter 



The vector manifestation (VM) is a new pattern for Wigner realization of chiral symmetry, in 
sharp contrast to the standard scenario of chiral symmetry restoration. In order to clarify the 
difference between the standard scenario and the VM, we consider the chiral representations 
of the low-lying mesons. 

In the broken phase, the eigenstate of the chiral representation under SU(3)i x SU(3)r 
does not generally agree with the mass eigenstate due to the existence of the Nambu-Goldstone 
bosons: There exists a representation mixing. Then the scalar, pseudoscalar, vector and axial- 
vector mesons belong to the following representations |8T?1 IHT] : 

\s) = |(3,3*) ©(3*, 3)), 

|tt) = |(3,3*)©(3*,3))sin^+|(l,8)©(8,l))cos^, 

\p) = |(1, 8) ©(8,1)), 

= |(3, 3*) ©(3*, 3)) cos^- |(1, 8) ©(8,1)) sin V, (5.1) 

where ip denotes the mixing angle and is determined as ip ~ 45°. 

Now we consider the chiral symmetry restoration, where it is expected that the above 
representation mixing is dissolved. From Eq. (|5.1|) . one can easily see that there are two 
possibilities for pattern of chiral symmetry restoration. One possible pattern is the case where 
cos ip — > when we approach the critical point. In this case, the pion belongs to | (3, 3*)ffi(3*, 3)) 
and becomes the chiral partner of the scalar meson: 

|tt) = |(3,3*) ©(3*, 3)), 

\ s ) = |(3,3*) © (3*, 3)) for cos^ -> 0. (5.2) 

The vector and axial-vector mesons are in the same multiplet |(1,8) © (8,1)). This is the 
standard scenario of chiral symmetry restoration. 
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Another possibility is the case where sinijj — > when we approach the critical point |21j . 
In this case, the pion belongs to pure |(1,8) © (8, 1)) and so its chiral partner is the vector 
meson: 

|tt) = 1(1,8)0(8,1)), 

\p) = |(1,8)©(8,1)) for sin^^O. (5.3) 

The scalar meson joins with the axial-vector meson in the same representation |(3, 3*)©(3*, 3)). 
This is nothing but the VM of chiral symmetry. 

In order to formulate the VM, we need a theory including both pions and vector mesons. 
One of such theories is the one based on the hidden local symmetry (HLS). In the following 
sections, we will show how the VM is formulated at the critical point in the framework of the 
HLS theory. We also study the predictions of the VM in hot matter. 



5.1 Conditions for Bare Parameters 

In this section, we summarize the conditions for the bare parameters obtained in Ref. [T4*] 
through the Wilsonian matching at the critical temperature. 



5.1.1 Case with Lorentz invariance 

We consider the Wilsonian matching near the chiral symmetry restoration point. Here we 
assume that the order of the chiral phase transition is second or weakly first order, and thus 
the quark condensate becomes zero continuously for T — > T c . First, note that the Wilsonian 
matching condition ()4.9|) provides 



*?(A;T C ; 
A 2 



8vr 2 



7l 



+ 



2tt 2 (fG, v G^) Tc 



A 4 



which implies that the matching with QCD dictates 



^0 



(5.4) 



J£(A;T C )^0 



(5.5) 



even at the critical temperature where the on-shell pion decay constant vanishes by adding the 
quantum corrections through the RGE including the quadratic divergence [HI] and hadronic 
thermal corrections, as we will show in section 15.31 As was shown in Ref. [22] for the VM in 
dense matter, the Lorentz non-invariant version of the VM conditions for the bare parameters 
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are obtained by the requirement of the equality between the axial-vector and vector current 
correlators in the HLS, which is valid also in hot matter (see next subsection): 

• rpt \ 2 / rps \ 2 



t / * cr,bare \ T^T C 1 s / -^bare \ T— >T C 1 

a barc=|-^ J ~» 1, «bare=lT^ J ~+ h (5-6) 

V 7r,bare/ \ 7r,bare/ 

^T.barc 0, g L Mre 0, (5.7) 

where a bare , «-bare; <?T,bare an d (?L,bare are the extensions of the parameters abare and gbare in the 
bare Lagrangian with the Lorentz symmetry breaking effect included as in Appendix A of 
Ref. [221 

When we use the conditions for the parameters a t,s in Eq. (|5.6|) and the above result that 
the Lorentz symmetry violation between the bare pion decay constants F^ s hare is small, we 
can easily show that the Lorentz symmetry breaking effect between the temporal and spatial 
bare sigma decay constants is also small, F* baxe ~ -P^bare- While we cannot determine the 
ratio g l, ham / 9t, ham through the Wilsonian matching since the transverse mode of vector meson 
decouples near the critical temperature. * x However this implies that the transverse mode is 
irrelevant for the quantities studied in this thesis. Therefore in the present analysis, we set 
9L,have = 9T,hare fc> r simplicity and use the Lorentz invariant Lagrangian at bare level. In the 
low temperature region, the intrinsic temperature dependences are negligible, so that we also 
use the Lorentz invariant Lagrangian at bare level as in the analysis by the ordinary chiral 
Lagrangian in Ref. [TT] . 

At the critical temperature, the axial-vector and vector current correlators derived in the 
OPE agree with each other for any value of Q 2 . Thus we require that these current correlators 
in the HLS are equal at the critical temperature for any value of Q 2 around A 2 . As we discussed 
above, we start from the Lorentz invariant bare Lagrangian even in hot matter, and then the 
axial-vector current correlator G^ LS ^ and the vector current correlator G^ LS ^ are expressed 
by the same forms as those at zero temperature with the bare parameters having the intrinsic 
temperature dependences: 

Gf LS \Q 2 ;T) = S^-2z 2 (A;T), 



^ tl In Ref. |22j . the analysis including the Lorentz non-invariance at bare HLS theory was carried out. Due 
to the equality between axial-vector and vector current correlators, (gi.barc, 5T,barc) - > (0,0) is satisfied when 
we approach the critical point. This implies that at the bare level the longitudinal mode becomes the real NG 
boson and couples to the vector current correlator, while the transverse mode decouples. Furthermore <?l — > 
is a fixed point for the RGE ^2]- Thus in any energy scale the transverse mode decouples from the vector 
current correlator. For details, see next subsection. 
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G r»(Q-T) = ^(A;r)[i-2^(A;r), 3 (A;r)] 

v W ' ' M p 2 (A;T)+Q 2 V ' ' V ' 

By taking account of the fact F 2 (A-T C ) 7^ derived from the Wilsonian matching condition 
given in Eq. (|5.4jl . the requirement G^ ILS ' ) = G^ ILS ' ) is satisfied only if the following conditions 
are met [T4*] : 

<7(A;T) T ^ Tc 0, (5.9) 
a(A;T) = F 2 (A;T)/F 2 (A;T) r =?' 1, (5.10) 
Z!(A;T)-z 2 (A;T) ^ 0. (5.11) 

Note that the intrinsic thermal effects act on the parameters in such a way that they become 
the values of Eqs. (f5~3]) - (j5~TT]l . 

Through the Wilsonian matching at non-zero temperature mentioned above, the parame- 
ters appearing in the hadronic thermal corrections calculated in section 13.21 have the intrinsic 
temperature dependences: F n , a and g appearing there should be regarded as 

F w = F n {fi = 0;T) , 
a = a (/i = M p (T);T) , 

g = g(fi = M p (T);T) , (5.12) 

where M p is determined from the on-shell condition: 

M 2 p ee M 2 p (T) = aQi = M p , T)g 2 {» = M p - T)F 2 ^ = M p - T) . (5.13) 

From the RGEs for g and a in Eas. (iri29l and (ITT28j) . we find that (g, a) = (0, 1) is the fixed 
point. Therefore, Eqs. ()5.9|) and (|5.10|l imply that g and a at the on-shell of the vector meson 
take the same values: 

a(fi = M p (T);T) T -^ c 1 , 

g^ = M p (T)-T) T ^ c 0, (5.14) 

where the parametric vector meson mass M P (T) is determined from the condition ()5.13j) . The 
above conditions with Eq. ()5.13|) imply that M P (T) also vanishes: 

M P (T) T ^S C . (5.15) 
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5.1.2 Case without Lorentz invariance 

In this subsection, we start from the Lagrangian with Lorentz non-invariance, and requiring 
that the axial-vector current correlator be equal to the vector current correlator at the critical 
point, we obtain the conditions for the bare parameters. Then we show that the conditions 
are satisfied in any energy scale, which is protected by the fixed point of the RGEs. 

First we show the Lorentz non-invariant version of the conditions satisfied at the critical 
temperature for the bare parameters, following Ref. [22] where the conditions for the current 
correlators with the bare parameters in dense matter were presented. We consider the matching 
near the critical temperature. At the chiral phase transition point, the axial- vector and vector 
current correlators must agree with each other: G^hls) = Gy^ ULS ^ and G^( H ls) = ^v(hls)- 
These equalities are satisfied for any values of po and p around the matching scale only if the 
following conditions are met: 

fl baro — ¥ 1) a barc ~~ > 1) 

#L,bare ~» 0, #T,barc ~ ► for T —> T c . (5.16) 

This implies that at bare level the longitudinal mode of the vector meson becomes the real 
NG boson and couples to the vector current correlator, while the transverse mode decouples. 

Next, we show that the conditions for the bare parameters for T ~^T C are satisfied in any 
energy scale and that this is protected by the fixed point of the RGEs. 

It was shown that the HLS gauge coupling g = is a fixed point of the RGE for g at one- 
loop level jSOJEl]. The existence of the fixed point g = is guaranteed by the gauge invariance. 
This is easily understood from the fact that the gauge field is normalized as = gp^. In the 
present case without Lorentz symmetry, the gauge field is normalized by g^ as = fl^/fy and 
thus gi = becomes a fixed point of the RGE for g^. 

Provided that gi = is a fixed point, we can show that a 1 = a s = 1 is also a fixed point of 
the coupled RGEs for a 1 and a s as follows: We start from the bare theory defined at a scale A 
with a^ are = a bare = 1 (and gz — 0). The parameters a* and a s at (A — 5 A) are calculated by 
integrating out the modes in [A — 8 A, A]. They are obtained from the two-point functions of 
Afj, and V M , denoted by LT^ and ITj^. We decompose these functions into 

= u'vnij + W - u'vonij + iT n ln + p r n Ln> ( 5 - 17 ) 

where u^u v ', {g^ v —u fJ- u u ) , P£ v and P^ v denote the temporal, spatial, longitudinal and transverse 
projection operators, respectively. The parameters a* and a s are defined by a 1 = Iln/IT^ , a s = 
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Figure 5.1: Diagrams for contributions to H 1 ^ at one-loop level. 
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Figure 5.2: Diagrams for contributions to Ely at one-loop level. 

njj/n^ [TU] . These expressions are further reduced to 

1 



a* (A - 5 A) 
a s {A - 5A) 



bare 



''bare 



+ 



+ 



(Ft 

\ 7r,bare 
1 



) 2 L 



n UA-A-6A)-al 3Xe U\(A ] A-SA) 



F l F s 

7r,bare 7r,bare 



IIjj(A;A-5A) 



''bare 



ni(A;A-5A) , (5.18 



where I1^ S ||(A; A — 5 A) denotes the quantum correction obtained by integrating the modes 
out between [A — 5 A, A]. We show the diagrams for contributions to 11^ and Ilj^ at one- 
loop level in Figs. 15.11 and 15.21 The contributions (a) in Fig. 15.11 and (a) in Fig. 15.21 are 
proportional to g\ bare . The contributions (c) in Fig. 15.11 and (d) in Fig. 15.21 are proportional 
to (a barc — 1). Taking gz^bare = and a bare = a bare = 1, these contributions vanish. We note 
that a (i.e., longitudinal vector meson) is massless and the chiral partner of pion at the critical 
temperature. Then the contributions (b) and (c) in Fig. 15.21 have a symmetry factor 1/2 
respectively and are obviously equal to the contribution (b) in Fig. 15.11 i.e., 
Thus from Eq. ()5.18j) . we obtain 



n 



(b)+(c)tiu 



a*(A - SA) — a bare = 1, 
a s (A-5A) =< arc = 1. 



(5.19) 



This implies that a 1 and a s are not renormalized at the scale (A — 5 A). Similarly, we include 
the corrections below the scale (A — 5A) in turn, and find that a 1 and a s do not receive the 
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quantum corrections. Eventually we conclude that a* = a s = 1 is a fixed point of the RGEs 
for a* and a s . 

From the above, we find that (gi, a 1 , a s ) = (0, 1, 1) is a fixed point of the combined RGEs 
for gL, a 1 and a s . Thus the VM condition is given by 

9l -> 0, 

a* -> 1, a s -> 1 for T -> T c . (5.20) 

The vector meson mass is never generated at the critical temperature since the quantum 
correction to M 2 p is proportional to g\. Because of ql — ► 0, the transverse vector meson at 
the critical point, in any energy scale, decouples from the vector current correlator. The VM 
condition for a* and a s leads to the equality between the ir and a (i.e., longitudinal vector 
meson) velocities: 

(K/K) 4 = {F°Fl/F°Fl) 2 

= a'/a'^l. (5.21) 

This is easily understood from a point of view of the VM since the longitudinal vector meson 
becomes the chiral partner of pion. We note that this condition V a — V n holds independently of 
the value of the bare pion velocity which is to be determined through the Wilsonian matching. 

5.2 Vector Meson Mass 

In this section, we briefly review that the vector manifestation (VM) in hot matter can be 
formulated following Ref . ^H] • Including the intrinsic temperature dependences of the param- 
eters near the critical temperature determined in the previous section, we show that the vector 
meson mass vanishes at the critical temperature. 

Let us study the vector meson pole mass near the critical temperature. As shown in 
section l5~T| the parametric vector meson mass M p vanishes at the critical temperature, which 
is driven by the intrinsic effects. Then, near the critical temperature we should take M p <^ T 
in Eq. ()3.120j) instead of T ^ M p which was taken to reach the expression in Eq. (j3.121j) in 
the low temperature region. Thus, by noting that 

G 2 {M p -T) J 2 (T), 
J!(M p ,T) M -C° J 2 (T), 
M p 2 Fl{M p -M p -T) A/ C° 0, (5.22) 
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the pole mass of the vector meson at T < T c becomes 

m p 2 (T) = M/ + A^ 2 ±^/ 2 (T) 

= M^ + N fg ^-^f^. (5.23) 

In the vicinity of a ~ 1 the hadronic thermal effect gives a positive correction, and then the 
vector meson pole mass is actually larger than the parametric mass M p . However, the intrinsic 
temperature dependences of the parameters obtained in section I5~T1 lead to g — »• and M p — > 
for T — > T c . Then, from Eq. ()5.23|) we conclude that the pole mass of the vector meson m p 
also vanishes at the critical temperature: 

m p {T) ^ for T -> T c . (5.24) 

This implies that the VM is formulated at the critical temperature in the framework of the 
HLS theory, which is consistent with the picture shown in Refs. [TUl El EH H] • 



5.3 Pion Decay Constants 

In this section, we show how the two decay constants, /* and vanish at the critical tem- 
perature following Ref. [To] . 

As we discussed around Eq. (15.14)) . at the critical temperature the intrinsic thermal effects 
lead to (g,a) — > (0,1) which is a fixed point of the coupled RGEs. Then the RGE for F n 
becomes 



J**- ~ -^u 2 (5 25) 

This RGE is easily solved and the relation between F 7r (A; T c ) and F^O; T c ) is given by 

^ 2 (0; T c ) = F 2 (A; T c ) - ^pA 2 . (5.26) 

Finally we obtain the pion decay constants as follows: 

/*F = F 2 (0;T) + nl(p,p;T), 

/^ = F 2 (0;T) + nl(p,p;T), (5.27) 

where the second terms are the hadronic thermal effects. In the VM limit (g, a) — > (0, 1) the 
temperature dependent parts become 

nl(p" p; T) T ^ - ^-T c 2 , ITl(p- p; T) T ^ - ^T c 2 . (5.28) 
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From Eqs. (j3.134j) . ()5.27|) and ()5.28|) . the order parameter f n becomes 



flip; T) T ^ fl(p; T c ) = F*(0; T c ) - ^-T c 2 . (5.29) 



Since the order parameter /„. vanishes at the critical temperature, this implies 



Thus we obtain 



F 2 (0;T) ^ F^;T C ) = ^-T 2 . (5.30) 



(fJ^O, fiK^O. (5.31) 



From Eq. (j3.146|) or equivalently Eq. ()3.141|1 . the above results imply that the temporal and 
spatial pion decay constants vanish simultaneously at the critical temperature |15j : 

ft(T c ) = f°(T c ) = 0. (5.32) 

Comparing Eq. (j5.29|) with the expression in the low temperature region in Eq. (j3.147J) where 
the vector meson is decoupled, we find that the coefficient of T c 2 is different by the factor 
|. This is the contribution from the cr-loop (longitudinal p-loop). In the low temperature 
region the 7r-loop effects give the dominant contributions to f n (T) and the p-loop effects are 
negligible. However by the vector manifestation the p contributions become essentially equal to 
the one of it, are then incorporated into fir(T) near the critical temperature. 



5.4 Predictions of the Vector Manifestation 

In this section, we summarize the predictions of the VM in hot matter studied in Refs. [T^ 

[IS EH EH EE]. 



5.4.1 Critical temperature 

In this subsection we estimate the value of the critical temperature T c where the order param- 
eter / 2 vanishes. 

We first determine T c by naively extending the expression ()3. 147)1 to the higher temperature 
region to get j , c ( hadron ) — igo MeV for Nf = 3. However this naive extension is inconsistent with 
the chiral restoration in QCD since the axial vector and vector current correlators do not agree 
with each other at that temperature. As is stressed in Ref. (Hj, the disagreement between two 
correlators is cured by including the intrinsic thermal effect. As can be seen from Eq. (J4.9)) . 
the intrinsic temperature dependence of the parameter F n is determined from {^G ^G^) T 
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Aqcd 


0.30 


0.35 


A 


0.8 


0.9 


1.0 


1.1 


0.8 


0.9 


1.0 


1.1 


T c 


0.20 


0.20 


0.22 


0.23 


0.20 


0.21 


0.22 


0.24 


Aocn 


0.40 


0.45 


A 


0.8 


0.9 


1.0 


1.1 


0.8 


0.9 


1.0 


1.1 


T 


0.21 


0.22 


0.23 


0.25 


0.22 


0.23 


0.24 


0.25 



Table 5.1: Values of the critical temperature for several choices of Aqcd and A. The units of 
Aqcd, A and T c are GeV. 

and (qq) T , and gives only a small contribution compared with the main term 1 + — . However 
it is important that the parameters in the hadronic corrections have the intrinsic temperature 
dependences as (a,g) —* (1,0) for T — > T c , which carry the information of QCD. Actually the 
inclusion of the intrinsic thermal effects provides the formula ()5.29j) for the pion decay constant 
at the critical temperature, in which the second term has an extra factor of 1/2 compared with 
the one in Eq. (13.1 4711 . 

In Ref. jTJj we determined the critical temperature from f n (T c ) = and estimated the 
value which is dependent on the matching scale A: From Eq. ()5.30)1 we obtain 

/ 94 

T C =J W F^-T C ). (5.33) 
V iV / 



Using Eqs.(JE2J and (jQBj) we get [Hj 



T r 



1+° 



2tt 2 (fG^G^Tc N 



A V N f 7T 2 L Tr 3 A 4 4 J 



(5.34) 



We would like to stress that the critical temperature is expressed in terms of the parameters 
appearing in the OPE. We evaluate the critical temperature for Nf = 3. The gluonic conden- 
sate at T c is about half of the value at T = |HSl H and we use (^G^G^) Tc = 0.006 GeV 4 . 
We show the predicted values of T c for several choices of Aqcd and A in Table 15.11 Note 
that the Wilsonian matching describes the experimental results very well for Aqcd — 0.4 GeV 
and A = 1.1 GeV at T = ^3]. At non-zero temperature, however, the matching scale A 
may be dependent on temperature. The smallest A in Table 15.11 is determined by requiring 
(2tt 2 /3) (f G^G^}/A 4 < 0.1. 

We should note that the above values may be changed when we adopt a different way to 
estimate the matrix elements of operators in the OPE side, e.g., an estimation with the dilute 
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pion gas approximation [75] or that by the lattice QCD calculation jSl]. Even when we choose 
one way to estimate the matrix elements in the OPE, some temperature effects are supposed 
to be left out due to the truncation of the OPE to neglect higher order operators, inclusion of 
which will cause a small change of the above values of the critical temperature. The important 
point is that as a result of the Wilsonian matching, T c is obtained as in Eq. (J5.34J1 in terms of 
the quark and gluonic condensates, not hadronic degrees of freedom. It is expected that the 
value of T c may become smaller than that obtained in this paper by including the higher order 
corrections. 



5.4.2 Axial- vector and vector charge susceptibilities 

In this subsection, we address the issue of what the relevant degrees of freedom can be at the 
chiral transition induced by high temperature and their possible implications on observables 
in heavy-ion physics. In doing this, we focus on the vector and axial-vector susceptibilities 
very near the critical temperature T c . The issue of what happens at high density is discussed 
in |BSj. 

Susceptibility is related with the fluctuation of a conserved quantity. The statistical expec- 
tation value of a conserved operator O is given by 

(O) = -^-^ 1, (5.35) 

where 7i denotes the Hamiltonian of system and \x is the chemical potential associated with 
O and we define the partition function Z as 

Z = Tr[e-^°)/ T ]. (5.36) 

The mean square deviation of O is 

(50) 2 = {0 2 }-(0} 2 
djj, 



T«£>. (5.37) 



Then we define the susceptibility of O as follows: 

8(0), 



X(T) 



dfx 

-l/T 



1/14=0 



J dr J d 3 x(O{r,x)O{0,6)). (5.38) 



Consider the vector isospin susceptibility (VSUS) xv and the axial-vector isospin suscep- 
tibility (ASUS) xa defined in terms of the vector charge density Vq{x) and the axial- vector 
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charge density A^(x) by the Euclidean correlators: 

-l/T 



S ab xv 
S ab XA 



dr / d 6 x(V o a (T,x)V o b (O,O)) , 



l/T 



dr / d 6 x (A%(t, x)A b (0, 0)) p 



where ( )p denotes thermal average and 



T 



(5.39) 
(5.40) 

(5.41) 



with the quark field ip and the r a Pauli matrix the generator of the flavor SU(2). 

The axial- vector susceptibility Xa(T) and the vector susceptibility Xv(T) for non-singlet 
currents * 2 are given by the 00-component of the axial-vector and vector current correlators 
in the static-low-momentum limit: 



XA (T) =2N f \im \im [G™(p ,p;T)] , 
Xv (T) = 2N f \im Llim [G$(p ,fiT)] , 



(5.42) 



where we have included the normalization factor of 2Nf. Using the current correlators given 
in Eqs. (j3.1()()|) and (j3.11()j) and noting that \im Po ^ Q Pl° = lim po ^ p 2 /p 2 = — 1, we can express 
Xa(T) and Xv(T) as 



X a(T) = -2N f lim lim [lli(p ,p;T) - Tl\(p^p,T 

p^Opo^O 



Xv(T) 



2Nf lim lim 



p^Opo^O 



% (n£ + 2n£, 



+ n 



(5.43) 



where for simplicity of notation, we have suppressed the argument (po, p; T) in the right-hand- 
side of the expression for xv(T). In HLS theory at one-loop level, the susceptibilities read 



XA (T) = 2N f 
Xv (T) = -2N f 1: 



F 2 (0) + lim lim {nl(p ,p; T) - ni(p ,p; T)} 

p^Opo^O 



p^Opo^O 



(a(Q)F^o) + W v ) (n£ + 2n 



+ n 



(5.44) 



a(0)i^(0) + fl* v -n^ 
where the parameter a(0) is defined by (see section 15. 4. 3|) 

rn , = 4 vac)t ( Po = o,p = o) 4 vac) -( Po = 0,p = 0) 

01 J ^ 2 (0) F*(0) 

We show the tree and one-loop contributions to xa and xv in Figs. 15.31 and 15.41 It follows 



(5.45) 
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Figure 5.3: Diagrams for contributions to xa- 




Figure 5.4: Diagrams for contributions to xv- 



from the static-low- momentum limit of (IIj_ — 11^) given in Eq. 
susceptibility Xa(T) takes the form 



HI that the axial-vector 



Xa{T) 



2N f 



F 2 (0) - N f Jl(0; T) + Nfa .]\{M p - T) 



iY/ ^|{ j!l(M ^ T) " j!l(0;T) 



(5.46) 



Near the critical temperature (T — > T c ), we have M p — »• 0, a — > 1 due to the intrinsic tempera- 
ture dependence in the VM in hot matter ^3]. Furthermore, from Eq. (J5.3()j) . we see that the 
parameter -F^(O) approaches ^T 2 for T — > T c . Substituting these conditions into Eq. (J5.46I) 
and noting that 



lim 

M p ^0 



we obtain 



-^-\J\(M P -T)-J\{Q-T) 

p 



N 2 f 

Xa(T c ) = -fT 2 



-Jf(0;T) = —T 2 
2 iv > i 24 



(5.47) 



(5.48) 



To obtain the vector susceptibility near the critical temperature, we first consider a(0)F 2 (0) + 

fly appearing in the numerator of the first term in the right-hand-side of Eq. ()5.44j) . Using 

# 2 We will confine ourselves to non-singlet (that is, isovector) susceptibilities, so we won't specify the isospin 
structure from here on. 



CHAPTER 5. VECTOR MANIFESTATION IN HOT MATTER 



69 



Eq. (jG.9|) . we get for the static-low-momentum limit of a(0)F 2 (0) + Yly as 
lim lim [a(0)F 2 (0) + U^p; T)] 

p^Opo^O 

Nf r ~ ~ ~ i 

= a(0)F 2 (0) - -f \2J\{M P - T) - J*(M p] T) + a 2 J 2 (0; T)J . (5.49) 

From Eq. (I3~T5%1) we can see that a(0) -> 1 as T -> T c since F 2 (M p ) -> F 2 (0) and M p -> 0. 
Furthermore, F 2 (0) -> ^-T 2 as we have shown in Eq. (|5.30j) . Then, the first term of Eq. ()5.49|) 
approaches -^T 2 . The second term, on the other hand, approaches ~-^T 2 as M p — > and 
a — > 1 for T — > T c . Thus, we have 

lim lim U(0)F 2 (0) + fl^(po,p; T)l — >• . (5.50) 

p^Opo^O L v J T-»T C 

This implies that only the second term in the right-hand-side of Eq. ()5.44j) contributes to 
the vector susceptibility near the critical temperature. Thus, taking M p —> and a — > 1, in 
Eq. (RTT0|) . we obtain 

Xv(T c ) = -fTl , (5.51) 

which agrees with the axial- vector susceptibility in Eq. (|5.48|) . This is a prediction, not an 
input condition, of the theory. For Nf = 2, we have 

Xa(T c )= X v(T c ) = ^T 2 c . (5.52) 

The result xv(T c ) = f^ 2 is consistent with the lattice result as interpreted in It is 
interesting to note that the RPA result obtained in jH] in NJL model in terms of a quasi- 
quark-quasi-antiquark bubble is reproduced quantitatively by the one-loop graphs in HLS 
with the VM. 

One might think that the VSUS diverges at T c because of the massless p meson pole. 
However the resultant VSUS is finite, and it will imply that the screening mass of p meson is 
also finite. 

It should be noticed that the pion pole effect does not contribute to the ASUS in Eq. (|5.48|) 
since the pion decay constant /* vanishes at the critical temperature as we have shown in 
section 15. 3[ and that the contribution to the ASUS comes from the non-pole contribution ex- 
pressed in Fig. 15.31 In three diagrams, the third diagram in Fig. l5.3r c) is proportional to (1 — a) 
and then it vanishes at the critical temperature due to the VM. Similarly, since the transverse 
p decouples at the critical point in the VM [^2 E] > the first diagram in Fig. I5.3f a) does not 
contribute. Thus, the above result for the ASUS in Eq. (J5.48)) comes from only the contribu- 
tion generated via the (longitudinal) vector meson plus pion loop [see Fig. 15.3( b)] . Similarly, 
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the vector meson pole effect to the VSUS vanishes at the critical temperature as shown in 
Eq. (|5.50[) . and the contribution to the VSUS is generated via the pion loop [Fig. 15.4( c)] and 
the longitudinal vector meson loop [Fig. 15.4( b)]. * 3 Since the longitudinal vector meson be- 
comes massless, degenerate with the pion as the chiral partner in the VM, loop contribution to 
the ASUS becomes identical to that to the VSUS. Thus, the massless vector meson predicted 
by the VM fixed point plays an essential role to obtain the above equality between the ASUS 
and the VSUS. 

In the present analysis, our aim is to show the qualitative structure of the ASUS and 
the VSUS in the VM, i.e., the equality between them is predicted by the VM. In order to 
compare our qualitative results with the lattice result, we need to go beyond the one-loop 
approximation. We note here that there is a result from a hard thermal loop calculation 
which gives xv(T c ) ~ 1-3T C 2 [HE]- However this result cannot be compared to ours for two 
reasons. First we need to sum higher loops in our formalism which may be done in random 
phase approximation as in [ST]. Second, the perturbative QCD with a hard thermal loop 
approximation may not be valid in the temperature regime we are considering. Even at 
T T c , the situation is not clear as pointed out in Ref. [HEj- 

5.4.3 Violation of vector dominance 

In this subsection, we shed some light on the validity of the vector dominance (VD) of elec- 
tromagnetic form factor of the pion near the critical temperature. In several analyses such 
as the one on the dilepton spectra in hot matter carried out in Ref. [Hj, the VD is assumed 
to be held even in the high temperature region. There are several analyses resulting in the 
dropping mass consistent with the VD, as shown in Refs. [EU EDI- On the other hand, the 
analysis done in Ref. [51] shows that the thermal vector meson mass goes up if the VD holds. 
Thus, it is interesting to study what the VM predicts on the VD. In the present analysis we 
present a new prediction of the VM in hot matter on the direct photon-7r-7r coupling which 
measures the validity of the VD of the electromagnetic form factor of the pion. We find that 
the VM predicts a large violation of the VD at the critical temperature. This indicates that the 
assumption of the VD may need to be weakened, at least in some amounts, for consistently 
including the effect of the dropping mass of the vector meson. 

In Ref. [72] it has been shown that VD is accidentally satisfied in Nf = 3 QCD at zero 
temperature and zero density, and that it is largely violated in large Nf QCD when the 

# 3 Note that the contribution from Fig. 15. 4f a) vanishes since the transverse p decouples and that the one from 
Fig.JOJd) also vanishes since it is proportional to (1 — a). 
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VM occurs. At non-zero temperature there exists the hadronic thermal correction to the 
parameters. Thus it is nontrivial whether or not the VD is realized in hot matter, especially 
near the critical temperature. Here we will show that the intrinsic temperature dependences of 
the parameters of the HLS Lagrangian play essential roles, and then the VD is largely violated 
near the critical temperature. 

Let us study the validity of the VD near the critical temperature. As we have shown in 
section 15. 3[ near the critical temperature II j_ and U s ± in Eqs. (I3.164J) and ()3.165jl approach the 
following expressions: 

TL ± (q,q;T) -> - ^j-T , 

nl(g,g;T) T =^ - ^-T 2 . (5.53) 

On the other hand, the functions ITjj and IT| in Eq. ()3.166|) in the limit of M p /T — > and 
a — > 1 become 

nJ | (0,0;T)=n^(0,0;T) - —±h{T) = -^T 2 . (5.54) 

Furthermore, from Eq. ()3.158|) . the parameter a(0) approaches 1 in the limit of M p —>■ and 

F 2 AM p )/F^) - 1: 

a(0) -> 1 . (5.55) 

From the above limits in Eqs. (|5.53j) . ()5.54j) and (|5.55|l . the numerators of a*(g; T) and a s (q; T) 
in Eqs. ()3.164|) and ()3.165|) behave as 

n;(0,0;T)-a(0)nl(g,g;T)^0, 

TTJ(0, 0; T) - a{0)TT ± {q, q; T) 0. (5.56) 

Thus we obtain 

a*(?;T),a s (£;T) T ^l , (5.57) 
namely, the vector dominance is largely violated near the critical temperature. 

5.4.4 Pion velocity 

In this subsection, we focus on the pion velocity at the critical temperature and study the 
quantum and hadronic thermal effects based on the VM. The pion velocity is one of the 
important quantities since it is a dynamical object, which controls the pion propagation in 
medium through the dispersion relation. 
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Non-renormalization theorem on pion velocity at T c 

As we mentioned in chapter the intrinsic temperature dependence generates the effect of 
Lorentz symmetry breaking at bare level. Then how does the Lorentz non-invariance at bare 
level influence physical quantities? In order to make it clear, in this section we study the pion 
decay constants and the pion velocity near the critical temperature. 

Following subsection 13.4. 2\ we define the on-shell of the pion from the pole of the longi- 
tudinal component G\ of the axial- vector current correlator. This pole structure is expressed 
by temporal and spatial components of the two-point function W^. The temporal and spatial 
pion decay constants are expressed as follows [T51 EH] : 

(^(p;T)) 2 = n* ± (V;p,p;T), 
fl{p-T)f°{p-T) = Ul(V n p,p;T), (5.58) 

where the on-shell condition po —> V n p was taken. We divide the two-point function IT^ into 
two parts, zero temperature (vacuum) and non-zero temperature parts, as = n^ ac ^+IP_^. 
The quantum correction is included in the vacuum part and the hadronic thermal 

correction is in IT^. In the present perturbative analysis, we obtain the pion velocity as ^Hj 

= v , +n±(Km + m^Mih^mhm, (5 . 59) 

where H±(V K p,p) denotes a possible finite renormalization effect. Note that the renormalization 
condition on V n is determined as Ii.±(V n p,p)\p=o = 0. 

In the following, we study the quantum and hadronic corrections to the pion velocity for 
T — > T c , on the assumption of the VM conditions (|5.2()|) . As we defined above, the two-point 
function associated with the pion velocity v n (p; T) is IT^po,^; T). The diagrams contributing 
to are shown in Fig. 15.11 As mentioned in subsection 15.1.21 diagram (a) is proportional 
to gL and diagram (c) has the factor (a* — 1). Then these contributions vanish at the critical 
point. We consider the contribution from diagram (b) only. 

Quantum correction at T c 

First we evaluate the quantum correction to the vacuum part n^ ac ^ b ^ At!/ . This is expressed 

as 

d n k T^(k;p)T u (-k;-p) 



nf c)(b) ^(po,p)=iV/| 



(5.60) 
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where r M denotes the Ana vertex as 

T^{k;p) = -Va 1 g^u 9 + Vfig 1 * - (2Jfc - p) v . (5.61) 



n • ~ dpi 

'2 



We note that the spatial component of this vertex P has an extra-factor V 2 as compared with 
the temporal one. In the present analysis it is important to include the quadratic divergences 
to obtain the RGEs in the Wilsonian sense as we have already seen in chapter 0] and In this 
section, when we evaluate four dimensional integral, we first integrate over fc from — oo to oo. 
Then we carry out the integral over three-dimensional momentum k with three-dimensional 
cutoff A3. In order to be consistent with ordinary regularization in four dimension [BK j I64 | IT3*] . 
we use the following replacement associated with quadratic divergence: 

Aa - T2 A * = 71 A ' (5 ' 62) 

d n ~ l k 1 A 2 f d n ~ l k k l W „., A 2 

-o J —r. (5.63) 



(2-n) n - 1 k 8vr 2 ' J (2w) n - 1 k 3 8tt 2 ' 
When we make these replacement, the present method of integral preserves chiral symmetry. 

As shown in Appendix A, n^ ac ^ and n^ ac ^ s are independent of the external momentum. 
Accordingly, the finite renormalization effect n_i_ is also independent of the external momentum 
and then vanishes: 

Il±(V n p,p) = 0. (5.64) 

Thus in the following, we take the external momentum as zero. In that case, the temporal and 
spatial components of jjC wac )A«' are expressed as n^ ac ' )t = n^ ac - )00 and n^ ac ' )s = — (5^/3) II 
Taking the VM limit (M p — > and V a — > V^), these components become 

hmn^>o = P = o) - Nf f dkodn ~ lk 4k ° 



vm ± A J i(2vr) n [-kl + V 2 k 2 } 2 



4 J (27r) n ~ l V n k 
N f 1 A 2 



4 K87T 2 ' 

N f^^2 f dk Q d n - x k AkW 



i(2%) n [-k 2 + V 2 k 2 } 2 



N f 4 f d n ~ v k tk j 
j y 1 



4 n J (27T)' 1 - 1 (V n k) s 



2 



4 '"" 8tt 2 ' 
Thus we obtain the temporal and spatial parts as 

limni vac)t (p =p = 0) - - ^ JL^l 



N f •• A 

f V n 5 l >—. (5.65) 



VM 4 
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Figure 5.5: Quadratic running of {F^) 2 and F*F* at T c . 

limn^> 0=p - = 0) = -^V;^. (5.66) 

These quadratic divergences are renormalized by (F* bare ) 2 and F* bare i^ bare , respectively. Then 
RGEs for the parameters (i 7 ^) 2 and F*F£ are expressed as 

Both (i 7 ^) 2 and i^-F* scale following the quadratic running /j?. However, the coefficient of ^ 2 
in the RGE for (F*) 2 is different from that for F*F£ (see Fig. 15. 5j) . 

When we use these RGEs, the scale dependence of the parametric pion velocity is 

dv 2 AKKKK) 2 ) 



— = v- , 

dfi an 
1 Nf 



VJF 1 ) 2 - F l F s — 

V T^\ ± IT J TT 7T T 7" 

V TT 



(F*) 4 (4tt)* 

= 0. (5.69) 

This implies that the parametric pion velocity at the critical temperature does not scale. In 
other words, the Lorentz non-invariance at bare level is not enhanced through the RGEs as 
long as we consider the pion velocity. As we noted below Eq. (|5.61j) . the factor V 2 is in the 
spatial component of the vertex T^. If V a were not equal to V^, the coefficients of running in 
the right-hand-side of Eqs. ()5.67|) and ()5.68|) would change. However, since the VM conditions 
do guarantee V a = V n , the quadratic running caused from A 2 in (i 7 ^) 2 and F^F* are exactly 
canceled in the second line of Eq. ()5.69j) . 
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Hadronic thermal correction at T c 

Next we study the hadronic thermal correction to the pion velocity at the critical temper- 
ature. The temporal and spatial parts of the hadronic thermal correction fi^ contribute to 
the pion velocity, which have the same structure as those of the quantum correction ny ac )'" / j 
except for a Bose-Einstein distribution function. Thus by the replacement of A 2 /(47r) 2 with 
T 2 /12 in TT^ ac ^' s , hadronic corrections to the temporal and spatial parts of ft^ are obtained 
as follows: 

H S n<( P0 ,,;T) = -§±e 

hmni(p ,p;T) = -^-KT c 2 . (5.70) 

Substituting Eq. ()5.70|) into Eq. ()5.59|) with = as shown in Eq. ()5.64|) . we obtain the 
physical pion velocity in the VM as 

2 t^ Tc 2 ni(y^,p;T c )-v; 2 ni(Kp,p;T c ) 
^ J ~* K TO 2 

= V n 2 . (5.71) 

Since the parametric pion velocity in the VM does not scale with energy [see Eq. (|5.fi9j) ]. 
in the above expression is equivalent to the bare pion velocity: 

v n (p; T) = K,barc(T) for T T c . (5.72) 

This implies that the pion velocity in the limit T — > T c receives neither hadronic nor quantum 
corrections due to the protection by the VM. 

In order to clarify the reason of this non-renormalization property, let us recall the fact that 
only the diagram (b) in Fig. 15.11 contributes to the pion velocity at the critical temperature. 
Away from the critical temperature, the contribution of the massive a (i.e., the longitudinal 
mode of massive vector meson) is suppressed owing to the Boltzmann factor exp[— M p /T], and 
then only the pion loop contributes to the pion velocity. Then there exists the 0(T 4 ) correction 
to the pion velocity fB] (see subsection IH.4.2J1 . Near the critical temperature, on the other 
hand, a becomes massless due to the VM since a (i.e., the longitudinal vector meson) becomes 
the chiral partner of the pion. Thus the absence of the hadronic corrections in the pion velocity 
at the critical temperature is due to the exact cancellation between the contribution of pion 
and that of its chiral partner a . Similarly the quantum correction generated from the pion 
loop is exactly canceled by that from the a loop. 
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Non-renormalization property 

Now we consider the meaning of our result ([5.72)1 . Based on the point of view that the 
bare HLS theory is defined from QCD, we presented the VM conditions realizing the chiral 
symmetry in QCD consistently, i.e, (gi, a*, a s ) — > (0, 1, 1) for T — > T c . This is the fixed point 
of the RGEs for the parameters g^, a* and a s . Although both pion decay constants (i 7 ^) 2 and 
F*F£ scale following the quadratic running, (-F*) 2 and F*F* show a different running since 
the coefficient of /i 2 in Eq. 1)5.67)1 is different from that in Eq. (J5.68)) . Nevertheless in the pion 
velocity at the critical temperature, the quadratic running in (F*) 2 is exactly canceled by that 
in F*F£ [see second line of Eq. ()5.69j) ]. There it was crucial for intricate cancellation of the 
quadratic running that the velocity of a (i.e., longitudinal vector meson) is equal to its chiral 
partner, i.e., V a — > V n for T — > T c . Note that this is not an extra condition but a consequence 
from the VM conditions for a* and a s ; we started simply from the VM conditions alone and 
found that V w does not receive quantum corrections at the restoration point. As we showed 
in Eq. ()5.70)1 . the hadronic correction to (F*) 2 is different from that to F*F£. In the pion 
velocity, however, the hadronic correction from (i 7 ^) 2 is exactly cancelled by that from F^F^ 
[see second line of Eq. (15.71)) ] . The VM conditions guarantee these exact cancellations of the 
quantum and hadronic corrections. This implies that (gi, a*, a s , V n ) = (0,1,1, any) forms a 
fixed line for four RGEs of gL,a l ,a s and V n . When one point on this fixed line is selected 
through the matching procedure as done in Ref . [H] , namely the value of V^bare is fixed, the 
present result implies that the point does not move in a subspace of the parameters. This is 
likely the manifestation of a new fixed point in the Lorentz non-invariant formulation of the 
VM. Approaching the restoration point of chiral symmetry, the physical pion velocity itself 
would flow into the fixed point. 

We should distinguish the consequences within HLS/VM from those beyond HLS/VM. 
Clearly the determination of the definite value of the bare pion velocity is done outside 
HLS/VM. On the other hand, our main result (j5.72j) holds independently of the value of 
the bare pion velocity itself. Applying this result to the case where one starts from the bare 
HLS theory with Lorentz invariance, i.e., V^bare = 1, one finds that the pion velocity at T c 
becomes the speed of light since v n = V^bare = 1, as obtained in Ref. [TB*] . 

As a consequence of the relation (|5.72|) . we can determine the temporal and spatial pion 
decay constants at the critical temperature when we take the bare pion velocity as finite. 
In the following, we study these decay constants and discuss their determinations based on 
Eq. (j5.72j) . Using Eq. (J5.69)) . we solve the RGEs ()5.67)) and ()5.68)) and obtain a relation 
between two parametric pion decay constants as F*(0; T c )F*(0;T c ) = V*(F*(0; T)) 2 . From 
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this and (|5.70|) . the temporal and spatial pion decay constants with the quantum and hadronic 
corrections are obtained as 

Ul? = (^(0;r e )) 2 -^^7?, 

Nt 

fJl = Ft(0;T c )F°(0-T c ) - -^V n T 2 

= VMf- (5-73) 

Since the order parameter (/*/*) vanishes as expected at the critical temperature, we find that 
= V*(fly = . Multiplying both side by v\ = V%, the above expression is reduced to 

(/*) 2 = VM) 2 = 0. (5.74) 

Now, the spatial pion decay constant vanishes at the critical temperature, f*{T c ) = . In the 
case of a vanishing pion velocity, /* can be finite at the restoration point. On the other hand, 
when V n is finite, Eq. (15.74)) leads to /*(T C ) = 0. Thus we find that both temporal and spatial 
pion decay constants vanish simultaneously at the critical temperature when the bare pion 
velocity is determined as finite. 

In order to know the value of the (bare) pion velocity, we need to specify a method that 
determines the bare parameters. 

As we will show below, the analysis performed on the basis of a Wilsonian matching gives 
the finie bare pion velocity at the critical temperature, i.e., K-,barc ^ 0. Thus, by combining 
Eq. (15.72)1 with estimation of V^barei the value of the physical pion velocity v n (T) at the critical 
temperature is obtained to be finite )18j . 



Estimation 

Now we evaluate the physical pion velocity at the critical temperature T c starting from the 
Lorentz non-invariant bare Lagrangian. According to the non-renormalization theorem, the 
bare velocity so calculated should correspond to the physical pion velocity at the chiral tran- 
sition point. Now, in the VM, bare parameters are determined by matching the HLS to QCD 
at the matching scale A and at temperature T = T c . 

We begin with a summary of the pion velocity found in the HLS / VM theory with Lorentz 
invariance ^SJEIES]- The pion velocity is given by Eq. (J3.135)) : 

vl(P; T) = 1 + [Rem T) - Reflip; T)l . (5.75) 

The VM dictates that if one ignores Lorentz non-invariance in the bare Lagrangian in medium, 
the pion velocity approaches the speed of light as T — > T c [13 EE] . 
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In the following, we extend the matching condition valid at low temperature, Eq. ([4.31)1 . 
to near the critical temperature, and determine the bare pion velocity at T c . We should in 
principle evaluate the matrix elements in terms of QCD variables only in order for performing 
the Wilsonian matching, which is as yet unavailable from model-independent QCD calcula- 
tions. Therefore, we make an estimation by extending the dilute gas approximation adopted 
in the QCD sum-rule analysis in the low-temperature region to the critical temperature with 
including all the light degrees of freedom expected in the VM. In the HLS/VM theory, both 
the longitudinal and transverse p mesons become massless at the critical temperature since the 
HLS gauge coupling constant g vanishes. At the critical point, the longitudinal p meson which 
becomes the NG boson a couples to the vector current whereas the transverse p mesons decou- 
ple from the theory because of the vanishing g. Thus we assume that thermal fluctuations of 
the system are dominated near T c not only by the pions but also by the longitudinal p mesons. 
In evaluating the thermal matrix elements of the non-scalar operators in the OPE, we extend 
the thermal pion gas approximation employed in Ref. [73] to the longitudinal p mesons that 
figure in our approach jTH] - This is feasible since at the critical temperature, we expect the 
equality Al(T c ) = A1(T C ) to hold as the massless p meson is the chiral partner of the pion in the 
VM * 4 . It should be noted that, although we use the dilute gas approximation, the treatment 
here is already beyond the low-temperature approximation adopted in Eq. ()4.22j) because the 
contribution from p meson is negligible in the low-temperature region. Since we treat the pion 
as a massless particle in the present analysis, it is reasonable to take AJ(T) ~ A\{T = 0). We 
therefore use 

A P A (T) ~ Al(T) ~ Al{T = 0) for T ~T C . (5.76) 

The properties of the scalar operators giving rise to the condensates are fairly well under- 
stood at chiral restoration. We know that the quark condensate must be zero at the critical 
temperature. Furthermore the value of the gluon condensate at the phase transition is known 
from lattice calculations to be roughly half of the one in the free space jHHj- We therefore can 
use in what follows the following values at T = T c : 

(qq) T = 0, (—G 2 ) T ~ 0.006GeV 4 . (5.77) 

7T 

Including the contributions from both pions and massless p mesons, Eq. (j4.26J) can be expressed 
as 



r ,(OPE)i(i) _ 32 4 T 6 



M 



(5.78) 



#4 We observe from Refs. jHSEI that Al (u+d] '(/x = 2.4 GcV) - A p 4 [u+d) {^ = 2.4 GeV) even at zero tempera- 
ture. 
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Therefore from Eq. (j4.31|) . we obtain the deviation <5 bare as 



5bare-l V* haxe - ^ ^ ^ 



1 32 j\ T 

(5.79) 



This is the matching condition to be used for determining the value of the bare pion velocity 
near the critical temperature. 

To make a rough estimate of <5bare, we use Al^ u+d \fi = 1 GeV) = 0.255 This value is 
arrived at by following Appendix B of [IS]- An^ is defined by 

{^q^D^.D^q^) = H) n_1 (PMi-^ - traces)^ (jm) , (5.80) 



where 



A^(fM) = 2 / dxx^iqix^) + (-l) n q(x,v)} ■ (5.81) 



For any charge state of the pion (tt°, tt + , tt~ -), Al {u+d) (n = 2,4) can be written in terms of the 
nth moment of valence quark distribution V^(fi) and sea quark distribution S^(fi) [73], 

A^ +d \ f ,) = W:( f i) + 8S:( f i) , (5.82) 

where 

V? = [ fa" -1 /^^), 
Jo 

SI = [ dxx^s^ix,^). (5.83) 



Simple parameterizations of the valence distribution v n (x,fi) and the sea distribution s n (x,fi) 
can be found in Ref. [03] — see Ref. [OH] for the updated results - where the parton distri- 
butions in the pions are determined through the 7r-N Drell-Yan and direct photon produc- 
tion processes. With the leading-order parton distribution functions given in [HI], we obtain 
^w{u+d) _ g 255 a t ^ = 1 GeV [73]. For the purpose of comparison with the lattice QCD 
result [02] > we nee d to calculate the value at fi = 2.4 GeV; it comes out to be Al {u+d) = 0.18. 
The value A^ u+d * = 0.18 is slightly bigger than the ~ 0.13 calculated by the lattice QCD |92j . 
while it is a bit smaller than the ~ 0.22 of Ref. [OB!- Note that the value ~ 0.18 agrees with 



the one determined by lattice QCD |92j within quoted errors. Using Af u+d) = 0.972 and 



A4 = 0.255 and for the range of matching scale (A = 0.8 — 1.1 GeV), that of QCD scale 
(A QCD = 0.30 - 0.45 GeV) and critical temperature (T c = 0.15 - 0.20 GeV), we get 



-'bare 



(T c ) = 0.0061 - 0.29 , (5.84) 
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where the A dependence of A\ \ is ignored as it is expected to be suppressed by more than 
1/A 6 . Thus we find the bare pion velocity to be close to the speed of light: 

K,barc(T c ) = 0.83 - 0.99 . (5.85) 

Thanks to the non-renormalization theorem, i.e., tv^c) = V^barc^c); we arrive at the physical 
pion velocity at chiral restoration [THj : 

v n (T c ) = 0.83 - 0.99. (5.86) 

This is in contrast to the result obtained from the chiral theory j^Jl, where the relevant 
degree of freedom near T c is only the pion. Their result is that the pion velocity becomes zero 
for T —* T c . Therefore by experiment and lattice analysis, we may be able to distinguish which 
theory is appropriate to describe the chiral phase transition. 



5.5 Critical Behavior of Vector Meson Mass 

We showed that the vector meson pole mass m p (T) goes to zero in the limit T — > T c and the 
chiral symmetry is restored as the VM in the framework of the HLS theory. In this section, 
we study how m p (T) falls in the limit T — > T c . 

/"TIT C\ 

We expand the vector current correlator G v at the matching scale around T c taking the 
limits g(A;T) < 1, a(A; T) - 1 < 1 and M p (A;T)/A < 1. Then the quantity G A - G v in the 
HLS sector is obtained as 

Gf LS \A;T)-G { y LS) (A;T) 

. 92(AiT) (M^)) ! _( a(A;r) ^)W), (5 . 87) 

where we neglected the higher order coefficients. The same quantity using the OPE has the 
following form: 

gt e hat) - G r E \A-,T) = *w-o°y , (5 . 88) 

Here we assume that all the terms of Eq. ()5.87|) have same scaling behaviors near T c . From 
Eqs. ()5.87j) and ()5.88j) . we obtain the scaling behavior on the gauge coupling constant as 

g{A;T) ~ (qq) T . (5.89) 

As shown in section IS~21 the vector meson pole mass m p (T) near T c is 

ml{T)^g\T)(a{T)Ft + 5T 2 ), (5.90) 
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where 5T 2 denotes the hadronic thermal correction and gives a positive contribution [see 
Eq. (|5.23j) ]. Thus m p (T) vanishes as (qq)r- 

m p (T)~(qq) T . (5.91) 

If (qq)\ falls as f%(T) near T c , then the p pole mass rrit(T) vanishes as f%{T). In such a case 
the scaling property in the VM may be consistent with the Brown-Rho scaling m p (T)/m p (0) ~ 

urn/wo) ma. 



Chapter 6 

Chiral Doubling of Heavy-Light 
Mesons 



Based on the manifestation of chiral symmetry a la linear sigma model, it was predicted 

a decade ago [EH EH! that the mass splitting AM between the M(0~,l~) and M(0 + ,1 + ) 

mesons where A4 denotes a heavy-light meson consisting of heavy quark Q and light antiquark 

q should be of the size of the constituent quark mass. Recently, BaBar p££j . CLEO [35] and 

subsequently the Belle collaboration jSSJ, discovered new D mesons with Q = c, c being 

charm quark, which most likely have spin-party + and 1 + and the mass difference to the 

D(0~,1~) is in fair agreement with the prediction of [§Tl I98j. In Ref. [EH], it was proposed 

that the splitting of M. and M. mesons could carry direct information on the property of chiral 

symmetry at some critical density or temperature at which the symmetry is restored * 1 . In 

this chapter, we pick up this idea and make a first step in consolidating the proposal of Ref. [EE] 

following Ref. [33]. I n doing this, we shall take the reverse direction: Instead of starting with 

a Lagrangian defined in the chiral-symmetry broken phase and then driving the system to the 

chiral symmetry restoration point by an external disturbance, we will start from an assumed 

structure of chiral symmetry at its restoration point and then make a prediction as to what 

happens to the splitting in the broken phase. We find that the splitting is directly proportional 

to the light-quark condensate and comes out to be of the size of the constituent quark mass 

consistent with the prediction of Refs. |HZl I0HJ - We shall associate this result as giving a link 

between the assumed structure of the chiral restoration point and the broken phase. 

Our procedure is following: We assume that the heavy-light hadrons are described by a VM- 

^ tl In what follows, unless otherwise specified, we will refer to heavy-light mesons generically as D but the 
arguments should apply better to heavier-quark mesons. Numerical estimates will however be made solely for 
the (open charm) D mesons. 
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fixed point theory at the chiral restoration point generically denoted C x (critical temperature 
T c [Tl] or density n c [22] or number of flavors |21j ) and by introducing the simplest form 
of the VM breaking terms, we compute the mass splitting of the chiral doublers in matter-free 
space in terms of the quantities that figure in the QCD correlators # 2 . 



6.1 Heavy Quark Symmetry 

Like the approximate light quark flavor symmetry (chiral symmetry) in the light quark sector 
of QCD, the heavy quark symmetry is also approximate in the situation where the heavy quark 
mass m,Q is much larger than Aqcd ; Aqcd/^q <S 1- The heavy quark symmetry arises in the 
limit of QCD where the heavy quark, charm and bottom, masses tuq are taken to be infinity. 
In this limit, dynamics of the heavy quark Q is independent of its mass and spin |l()()j . In 
this section, we give a brief review of the heavy quark symmetry and the multiplet of hadrons 
containing a single heavy quark. 

In the limit where the heavy quark mass goes to infinity, rriQ — > oo, the four velocity of 
the heavy quark is fixed. The size of the meson is of order 1/Aqcd, and then the typical 
momentum scale of the light degrees of freedom is of order Aqcd- We decompose the four 
momentum of the heavy quark into 

Pn = m Q v^ + kp, (6.1) 

where is called the residual momentum and is of order Aqcd- The propagator of the heavy 
quark is given by 

-(p^ + m Q ) ^ 2 ^ 



p 2 — vtlq + ie 

Substituting Eq. (|6.1|) into the propagator, we obtain the leading form as follows: 

— 7~~~ ' (6-3) 

V ■ K + It 

where we moved the projection operator + l)/2 taking out of the large component. We 

should note that it is independent of the heavy quark mass vtlq and the gamma matrices 

are completely disappeared. Similarly the vertices are also independent of vtlq and gamma 

matrices. The heavy quark is sitting with infinite mass rriq — > oo, and thus the spin of the 

heavy quark does not flip at all by the interaction of the light degrees of freedom which is of 

order Aqcd- Thus the absence of gamma matrices indicates that the spin of heavy quark is 

^Introducing vector mesons in the light-quark sector of heavy-light mesons was considered in Ref. |97j but 
without the matching to QCD and hence without the VM fixed point. 
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conserved. Hence if there are Nf heavy quarks with the same four velocity, the effective heavy 
quark theory has a SU(2Nf) spin-flavor symmetry |lUUj . 

Let us consider the heavy-light meson multiplet. For the hadrons containing a single heavy 
quark Q, there is a SU (2) spin symmetry between Q with up-spin, Q(T) an d Q with down-spin, 
Q(l). Thus the hadrons with the spin combined f and j. with the spin of the light degrees of 
freedom sii ght are the spin partners each other. In general for each Slight, there is a degenerate 
doublet with the total spin s + and s_: 

s± = Slight ± - • (6.4) 

Hadron in the ground state has Sn g ht = 1/2 and negative parity and then the hadron with 
J p = 0~ is the spin partner of the one with J p = 1~. 

6.2 An Effective Field Theory 

In this section we give our reasoning that leads to the Lagrangian that defines our approach. 
Here we construct the Lagrangian using the approximate chiral S77(3)l x S77(3)r symmetry in 
the light-quark sector and the heavy quark symmetry in the heavy-quark sector. We will start 
from the Lagrangian given at the vector manifestation (VM) fixed point. We first describe 
how to construct the fixed point Lagrangian based on the VM. Then, we account for the effect 
of spontaneous chiral symmetry breaking by adding a bare parameter for the mass splitting in 
the heavy sector and including the deviation of the HLS parameters from the values at the VM 
fixed point. The explicit form of the Lagrangian so constructed is shown in subsection 1(12.81 

6.2.1 The fixed point Lagrangian 

Let us consider the VM at the point at which chiral symmetry is restored (in the chiral limit). 
At the VM at its fixed point characterized by (g, a) = (0, 1), the two 1-forms become 

"J* = ^7 • ££ - 0*6. ■ • (6-5) 

Note that the above ay and a±^ do not contain the HLS gauge field since the gauge coupling 
g vanishes at the VM fixed point. It is convenient to define the (L,R) 1-forms: 

oir^ = ay + a ±fl = -d^n ■ £r , 

fti/i = ay - a± fi = -<9 m £l • ^1 > (6-6) 
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which can be regarded as belonging to the chiral representation (1,8) and (8, 1), respectively, 
transforming under chiral S77(3)l x SU(3)n as 

«l m -* QhOiLnOi • (6.7) 
By using these 1-forms, the HLS Lagrangian at the VM fixed point can be written as 



1 I 

A* g ht = 2 F * tr \ a R^ a R\ + 2 F l tr [ a Lpa£] , (6. 



where the * affixed to the Lagrangian denotes that it is a fixed-point Lagrangian, and F w 
denotes the bare pion decay constant. Note that the physical pion decay constant f n vanishes 
at the VM fixed point by the quadratic divergence although the bare one is non-zero ^3]. It 
should be stressed that the above fixed point Lagrangian is approached only as a limit of chiral 
symmetry restoration [T3] . 

Next we consider the fixed-point Lagrangian of the heavy meson sector at the chiral restora- 
tion point identified with the VM fixed point. Let us introduce two heavy- meson fields T~Cr 
and Hl transforming under chiral SU(3)r x SU(3)l as 

n R ^H R gi, H L ^H L g{. (6.9) 

By using these fields together with the light-meson 1-forms a£ R , the fixed point Lagrangian 
of the heavy mesons is expressed as * 3 

^heavy = ~tr [HrW^Hr] - tl [H L tV^H L ] 

+ m tr [HrH r + H l H l ] 



+ 2k tr 



(6.10) 



where is the velocity of heavy meson, mo represents the mass generated by the interaction 
between heavy quark and the "pion cloud" surrounding the heavy quark, and k is a real 
constant to be determined. 



6.2.2 Effects of spontaneous chiral symmetry breaking 

Next we consider what happens in the broken phase of chiral symmetry. In the real world at low 

temperature and low density, chiral symmetry is spontaneously broken by the non-vanishing 

# 3 We assign the right chirality to TLr, and the left chirality to Hl- Then the interaction term has the right 
and left projection operators. Note that the insertion of (1 ± 75) to kinetic and mass termes does not cause 
any difference. 
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quark condensate. In the scenario of chiral-symmetry manifestation a la linear sigma model, 
the effect of spontaneous chiral symmetry breaking is expressed by the vacuum expectation 
value (VEV) of the scalar fields. In the VM, on the other hand, it is signaled by the HLS 
Lagrangian departing from the VM fixed point: There the gauge coupling constant g 7^ * 4 
and we have the kinetic term of the HLS gauge bosons £ p km = — |t r [PfiuP^]- The derivatives 
in the HLS 1-forms become the covariant derivatives and then a^^ and a Rfl are covariantized: 

% = - igpp, 

otR^ — > a Rfl = a Rtl - gp^, 

-> «l m = <y Lfl - gp^. (6.11) 

These 1-forms transform as a R ( L ^ —> haRrn^h* with h G [<S77(3)v]iocai as shown in Eq. (j3.6|) . 

Although a = 1 at the VM fixed point, generally a ^ 1 in the broken phase. We therefore 
expect to have a term of the form |(a — l)F^tr[aL^cx R }. Thus the Lagrangian for the light 
mesons takes the following form: 

flight = ° ^ Fl%r[a. Rtl a. R + a L ^ L ] 

+ ^^F*tr[a^a£] + C pkin . (6.12) 

By using aiu and a± fl given in Eq. (J3.5j) . this Lagrangian is rewritten as 

Aight = F^tr[a_L M a!^] + Fjtrfdi^d^] + £ pkin , (6.13) 

which is nothing but the general HLS Lagrangian. 

We next consider the spontaneous breaking of chiral symmetry in the heavy-meson sector. 
One of the most important effects of the symmetry breaking is to generate the mass split- 
ting between the odd parity multiplet and the even parity multiplet jHZI- This effect can be 
represented by the Lagrangian of the form: 

£ xSB = ^AMtr[H L H R + H R H L ] , (6.14) 

where Hr(l) transforms under the HLS as Hr(l) —> 7~Lr(l) . Here AM is the bare parameter 
corresponding to the mass splitting between the two multiplets. An important point of our 
work is that the bare AM can be determined by matching the EFT with QCD as we will show 
in section IB~5l The matching actually shows that AM is proportional to the quark condensate: 

AM~(qq). (6.15) 

# 4 Actually, near the chiral restoration point, the Wilsonian matching between HLS and QCD dictates [T3] 
that (in the chiral limit) the HLS gauge coupling be proportional to the quark condensate: g ~ (qq)- Sec 
section 1531 
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6.2.3 Lagrangian in parity eigenfields 



In order to compute the mass splitting between M. and Ai, it is convenient to go to the 
corresponding fields in parity eigenstate, H (odd-parity) and G (even-parity) as defined, e.g., 
in Ref. |JE|; 

1 

i 

Here, the pseudoscalar meson P and the vector meson P* are included in the H field as 

1 + V/" 



[G-iH l5 ] , 



(6.16) 



H 



[, 15 p + y p;] 



2 L - MJ ' 
and the scalar meson Q* and the axial-vector meson Q M are in G as 

1 + iyy" 



G 



[Q* - i-f-rsQA 



(6.17) 



(6.18) 



In terms of the H and G fields, the heavy-meson Lagrangian off the VM fixed point is of the 
form 

l~- heavy ^kin ^int > (6.19) 



with 



£int 



tr [H (ziyD" - M H ) H] - tr [G (iv^D" - M G ) G] 



(6.20) 



k 



tr[^ 75 Q^H] - tv[Hv^H] 

+ ix[Gw B a£G]+tr[Gv lt a§G] 

- iti[Ga L ^ lb H] + iti[Ha L ^ lb G} 

- itv[Ga\\^H] + itv[Ha\\^G\ 

where the covariant derivatives acting on H and G are defined as 

DpH = (dp - igpj H , D^G = {dp - igp^) G . 



(6.21) 



(6.22) 



In the above expression, Mh and Mq denote the masses of the parity-odd multiplet H and 
the parity-even multiplet G, respectively. They are related to m and AM as 



M H = -m --AM, 

Mq = -mo + ^AM . 
The mass splitting between G and H is therefore given by 

Mq - M H = AM . 



(6.23) 



(6.24) 
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6.3 Matching to the Operator Product Expansion 



The bare parameter AMb arc which carries information on QCD should be determined by match- 
ing the EFT correltators to QCD ones. We are concerned with the pseudoscalar correlator Gp 
and the scalar correlator G$- In the EFT sector, the correlators at the matching scale are of 
the form * 5 



G P (Q 2 ) 
G S (Q 2 ) 



F D M D 



Q 2 ' 



Ml 
Ml + Q 2 ' 



(6.25) 



where Fp> (F D ) denotes the D-meson (Z)-meson) decay constant and the space-like momentum 
Q 2 = {Mo + Ajv/) 2 with Am being the matching scale. We note that the heavy quark limit 
Mo — > oo should be taken with Am kept fixed since Am must be smaller than the chiral 
symmetry breaking scale characterized by A x ~ 4:irf n . Then, Q 2 should be regarded as Q 2 ~ 
M D in the present framework based on the chiral and heavy quark symmetries. If we ignore 
the difference between F D and F D which can be justified by the QCD sum rule analysis |lUlj . 
then we get 

A SP (Q 2 ) = G S (Q 2 ) - G P (Q 2 ) ~ —°-2-AM D . (6.26) 



In the QCD sector, the correlators G$ and Gp are given by the operator product expansion 
(OPE) as 



G S (Q 2 



G P (Q" 



+ 



G(Q 2 ) 


pert 


m 2 H 




m 2 H + Q 2 


G(Q 2 ) 


pert 


m 2 H 





m H (qq) 



12vr 



m 2 H + Q 2 



m H {qq) + -^(G^G^) 
lZtx 



(6.27) 



where rrtu is the heavy-quark mass. To the accuracy we are aiming at, the OPE can be 
truncated at 0(1/ Q 2 ). The explicit expression for the perturbative contribution G(Q 2 )| pcrt is 
available in the literature but we do not need it since it drops out in the difference. From these 
correlators, the A$p becomes 



A SP (Q S 



2m 3 H 



m 



-(??)■ 



(6.28) 



^ 5 Hcre and in the rest of this chapter, the heavy meson is denoted D with the open charm heavy meson in 
mind. However the arguments (except for the numerical values) are generic for all heavy mesons M. 
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Figure 6.1: Diagram contributing to the mass difference. 

Equating Eq. ()6.26j) to Eq. (j6.28j) and neglecting the difference (m# — Md), we obtain the 
following matching condition: 

3F|, AM ~ —2(qq). (6.29) 
Thus at the matching scale, the splitting is 

AM baJe ~-^M (6.30) 

As announced, the bare splitting is indeed proportional to the light-quark condensate. The 
quantum corrections do not change the dependence on the quark condensate [see section l6~^j . 



6.4 Quantum Corrections and RGE 

Given the bare Lagrangian whose parameters are fixed at the matching scale Am, the next step 
is to decimate the theory a la Wilson to the scale at which AM is measured. This amounts to 
calculating quantum corrections to the mass difference AM in the framework of the present 
EFT. 

This calculation turns out to be surprisingly simple for a ~ 1. If one sets a = 1 which 
is the approximation we are adopting here, «l does not mix with an in the light sector, and 
then ai couples to only Hi and or to only Hr. As a result TLl(r) cannot connect to Hr^l) 
by the exchange of or ocr. Only the p-loop links between the fields with different chiralities 
as shown in Fig. 16.11 We have verified this approximation to be reliable since corrections to 
the result with a = 1 come only at higher loop orders [see the next paragraph]. The diagram 
shown in Fig. 16.11 contributes to the two-point function as 

div 



n 



LB 



2 

; AMC 2 (iV / )|^(l-2A;-A ; 2 )lnA, (6.31) 



where C2(Nf) is the second Casimir defined by (T a )ij(T a )ji = C2(Nf)5n with i, j and / denoting 
the flavor indices of the light quarks. This divergence is renormalized by the bare contribution 
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of the form n^bare = |AMb are . Thus the renormalization-group equation (RGE) takes the 



2 

form 



A^^T = C 2 (N f ) jL(l-2k- k 2 ) AM. (6.32) 



For an approximate estimate that we are interested in at this point, it seems reasonable to 
ignore the scale dependence in g and k. Then the solution is simple: 

AM = AM bare x C quantum , (6.33) 

where we define C quantum by 

9 2 (, „, , a \, A" 



Cquantum exp 



-C 2 (N f )^(l-2k-k 2 )ln-\ . (6.34) 



//J 

This shows unequivocally that the mass splitting is dictated by the "bare" splitting AM bare 
proportional to (qq) corrected by the quantum effect C quantum . 

Next we lift the condition a = 1 made in the above analysis. For this purpose, we compute 
the quantum effects to the masses of 0~ (P) and + (Q*) -D-mesons by calculating the one-loop 
corrections to the two-point functions of P and Q* denoted by Upp and Hq*q* [for the explicit 
calculation, see Appendix |HJ. We find that amazingly, the resultant form of the quantum 
correction exactly agrees with the previous one which was obtained by taking a = 1. To arrive 
at this result, it is essential that P (or P*) be the chiral partner of Q* (or Q^) as follows: 
The loop diagrams shown in Fig. lH.l| and Fig. lH.2| in Appendix |H] have power and logarithmic 
divergences. However all the divergences of the diagrams with pion loop are exactly canceled 
among themselves since the internal (or external) particles are chiral partners. In a similar 
way, the exact cancellation takes place in the diagrams with a loop. Finally, the logarithmic 
divergence from the p loop does contribute to the mass difference. This shows that the effect 
of spontaneous chiral symmetry breaking introduced as the deviation of a from 1 does not get 
transferred to the heavy sector. Thus even in the case of a ^ 1, the bare mass splitting is 
enhanced by only the vector meson loop, with the pions not figuring in the quantum corrections 
one- loop order. Solving the RGE ([H.lip . which is exactly same as Eq. ()6.32|) . we 
obtain exactly the same mass splitting as the one given in Eq. (J6.33J) . 



6.5 Mass Splitting 

In this section we make a numerical estimation of the mass splitting for the chiral doublers in 
the open charm system. (Here D denotes the open charm meson.) Since we are considering the 
chiral limit, strictly speaking, a precise comparison with experiments is not feasible particularly 
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if the light quark is strange, so what we obtain should be considered as semi-quantitative at 
best. This caveat should be kept in mind in what follows. 

Determining the bare mass splitting from the matching condition ()6.30)1 requires the quark 
condensate at that scale and the D-meson decay constant Fd- For the quark condensate, we 
shall use the so-called "standard value" [1021 (QQ) = -(225 ±25 MeV) 3 at 1 GeV. Extrapolated 
to the scale Am = 1.1 GeV we shall adopt here, this gives 

(qq) AM = -(228 ± 25 MeV) 3 . (6.35) 

Unfortunately this value is not firmly established, there being no consensus on it. The values 
found in the literature vary widely, even by a factor of ~ 2, some higher |lU4j and some 
lower |l()5j . (We will study the dependence on the values of the mass splitting of those of the 
quark condensate later.) Here we take the standard value as a median * 6 . As for the D-meson 
decay constant, we take as a typical value Fd = 0.205 ± 0.020 GeV obtained from the QCD 
sum rule analysis |101j . Plugging the above input values into Eq. ()6.3U|) we obtain 

AAf barc ~ 0.19 GeV . (6.36) 

By taking /i = m p = 771 MeV, A = A M = 1.1 GeV, g = g(m p ) = 6.27 determined via the 
Wilsonian matching for (Am, Aqcd) = (1.1,0.4) GeV in Ref. ^Hj and k ~ 0.59 extracted from 
the D* — > Dtt decay [see section 16.6. lj in Eq. 1)6.34)1 . we find for Nf = 3 

Cquantum = 1-6 . (6.37) 

This is a sizable quantum correction involving only the vector meson. If one takes into account 
the uncertainties involved in the condensate and the decay constant, the quantum-corrected 
splitting AM comes out to be 

AM ps 0.31 ± 0.12 GeV . (6.38) 

Despite the uncertainty involved, ()6.38j) is a pleasing result. It shows that the splitting is 

indeed of the size of the constituent quark mass of a chiral quark £ ~ m p /3 ~ 310 MeV and 

is directly proportional to the quark condensate. 

was shown in Ref. |105| that there is a strong Nf dependence on the quark condensate and the value 
of the quark condensate for QCD with three masslcss quarks is smaller than the value used in estimating the 
value of the mass splitting in section Tfi. 51 In the present analysis, we extract the value of the coupling constant 
k from the experiment. To be consistent, we need to use the quark condensate together with other parameters 
involved determined at the same scale from experimental and/or lattice data. This corresponds to the standard 
value of the condensate we are using here. 
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We emphasized in the above analysis that there is a great deal of uncertainty on the value 
of the quark condensate at the relevant matching scale A^. Here we list a few examples to 
show what sort of uncertainty we are faced with. We took (<?<?) 1GcV = —(225 ±25 MeV) 3 in the 
above analysis as a "standard value." For comparison, we shall take two other values quoted 
in Ref. |lU4j (without making any judgments on their validity). Consider therefore 

(^W = -(225 ±25 MeV) 3 , 
/- \ / -(273 ±19 MeV) 3 , 

{qq) ^ = 1 -(316 ±24 MeV) 3 . ^ 

Brought by RGE to the scale we are working at, Am = 1.1 GeV, and substituted into our 
formula for AM, we get the corresponding quantum corrected splitting 

0.31 ±0.12 GeV , 

0.43 ±0.12 GeV , (6.40) 
0.67 ±0.20 GeV . 

This result clearly shows that the splitting cannot be pinned down unless one has a confirmed 
quark condensate. 

We should stress other several caveats associated with ()6.38j) . Apart from the sensitivity 
to the quark condensate, if one naively plugs in the matching scale A M into the RGE solution, 
one finds the splitting is not insensitive as it should be to the scale change. This is neither 
surprising nor too disturbing since our RGE solution is obtained with the scale dependence in 
both g and k ignored. In order to eliminate this dependence on the matching scale, it will be 
necessary to solve the RGE with the full scale dependence taken into account - which is at 
the moment beyond our scope here. The best we can do within the scheme adopted is to pick 
the optimal Am determined phenomenologically from elsewhere ^H] arid this is what we have 
done above. 



6.6 Hadronic Decay Modes 

In this section we turn to the hadronic decay processes of the D mesons and make predictions 
of our scenario based on the vector manifestation (VM) of chiral symmetry. Here we adopt 
the notations D u d and D U ( i for the heavy ground-state mesons and heavy excited mesons 
composed of cu and cd, and D s and D s for those composed of cs. The spin-parity quantum 
numbers will be explicitly written as _D Uj( i(0~). For the heavy vector meson, we follow the 
notation adopted by the Particle Data Group (PDG) [2E] and write _D* d (l~) and -D*(l~). 
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Unless otherwise noted, the masses of the ground-state heavy mesons will be denoted as 
and those of the excited states as M5. 

6.6.1 D* — > D -\- 7r 

Before studying the decay processes of the excited heavy mesons, we first calculate the decay 
width of D* d — > D Uj d + 7r so as to determine the coupling constant k. The decay widths of the 
7T° and the n ± modes are given by 

r(^(i-)-JM°-) + *°) = dss:(*>£ 

where p n = \p n \ denotes the three-momentum of the pion in the rest frame of the decaying 
particle D* d (l~), and Mq the "heavy quark mass" introduced for correctly normalizing the 
heavy meson field. In the present analysis we use the following reduced mass for definiteness: 

M Q = - (m d{0 -) + 3M D .(i-)) = 1974 MeV . (6.42) 

The total width is not determined for D*(l~), although the branching fractions for both 
the 7T° and the ir + decay modes are known experimentally. For D d (l~) meson, on the other 
hand, the total width is also determined. Using the values listed in PDG table j2H], the partial 
decay widths are estimated to be 

T(D*(1-) -> D d (0~) + 7r°) = 29.5 ± 6.8 keV , 

T(D* d (r) D u (0~) + 7T+) = 65 ± 15 keV . (6.43) 

Here the 7r° mode will be used as an input to fix k. From the experimental masses M^*^-) = 
2010.1 MeV, M Dd (p-) = 1869.4 MeV and M n o = 134.9766 MeV together with the value of the 
pion decay constant F n = 92.42 ± 0.26 MeV, we obtain 

k = 0.59 ± 0.07 . (6.44) 

Note that the error is mainly from that of the D d (l~) — > Dd(0~) + ir° decay width. 

In the following analysis, we shall use the central value of k to make predictions for the 
decay widths of D mesons. Each prediction includes at least about 25% error from the 
value of k. For the masses of excited D mesons, we use M^^^ = 2317 MeV determined by 
BaBar [S3, M D ~ s(1+) = 2460 MeV by CLEO j3Sj and (M^ d(0+) , M 6u = (2308, 2427) MeV 
by Belle jHSj- Table l6~D summarizes the input parameters used in the present analysis. 
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D U) d meson masses 


M D M+) 




M^j* J1-) 


Mn h (o-) 


(MeV) 


2427 


2308 


2010 


1865 


D s meson masses 


1V1 D 3 (1+) 


J "D S (0+) 






(MeV) 


2460 


2317 


2112 


1969 


Light meson masses 


M f 


p 


M„ 


Ms 


(MeV) 


138.039 


771.1 


547.30 


1019.456 


7r°-?7 mixing 


An 


A 2 i 


LV, (MeV) 2 






0.71 


-0.52 


-4.25 x 10 3 


-1.06 x 10" 2 


<p-p mixing 


iy^- (MeV) 
3.11 x 10" 4 


T p ^ +7T - (MeV) 
149.2 







Table 6.1: The values of input parameters. We use the values of M i =,^ 0+ - ) [33], M^^+j jSHl and 
d ^ Q+ 1+ n The .D mesons in the ground state, light mesons and decay widths T(0, p) 

are the values listed by the PDG table [2E]. As for the parameters associated with the 7r°-r] 
mixing, we use the values given in Refs. |1()6^ 1107] . 

6.6.2 D — > D -\- 7T 

For the systems of cu and cd, the following decay processes of the D u ^ meson into the D u ^ 
meson and one pion are allowed by the spin and parity: 

D u4 (0 + ) -> D u>d {0-) + 7T D u , d {l + ) -> D* u>d {r) + 7T. (6.45) 

Their partial decay widths are given by 

TVn n , M Pn f k M Q \ 2 

r(i>„,^ D „, J + .») = |(A^) 2 , (6.46) 

where E n is the energy of the pion, and the reduced mass Mq is defined in Eq. (|6.42|) . With 
the input parameters given in Table l6~T| these decay widths come out to be 

£> M (0-)+7r°)=74±18MeV, 
D u , d {0~) + tt ± ) = 147 ± 37 MeV , 
D l,A l ~) + 7r °) = 57 ± 14 MeV , 

D * u ,d( l ~) + 7r± ) = 114 ±29 MeV , (6.47) 
Here we have put the 25% error from the value of k, which we expect the dominant one. 



r(A, rf (o+) 
r(A^(i + ) 
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For the system of cs there are two decay processes of the D s meson into the D s meson and 
one pion: 

D s {0+) - D a (0-) + 7T° D S {1+) - D a *(l") + 7T°. (6.48) 

These processes violate the isospin invariance, and hence are suppressed. In the present analysis 
we assume as in Ref. |l()8j that the isospin violation occurs dominantly through the ir -r] mixing. 
In other words, we assume that the D s meson decays into the D s meson and the virtual r\ meson 
which mixes with the ir° through the 7i°-r] mixing. Then, the decay width is given by 

T(D S ^D s + n°) = ^- (A^^A^) 2 , (6.49) 
where A^-o^ denotes the ir°-ri mixing and takes the following form ll()7j : 

K ° r > = ' w^wj ^ ~ K ^ M - o) (6 - 50) 

with II^Ojj and K^„ being the mass-type and kinetic-type ir°-ri mixing, respectively. An and 
A 2 \ are the components of the 77-77' mixing matrix in the two-mixing-angle scheme |lU6j . By 
using the values listed in Table 16. 1\ the 7r°-77 mixing is estimated as 

= _5.3 x 1(T 3 . (6.51) 

From this value, the predicted decay widths are estimated as 

r( J D s (0 + )^ J D s (0') + 7r°)~4keV, 

T(D S (1 + ) -> D*(l~) + vr°) ~ 4keV . (6.52) 
6.6.3 D(l+) -> D(0 + ) + 7T 

With the masses of D U ^{1 + ) and D u( i(0 + ) listed in Table 16.11 the intra-multiplet decay 
D u ,d{^ + ) — ► D U( i(0 + ) + Ti is not allowed kinematically. Since the experimental errors for 
the masses are large * 7 , this decay mode may still turn out to be possible. To show how large 
the possible decay width is, we use / 0+ > = 2272 MeV and M/j n+) — 2464 MeV together 
with the formulas 



71 



pl ( 


Mq 


—) 


24tt \ 




Fir J 


Pl 


( M Q 


\ 2 

~) 


127T 







T(D u4 (l+)^D u , d (0 + )+7r + ) = ^- [—«-—) . (6.53) 



# 7 The Belle collaboration QUI EDI EU gives M bu = 2427 ± 26 ± 20 ± 17MeV and M bu , 0+) 
2308 ± 17 ± 15 ±28 MeV. 
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The resultant decay widths are given by 

T(D U)d (l+) -> D U)d (0+) + tt°) = 0.7 ± 0.2 MeV , 

T(D u4 (l + ) A^(0+) + 7T+) = 1.5 ± 0.4 MeV , (6.54) 

where we put 25% error coming from the value of k. We expect that this is the main source 
of the entire error. They are smaller by the order of 10~ 2 than other one-pion modes [see 
Table IB~2"] . This is caused by the suppression from the phase space. 

With the present input values of D masses, the process D s (l + ) — > D s (0 + ) +tt° is kinemat- 
ically allowed. Similarly to the D s — > D s + 7r° decay, we assume that this decay is dominated 
by the process through the n°-ri mixing. Then, the decay width is estimated as 

rp,(i+)^(o+) + .°) = |(A_^.^) 2 

~ 2 x 10~ 3 keV. (6.55) 
This is very tiny due to the isospin violation and the phase-space suppression. 

6.6.4 D^D + 2tt 

There are several processes such as D — ► D + 7r ± 7r =F to which the light scalar mesons could give 

important contributions. In models based on the standard scenario of the chiral symmetry 

restoration in the light quark sector, the scalar-meson coupling to the heavy-quark system is 

related to the pion coupling, enabling one to compute the decay width. In our model based on 

the VM of the chiral symmetry restoration, on the other hand, it is the coupling constant of the 

vector meson to the heavy system that is related to the pion coupling constant: Here coupling 

of the scalar meson is not directly connected, at least in the present framework which contains 

no explicit scalar fields * 8 , to do that of the pion. So, while we cannot make firm predictions 

to processes for which scalar mesons might contribute, we can make definite predictions on 

certain decay widths for which scalar mesons do not figure. If one ignores isospin violation, 

the two-pion decay processes D u d — > D u4 + 7r ± 7r° receive no contributions from scalar mesons. 

We give predictions for these processes below. As for the two-pion decay modes of the D s 

meson, the scalar mesons could give a contribution. To have an idea, we shall also compute 

the vector-meson contribution to this process. 

# 8 Scalar excitations can of course be generated at high loop level to assure unitarity or with the account of 
QCD trace anomaly but we shall not attempt this extension in this paper. 
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First, consider D u ,d{® + ) — > D* d (l ) + 7i ± 'K . In this process, there are two contributions: 

D u , d (0 + ) D: i4 {l-) + ^ (direct) 

D l,A l ~) + (P ± -> ^tt ) (p-mediation) . (6.56) 

The decay width is given by 

r(D u>d (o + ) ^ D* u4 (n + 

Ml k- 



64(2vr) 3 M| F£ 



dm ln I dm L\ F DD 



mL - 4M 2 + -^(m^ _Ml -M 2 n - 2Ml + 2m 2 D J 

7T7T 7T 1 /j 71^2 ^ Jj L) 7T ' U7t / 



(6.57) 



with m 2 ^ = (p D + p n ) 2 and m 2 ^ = (p l7r + P2tt) 2 - The form factor Fp D is taken to be of the 
form 

^ = 1 + < 6 - 58 > 

The first term of the form factor comes from the direct contribution and the second from 
the p-mediation. Here we have neglected the p meson width in the propagator, since the 
maximum value of m nn is about 300 MeV with the input values listed in Table I6~T1 We can see 
that the form factor Ff) D vanishes in the limit of — > 0, which is a consequence of chiral 
symmetry * 9 . We note that fn n7T \ mSix ~ 300 MeV makes this decay width strongly suppressed 
due to the large cancellation between the direct and p-mediated contributions. Furthermore, 
since 300 MeV is close to the two-pion threshold, additional suppression comes from the phase 
space. Due to these two types of suppressions the predicted decay width is predicted to be 
very small, of the order of 10~ 2 keV. * 10 

Next we consider the process D u4 {l + ) D* d (l~) +mv. Again there are two contributions, 
direct and a p-mediated: 

D u4 {\ + ) D;,ir)+i¥ (direct) 

Dldi 1 ') + (P ± ^n ) (p-mediation) (6.59) 



# 9 It should be stressed that this cancellation occurs because the vector meson is included consistently with 
chiral symmetry, and that it is not a specific feature of the VM. The chiral symmetry restoration based on the 
VM implies that the coupling constant of the vector meson to the heavy system is equal to that of the pion. 
# 10 Note that the prediction on the decay width is very sensitive to the precise value of the mass of D(0 + ) 
meson. 
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The resultant decay width is given by 



M 2 



k 2 



96(27r) 3 M| 



dm 2 D% 



d™*n\ F DD\ 



4Mi + 



(mL -Ml -Ml 



Dn) 



(6.60) 



Similarly to T(D u ^(0 + ) — > _D* d (l~) + ^ir ), the width is again suppressed due to the large 
cancellation between the direct and p-mediated contributions. The suppression from the phase 
space, on the other hand, is not so large since mnn\ maiX — 420 MeV is not so close to the two-pion 
threshold. The resulting decay width is 



+ 7T 71 J 



12 ± 3 keV 



(6.61) 



where the error comes from the value of k. 
The decay width of the process 



D, 



D u4 (0-) + (p ± - 



(direct) 
(p-mediation) 



(6.62) 



is given by 



r(D 



u,d\ 



D u , d (0 
k 2 



M 2 

192(27r)3M?^ I 



dm 2 \F~ I 2 

a " l -KTv\ r DD\ 



m 2 n7T - 4M 2 + jr^{mL ~ M% - M 2 D - 2M 2 + 2m Dn 
b 



2 \2 



(6.63) 



In the present case, max ~ 560 MeV is much larger than the two-pion threshold and hence 
the width becomes larger than other two-pion processes. We find 



r(D M (l + ) -> D Uid (0-) + t^tt ) = 310 ± 80keV 



(6.64) 



where the error comes from the value of k. 

Finally we turn to the decay D s (l + ) — > D s (0~) + 7r + 7r~ which as mentioned could receive 
direct contributions from scalar excitations. Since we have not incorporated scalar degrees of 
freedom in the theory, we might not be able to make a reliable estimate even if were to go to 
higher-loop orders. Just to have an idea as to how important the vector meson contribution 
can be, we calculate the decay width in which the D s meson decays into two pions through 
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the (j> meson. This isospin violating decay can occur through the direct (p-n-n coupling and 
the (f)-p mixing: 



D s (l + ) 



D s (0- 
D s (0- 



7T + 7T~ 



) + (0 - P 



-7T+7T- 



(direct) 
(j)-p mixing) 



(6.65) 



Since the main contribution to the </> — > ixtx is expected to be given by the <p-p mixing, we shall 
neglect the direct 0-7r-7r-coupling contribution in the following. Then the decay width is given 
by 



r(£> fl (l+) ^£> s ((T)+7r + 7r- 
k 2 



Ml 



192(2tt) 3 M| F% 



mi - AMI + 



dm Dn / dm^ 



M 2 p Tl 4 



L(mL-M 2 )(m^-M2)J 



Ml -M 2 D - 2Ml + 2m 2 D J 



(6.66) 



where li^p denotes the 0-p mixing given by 



2\2 



Ui p = (Mr - Mi) 



7T 7T 



(6.67) 



with p v (X) being the three-momentum of pion in the rest frame of the decaying particle 
X = (ft, p. Using the values listed in Table 16.11 we have 



= 530 (MeV) 2 



so the decay width is predicted to be 



r(D,(l+) - d s (q- 



7T+7T-) ~ 2 x 10" 4 keV . 



(6.68) 



(6.69) 



The (f)-p mixing is caused by the isospin violation, and this process is highly suppressed. We 
conclude that should a measured width come out to be substantially greater than what we 
found here, it would mean that either scalars must figure importantly or the VM is invalid in 
its present form. 



6.6.5 Summary of hadronic decay modes 

Our predictions of the decay widths are summarized in Table 16.21 It should be stressed that 
the values obtained in this paper on the one-pion reflect only that the D meson is the chiral 
partner of the D meson. They are not specific to the VM. We therefore expect that as 
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Decaying particle Process Width (MeV) 

D Ujd + ^0- + 7T° (7.4 ± 1.8) x 10 1 

(1.5 ±0.4) x 10 2 
~ 2 x 10~ 5 
(5.7 ± 1.4) x 10 1 
(1.1 ±0.3) x 10 2 
(3.1 ±0.8) x 10- 1 
(1.2 ±0.3) x 10~ 2 
~4x 10~ 3 
~2x 10~ 6 
~4x 10~ 3 

through — > p° — »• 7r + 7r _ ) ~2x 10 -7 



+ - 


-> - 


i 
h 7T U 


A + 

+ - 


-> - 




+ - 


-> 1 - 


± 
h 7T 7T 


1 + 
1 


— » i 


r 7T 


1 + 
1 T - 


-i — 
-> 1 - 


± 

h 7T 


1 + 


— > u - 


h 71 7T 


1 + 
1 


— > 1 




0+ - 


-> o- - 


h7T° 


1+ - 


^0+- 


h7T° 


1+ - 


-> 1" - 


hTT 


1+ - 


-> o- - 


h 7T + 7T~ 



Table 6.2: The predicted values of the hadronic decay processes. 25% error comes from the 
value of k. "~" implies that the precise value also depends on parameters other than k, and 
we can give only the order estimation. 

far as the one-pion processes are concerned, there will be no essential differences between our 
predictions and those in Ref. |108j . However, in the two-pion decay processes in which the scalar 
meson does not mediate, our scenario based on the VM can make definite predictions which 
might be distinguished from that based on the standard picture. Especially for D u> d{l + ) — > 
D u ,d{0~) ± ^^7?°, we obtain a larger width than for other two-pion modes. Although the 
predicted width is still small - perhaps too small to be detected experimentally, it is important 
because of the following reason. In our approach, since the excited heavy meson multiplets of 
-D(0 + ) and -D(l + ) denoted by G are the chiral partners to the ground-state multiplets denote 
by H, the G-H-tt coupling is the same as the H-H-tt coupling [see the fifth and first terms of 
Eq. ([6.21)1 ] . Thus the width which is dependent on the strength of k is a good probe to test 
our scenario. The common k is also essential for the ratio of the widths of the two-pion modes 
to those of the one-pion modes, which has no k dependence. These are therefore are definite 
predictions of our scenario. From the values listed in Table l6~2"j we obtain 



-p(had) 

1 7T+27T 



where r^f^l is the sum of the widths of the one-pion and two-pion modes of the decaying 
D U J1 + ). 



Chapter 7 

Summary and Discussions 



In this thesis, we first showed the details of the hidden local symmetry (HLS) theory at finite 
temperature. Then in order to determine the parameters of the theory, we presented the 
extension of the Wilsonian matching at zero temperature to the version in hot matter. Next 
we showed the formulation of the vector manifestation (VM) in hot matter using the HLS 
theory. Based on the VM, we gave several predictions which are testable by experiments and 
lattice analysis. In the following, we summarize them and give discussions. 

HLS theory at finite temperature 

In chapter El we applied the HLS theory developed in Ref. ^3] to QCD at finite temperature. 
Both pions and vector mesons are introduced as the relevant degrees of freedom into the theory, 
and thus the HLS theory has a wider range of the validity than the ordinary chiral perturbation 
theory (ChPT) including only pions. Especially we can see it in studying the pion decay 
constants which was carried out in section EOl Two kinds of pion decay constants, temporal 
and spatial components (/* and appear due to the lack of Lorentz invariance in medium 
and in general they do not agree with each other. The splitting between j\ and appears at 
one-loop level in the HLS theory, which is generated by the thermal loop contribution of the 
vector meson. [The non-zero splitting also implies that the pion velocity is not equal to the 
speed of light.] On the other hand, the splitting of the pion decay constant is not generated 
at one-loop level in the ordinary ChPT. 

Further by using the Wilsonian matching studied in chapter the intrinsic temperature 
dependences of the bare parameters are included. We in fact determined the intrinsic thermal 
effects to the bare pion decay constants by performing the matching to the operator product 
expansion (OPE) in section POl 
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Formulation of the VM in hot matter 

In chapter we first presented the formulation of the VM in hot matter using the HLS theory 
combined with the Wilsonian matching. Throughout this thesis, we assume that the chiral 
phase transition is of second or weakly first order so that the axial-vector and vector current 
correlators agree with each other at the critical temperature. Thus the Wilsonian matching 
leads to the following conditions on the parameters of the bare HLS Lagrangian: 

3barc(T) r ^ Tc (or M pMm (T) T ^ c 0) , 

a b are(T) = (F fT , barc (T)/F^ barc (T)) 2 T => Tc 1 . (7.1) 

We should stress the following features: (i) The above conditions for the bare parameters are 
realized by the intrinsic thermal effect only, (ii) (g, a) = (0, 1) is a fixed point of RGEs. The 
agreement of both current correlators is satisfied only in the case that we take into account 
the intrinsic thermal effects as well as the hadronic ones. 

Vector meson mass : Near T c , the vector meson mass has a positive thermal contribution, 
which is proportional to the square of the HLS gauge coupling g, as studied in section l5~2l 
When we take the limit T — > T c , both the parametric mass M p and hadronic thermal 
correction go to zero following Eq. (|7.1|) . Thus we showed that the vector meson mass 
vanishes: 



m%T) g\T) a(T)F^ = 0; T) + 6T 



T ^0 



where 5T 2 denotes the thermal contribution and F n (fi = 0; T) is the on-shell parametric 
pion decay constant obtained in Eq. (|5.H()jl to be non-zero. 

Order parameter : If we write explicitly the quantum and hadronic thermal corrections to 
the parameter F n , the pion decay constant at T c takes the following form as studied in 
section 15. HI 

H(T) T =T' F|(A; T.) - - ^ = . 

This implies that the order parameter goes to zero by the cancellation among the non- 
zero bare pion decay constant, the quantum correction and the leading contribution of 
matter effect to F%. Here we should note that both the above quantum and hadronic 
contributions include the contribution of the massless vector meson in the same footing 
as the pion. 



CHAPTER 7. SUMMARY AND DISCUSSIONS 



103 



The conditions for the bare parameters (|7.1|) . which we call "the VM conditions in hot matter", 
are essential for the VM to take place in hot matter. Further even when we take into account 
the Lorentz non-invariance in the bare HLS theory, we obtained the similar conditions to 
Eq. (|7.1|) and showed that the conditions are still non-renormalized: There are no quantum 
corrections to the conditions. In section 1531 we studied the critical behavior of m p and showed 
that m p (T) falls as the quark condensate in the limit T — > T c : 

m P ( T ) ~ (qq)t ■ 

At present, there are no clear lattice data for the p pole mass in hot matter. Our result 
here will be checked by lattice analyses in future. |112j 

Predictions 

In section 15.41 we gave the following predictions of the VM in hot matter: 

• Critical temperature : From the order parameter f n including the intrinsic thermal ef- 
fect as well as the hadronic correction, we obtained T c = 200 - 250 MeV for several choices 
of the matching scale (A = 0.8- 1.1 GeV) and the scale of QCD (A QCD = 0.30-0.45 GeV). 

• Axial- vector and vector charge susceptibilities (ASUS and VSUS): The masslss 
vector meson becomes a relevant degree of freedom near T c and it contributes to the 
VSUS. The equality between ASUS and VSUS is predicted by the VM. 

• Violation of vector dominance : We showed that the vector dominance (VD) of 
the electromagnetic form factor of the pion is largely violated near T c . This indicates 
that the assumption of the VD may need to be weakened, at least in some amounts, for 
consistently including the effect of the dropping mass of the vector meson. 

• Pion velocity : We first showed that the pion velocity at T c never receive any quantum 
and hadronic thermal corrections protected by the VM. Using the non-renormalization 
property, we next performed the Wilsonian matching to the OPE and then obtained that 
the pion velocity is close to the speed of light. 

This is in contrast to the result obtained from the chiral theory [HE] , where the relevant 
degree of freedom near T c is only the pion. Their result is that the pion velocity becomes 
zero for T — > T c . Therefore by experiment and lattice analysis, we may be able to 
distinguish which theory is appropriate to describe the chiral phase transition. 
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We note that the above estimated values may be changed by including higher order effects. 
On the other hand, it is expected that the VM is governed by the fixed point and not changed 
by such higher order effects. 

Remnant of the VM in the real world 

In chapter El we studied the heavy meson system which consists of one heavy quark and one 
light anti-quark, which was motivated by the recent discovery of new D mesons ( J p = + and 
1 + ) in Babar, CLEO and Belle [HUES EH]- We showed that the mass splitting between new 
D mesons (denoted by D in the text) and the existing D mesons ( J p = 0~ and 1~) is directly 
proportional to the light quark condensate. The estimated value was in good agreement with 
the experiments. We also gave the predictions on the decay widths and the branching ratios 
for the characteristic processes. In the two-pion decay processes in which the scalar meson 
does not mediate, our scenario gives definite predictions, since the vector meson coupling to 
the heavy system is equivalent to the pion coupling due to the VM. Although the predicted 
values of widths are small, we hope that they are clarified in future experiment. 

One of the significant results of the analysis, is that the vector meson plays an important 
role in accounting for the splitting in the D and D mesons: The bare mass splitting determined 
through the matching is estimated as about 190 MeV, too small to explain the observed mass 
difference. However by including the quantum corrections through the hadronic loop, the 
bare mass splitting is enhanced by ~ 60%, where only the loop effect of the vector meson 
contributes to the running of the mass splitting. The contributions from the pion loop are 
completely canceled among themselves. This implies that the observed mass difference can 
not be understood if one takes only the pion as the relevant degree of freedom and that we 
need other degrees of freedom. In the VM, it is nothing but the vector meson. The situation 
here is much like in the calculation of pion velocity at the chiral restoration point: The pion 
velocity is zero if the pion is the only effective degree of freedom but approaches 1 if the vector 
meson with the VM is included [TBI ITS] . 

In the construction of the EFT, we assumed the VM as the chiral symmetry restoration 
of QCD in the light sector, and then we introduced the effects of the spontaneous chiral 
symmetry breaking by putting the mass term as soft breaking term. In the analysis, we 
assumed that the coefficients of the interaction such as DDn and DDp denoted by k are 
universal, which are the remnant of the VM. Strictly speaking, each interaction term may 
have its own coefficient different from others. However, we expect that the effect of these 
interaction terms is suppressed by the factor 1/A and as a result the contribution to AM is 
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small since the dimension of them is higher than that of the mass term. 

Furthermore, the result is independent of the deviation of a from the fixed point value 1, 
one-loop level. This implies that the deviation of a from 1 which reflects the effect 
of spontaneous chiral symmetry breaking in the light quark sector does not get transmitted 
to the heavy sector. This strongly suggests that the deviation from a = 1 involves physics 
that is not as primary as the non-vanishing gauge coupling g ^ in the description of the 
broken phase: The deviation seems to be a "secondary" phenomenon, which is generated from 
g ^ as expected in Refs. [52j . [Although a = 1 is the fixed point of the RGE at one-loop 
level, the deviation of a from 1 is generated by the finite renormalization part once we allow 
the deviation of the gauge coupling g from (see Appendix E}.] In fact, even when we start 
from the bare HLS theory with g ^ and a = 1, the physical quantities obtained through 
the Wilsonian matching are in good agreement with experimental results as discussed in 
This observation supports the above argument. 

The consistency of the physics taking abare = 1 and the universal coupling k with the one in 
the real world remarkably indicates that the VM is realized as the chiral symmetry restoration 
in QCD. 

Several comments are in order: 

Although we concentrated on the hot matter calculation in this thesis, the present approach 
can be applied to the general hot and/or dense matter calculation. As mentioned in chapter 
the vector meson is one of the most important degrees of freedom in hot and/or dense matter 
since it will be lighter in medium than in the vacuum due to the (partial) chiral symmetry 
restoration. The HLS theory is an efficient implement to study hadron properties in hot and/or 
dense matter as well as the chiral phase transition. 

One might think that the VM is same as the Georgi's vector realization [02], in which the 
order parameter f n is non-zero although the chiral symmetry is unbroken. However in the VM 
the order parameter certainly becomes zero and the chiral symmetry is restored accompanied 
by massless vector meson. Therefore the VM is consistent with the Ward-Takahashi identity 
since it is the Wigner realization |11H| . 

As shown in Ref. [21], in the VM only the longitudinal p couples to the vector current near 
the critical point, and the transverse p is decoupled from it. The A± in the VM is resolved 
and/or decoupled from the axial- vector current near T c since there is no contribution in the 
vector current correlator to be matched with the axial-vector correlator. We expect that the 
scalar meson is also resolved and/or decoupled near T c since it in the VM is in the same 
representation as the Ai is in. We also expect that excited mesons are also resolved and/or 
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decoupled. 

The VSUS was obtained to be finite, which will imply that the screening mass of the vector 
meson is also finite at T c . 

In this thesis, we performed our analysis at the chiral limit. We need to include the explicit 
chiral symmetry breaking effect from the current quark masses when we apply the present 
analysis to the real QCD. In such a case, we need the Wilsonian matching conditions with 
including non-zero quark mass which have not yet been established. Here we expect that the 
qualitative structure obtained in the present analysis will not be changed by the inclusion of 
the current quark masses. 



Acknowledgment 



The author is sincerely grateful to Professor Masayasu Harada, who is her 
adviser, and Professor Mannque Rho for plenty of useful discussions and continuous 
encouragements throughout her research. She would like to thank Professor Koichi 
Yamawaki for valuable discussions and encouragement. She would like to express her 
thanks to Doctor Youngman Kim, who is a collaborator, for many fruitful discussions 
and exchange of their opinions. 

She would like to express her sincere gratitude for all the helps and encourage- 
ments she has received. She would never complete this thesis without their helps. 

Support by the 21st Century COE Program of Nagoya University provided by 
Japan Society for the Promotion of Science (15COEG01) is much appreciated. 



107 



Appendix A 

Polarization Tensors at Finite 
Temperature 



At non-zero temperature there exist four independent polarization tensors, u^u v ', {g^ u — u^u v \ 
P£ u and P£ v ■ The rest frame of medium is shown by w M = (1,0). Pf* v and P^ v are given 
by H3 



p fx PiPj\ 

t^Tiiv - 9fj,i\0ij ~pT) g i v 



Pl^v = —\9ijv — ~!r ) ~~ Pthv, (a.i) 

where we define p M = {po , P) and p = \p\. They have the following properties: 

Pi paf p f p p<w p f 

Px^^r = 0. (A.2) 
Let us decompose a tensor IF"'(po,p) into 

if" = uvn* + {gr - u»u v )w + pj^u L + p^ u u t . (a.3) 

Each component is obtained as 

jjt = n oo + ^ n i 0; 

Po 

p p p 2 
n L -n s = ^n^ + ^ir°, 

P P Po 

it t -it = \p t ^ v . (A.4) 
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Appendix B 

Loop Integrals at Finite Temperature 



In this appendix we list the explicit forms of the functions appearing in the hadronic thermal 
corrections, A , B and B^ [see Eqs. (I3.68|) - (J3.73|) for definitions] in various limits relevant 
to the present analysis. 

The functions A , which are independent of external momentum p p , is given by 



A (M p ;T) = Jl{M p] T) 
3 (0;T) = J 2 (T) , 



(B.l) 
(B.2) 



where J\ and 1 2 are defined in Eqs. (jD.ljl and (|D.2jl . 

We list the relevant limits of the functions B (p ,p; M p , 0; T) and B^ u (p , p; M p , 0; T) which 
appear in the two-point function of A p . When the pion momentum is taken as its on-shell, 
the function Bq becomes 

B (Po =P + i^P;M p ,0;T) 
d 3 k 



(2tt)3 [2u p e-p/ T -l{(uj p - p) 2 - {uiy - 2te(u p - p) (u p + pf - (u%) 2 + 2ie(u p + p) 



+ 



+ 



2ujI e^l T - 1 1 (u£ - p) 2 - lu 2 - 2ie(w? - p) (c£ + p) 2 - u 2 + 2ie(u% + p) 



(B.3) 



where we put e — > +0 to make the analytic continuation of the frequency po = i2irnT to the 
Minkowski variable, and we defined 



k\ 2 + Ml , (Jf 



\k — p\ 



(B.4) 



Two components B t and B s of B^ in the same limit are given by 



B (p =p + ie,p;M p ,0;T) 
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d 3 k r-i 



(2uj p - pf 



+ 



(2u p +pf 



(2tt) 3 [2u p e^' T -\\{u p - p) 2 - (c^) 2 - 2ie(u p - p) {u p + p) 2 - (ul) 2 + 2ie(u p + p) 



2u) p -p 



2ui p + p 



+ 



p • (2k — p) 

p 

-1 1 



(cj p - p) 2 - {ul) 2 - 2ie(u p - p) (cj p + p) 2 - (cj?) 2 + 2ie(cu p + p) 
(2^-p) 2 (2^+p) 2 



+ 



+ 



p • (2/c — p) 



P 



2<-p 



2^+p 



'wj? - p) 2 - cj 2 - 2ze(cj£ - p) (cj£ + p) 2 - uj 2 + 2ie(u;£ + p) 



(B.5) 



5 S (p =p + *e,p;M p ,0;T) 



(27T) 



2ui n e"e/ T - 1 



(2k-p-p 2 f 
p ■ {2k — p) 



+ 



1 



P 



{uj p - p) 2 - {ul) 2 - 2ie(u p - p) (uj p + p) 2 - (ul) 2 + 2ie(cu p + p) 
2w„ -p 



^ 2 - [ooiy 

2u p + p 



(u p - p) 2 - (ul) 2 - 2ie(u p - p) (uj p + p) 2 - (ul) 2 + 2ze(cj p + p) 
1 



2^e^/ T -l 
(2k-p-p 2 ) 2 



p 2 



+ 



p • (2A; — p) 



P 



(u% - p) 2 - u; 2 - 2%e{ul - p) (u% + p) 2 - u; 2 + 2ie{ul + p) 
2cu£ - p 2ujI + p 



(cu? - p) 2 - cj 2 - 2ie(cu? - p) (u;£ + p) 2 - uj 2 + 2ie(w? + p) 



(B.6) 



We further take the M p — »• limits of the above expressions. The limit of i?o(po = P + 
ie,p; M p ,0;T) includes the infrared logarithmic divergence InM 2 . However, it appears multi- 
plied by M 2 in lTj_ and IT|_, and the product M 2 B (p = p+ie,p; M p , 0; T) vanishes at M p — > 
limit: 



lim M 2 p B (p =p + ie, p; M p , 0; T) = 



(B.7) 



As for £> ' (po = P + ze, p; M p , 0; T), we obtain 



lim B (p = p + z'e, p; M p , 0; T) 

A/ p -»0 



-2/ 2 (T) , 

\imW(p =p + i6,P;M p ,0-,T) = -2/ 2 (T) . (B.8) 
In the static limit (p — > 0), the functions M 2 B (p Q ,p; M p , 0; T) and _B*(p ,p; M p , 0; T) — 



i? (po, p; M p , 0; T) become 



lim, Jl^So(po,p;M p ,0;T) 



Po^O 
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Ml 



d 3 k 



-1 



(2vr) 3 u* - (u p p ) 2 [uj p e<V T - 1 u p p e^/ T - 1 



lim 5 (p ,p-M p ,M p ;T) 

Po^O 

d 3 k 



1 



Wp- (2k -p) [uj p e^/ T -I 
1 1 

"^e-?/r_i 



4u; p - p- (2k -p) 



4(u p ) 2 + p ■ (2k - p) 



Taking the low momentum limit (p — * 0), these expressions become 

hmM 2 B (p = 0,p; M p , M p ; T) = \ \j°_ x (M p ; T) - J 2 (M p , T) 



lim B L (p = Q,p; M p , M p ; T) = - 2 

p— >o 

In the VM limit (M p — > 0), these are reduced to 



M 2 J°(M P ,T) + 2J 2 (M P ,T) 



limM 2 B (p = 0,p;M p ,M p ,T) ^ 0, 



lim B (p = 0,p;M p ,M p ;T) 
p— »o 



-47 2 (T). 



(B.17) 



(B.18) 



(B.19) 
(B.20) 

(B.21) 
(B.22) 



l s <j> s 

On the other hand, the functions Bq, B — B and B — B in the low momentum limit 



(p — > 0) are given by 



1 ~- 



lim S (po,p; M p , M p ; T) = -F 3 2 (p ; M p ; T) 



lim 



B (p ,p; M p , M p , T)-B ( Po ,p; M p , M p , T) 
B T (p ,p; M p , M p - T) - B s (p ,p; M„, M pi T) 

where the function F£ is defined in Appendix iDl 



lim 

P^oL 



(B.23) 

p*(p ;M p ,T), (B.24) 
pi(p ;M p ,T), (B.25) 



Appendix C 
Quantum Corrections 



In this appendix we briefly summarize the quantum corrections to the two-point functions. 
Let us define the functions £?Q Vac \ i^( vac )M l/ an d by the following integrals [T3] : 

R w^,„. M M i f d " k {2k- P n2k- P r 

B (p,M u M 2 ) - j i(2x) 4 [M? _ A2][M |_ (A _ p)2] ■ ( &3 > 

In the present analysis it is important to include the quadratic divergences to obtain the 
RGEs in the Wilsonian sense. Here, following Refs. [03 EH CSj we adopt the dimensional 
regularization and identify the quadratic divergences with the presence of poles of ultraviolet 
origin at n = 2 (HHl- This can be done by the following replacement in the Feynman integrals: 



d n k 1 A 2 f d n k k^k v A 2 



9,u ■ (C.4) 



i(2vr)" -k 2 (4vr) 2 ' J i(2n) n {-k 2 } 2 2(4vr 
On the other hand, the logarithmic divergence is identified with the pole at n = 4: 

^ + l^lnA 2 , (C.5) 

where 

- = j?— -1e + ln(47r) , (C.6) 
e 4 — n 

with 7s being the Euler constant. 

By using the replacements in Eqs. (|C.4J) and (|C.5jl . the integrals in Eqs. ()C.2J) and 

3l) are evaluated as 

A 2 M 2 A 2 
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Bt c) (p 2 ;M u M 2 ) 

B {vac) ^(p;M l ,M 2 ) 
1 



[In A 2 - 1- F (p 2 ;M ll M 2 



(C. 



-g 



2A 2 - Ml In - M 2 2 In - (M 2 - M 2 2 )F A {p 2 - M 1} M 2 ) 



(4vr) 2 



(4vr) 2 

where F , F4 and F 3 are defined by 



- In A 2 - F (p 2 ; M u M 2 ) + 4F 3 (p 2 ; M u M 2 ) 



(C.9) 



F (s; M l5 M 2 ) = / obln [(1 - x)M? + xM| - x(l - x)s] , 
Jo 

F A (s;Mi,M 2 ) = / ob(l-2x) In [(1 - x)M 2 + 2M 2 - x(l - , 
Jo 

F 3 (s;M l ,M 2 ) = / dxx{\ - x) In [(1 - x)M 2 + xM 2 - x(l - x)s] . 

JO 



(CIO) 



We consider the quantum corrections denoted by n( vac ) /u ' [see Eq. (|3.76j) ]. At tree level the 
two-point functions of A^, V M and V M are given by 

r (tree)/*v/ \ _ r 2 



n 



nf ree)A "» 



F 2 barc ^^ + 2^ barc (p 2 ^-pV) , 
F 2 ^ + 2z 1>bare (p 2 ^-^V) , 



9L 



n 



(tree)^ti' 

VII 



(P) 



^ 2 barc^ + ^3,bare(pV^-pV 



(C.11) 



Thus the one-loop contributions to H^ ac ' tlv give the quantum corrections to F 2 and 2:2- Simi- 



larly, each of the one-loop contributions to TT|| VtiL;Ail/ , jj^'v anc i n^[j' WAil/ includes the quantum 



j-(vac) \lv jj(vac)/ii/ an( j jj(vac)/u/ 

corrections to two parameters shown above. For distinguishing the quantum corrections to two 
parameters included in the two-point function, it is convenient to decompose each two-point 
function as 

jjCvac)/!./^ = u (vac)s (p)g^ + n (vac)LT (p)(^V - (C.12) 

It should be noticed that, since we use the Lagrangian with Lorentz invariance, the form of the 
quantum corrections is Lorentz invariant. Then, the relation between four components given 
in Eq. ()5.17|) and two components shown above are given by 



pj(vac)t pj(vac)s pj(vac)5 

■Q(vac)L -Q(vac)T jj (vac) LT 



(C.13) 



Using the decomposition in Eq. (|C.12J) . we identify n^f^ op ){p 2 ) with the quantum correction 



to F 2 , n^ ^L^ 2 ) with that to z 2 , and so on. It should be noticed that the following relation 
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is satisfied jHH EB] : 



n (vac)5 (p2 



n 



(vac) 5 



(P 2 



(C.14) 



Then the quantum correction to F 2 can be extracted from any of nf^, Uy" and IlyJji. 

We note that in Refs. (EH EH the finite corrections of are neglected by assuming 

that they are small. In this paper we include these finite contributions in addition to the 
divergent corrections. As in Refs. [HUE]; we adopt the on-shell renormalization condition. 
They are expressed as * x : 



Re 



n 



(vac)S 



(p 2 = o) 



Re 
Re 



n(vac)S/ 2 
V (P 



Mi 



*?(A* = 0) 
F CT 2 (/i = M t 



pi- 



ll 



(vac) LT ( 2 
V 



(P 2 = M) 



g\ti = M p 



(C.15) 
(C.16) 
(C.17) 



where \i denotes the renormalization point. Using the above renormalization conditions, we 
can write the two-point functions as 



n 



(vac)S/ 2 



(P 2 



r4 vac) V 



n 



(vac) LT ( 2 
V 



(P 2 



^ 2 (o) + ni(p 2 ), 
f 2 (m p ) + u v (p 2 ), 

1 + nf ( P 2 ) 



9 2 (M P 



(C.18) 



where 11^ (p 2 ), H v (p 2 ) and H v T (p 2 ) denote the finite renormalization effects. Their explicit 
forms are listed below. From the above renormalization conditions they satisfy 



mP = 0) = RelW = Mi) = RelRV (p 



Mi 



. 



(C.19) 



Thus in the present renormalization scheme, all the quantum effects for the on-shell parameters 
at leading order are included through the renormalization group equations. 



^In the framework of the ChPT with HLS, the renormalization point fi should be taken as fj, > M p since the 
vector meson decouples at [i — M p . Below the scale M p the parameter F„, which is expressed as F„\fi) in 
Refs. j^El, runs by the quantum correction from the pion loop alone. Then F%(/j, = 0) in the right-hand-side 
of the renormalization condition Eq. (|C.15|) is defined by [F„\(j, = 0)] 2 . Due to the presence of the quadratic 
divergence F 7r {M p ) is not connected smoothly to F^\M p ). The relation between them is expressed as |64II13) 
[F^ ] {M p )] 2 = F%{M p ) + [N f /(4n) 2 ][a(M p )/2}M^ where a(ji) for fi > M p is defined as a(ji) = F*(jjl)/F*{h). 
By using this relation, the renormalization condition (|C.15|l determines the condition for F 2 {M p ). Strictly 
speaking, we should use F. K {M p ) in the calculations in the present analysis. However, the difference between 
^(0) and F. K {M p ) inside the loop correction coming from the finite renormalization effect is of higher order, 
and we use ^(0) inside one- loop corrections in this paper. 
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In the following, using the above functions, we summarize the quantum corrections to 
two components n( vac ) 5 and n^ vac ^ LT of the two-point functions n^ ac - )Ml/ , ~\j(™c)vv ^ jjjyac)^ 
11^" which are defined in Eq. (|C.12J) . For A^-Ay two-point function, we obtain 



ii ±(l-loop)^ I 



N f 



n 



(vac)LT / 2 



(P 2 ) 



(4vr) 2 
N f a 



— A 2 + -a M 2 In A 2 + -a M 2 In M 2 

2 4 p 4 p p 



±(l-loop)u- / 4 

Corrections to V^-V u two-point function are given by 



a M 2 { 1 + F (p 2 ; M p , 0) + -F A (p 2 ; M p , 0) 



- In A 2 - F (p 2 ; M p , 0) + 4F 3 (p 2 ; M p , 0) 



(C.20) 



L ||(l-loop)< 



(47T) 



A 2 + 7 M 2 InM 2 + M 2 <| 7 InM 2 - 1 - F (p 2 ; M p , M P/ 



n (vac)Lr ^2 
li ||(l-loop) 



(p 2 ) 



Nf 1 



(4vr) 2 8L 



a 2 - 4a + 5 



lnA 2 -F (p 2 ;M p ,M p , 



+ 4F 3 (p 2 ; M p , Mp) - 4(2 - a) 2 InM 2 



(C.21) 



As for Vn-V u we have 



n^,, i—i^ ) — ii ||(i-ioop)Vi' )■> 



n 



L y (i-ioop) ^ 

(vac)LT / 2 



V(l-loop) 



G* 2 ) 



(47T) 



7 lnA 2 + 4 + ^Fo(p 2 ;M p ,M p ) 



24 
9 



F 3 (p 2 ; Mp, Mp) - - F (p 2 ; 0, 0) - 4F 3 (p 2 ; 0, 0) (C.22) 



Finally, corrections to V p-Vv two-point function are expressed as 



ii y||(l-loop)^ I 



_ n (vac)S , 2 \ 
ii ||(l-loop)^ li 



n 



(vac) LT ^2 
V||(l-loop) 



1 



(4tt) 



1 + 2a - a" 



lnA 2 -a(2-a)lnM 



F (p 2 - Mp, Mp) + 4F 3 (p 2 ; M p , M p ) . (C.23) 



The renormalization conditions in Eqs. (jC.15|) - (|C.17|) lead to the following relations among 
the bare and renormalized parameters: 



,bare 



4(4vr) 



■[2(2-a)A 2 + 3aM 2 lnA i 



^-^aM^M^ + l), 
^[(a 2 + l)A 2 + 3M 2 lnA 2 ] 



(C.24) 



N 



F 2 

cr,bare ~\2 



4(4vr 
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Appendix D 



Functions 



In this appendix, we list the integral forms of the functions which appear in the expressions of 
the physical quantities and the several limits of the loop integrals shown in Appendix El The 
functions I n (T) and J^(M;T) (n, m: integers) are given by 

I(T)= fJ^jk^ = ^ iT n 

I, = T , h= T5 , (D.D 



d 3 k 1 \k 



n—2 



where 

uj = Vk 2 + M 2 . (D.3) 

The functions F.^(p ; M;T) and G n (pQ]T), which appear in the vector meson pole mass in 
section EUl are defined as 

d 3 k 1 4\k\ n ~ 2 

F»(p ;M;T) 



(2tt) 3 e"/ T - 1 u(Au 2 - p 2 ) ' 
i 3 k \k\ n ~ 3 4\k\ 2 

(27T)3 e fc /r _ l 4 |fc|2 _ po 2 



d 3 k \k\ n -' 3 A\k\ 2 frs s 

Gn( Po ;T) = / — ' ' -' 1 - . (D.4) 
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Appendix E 

Temporal and Spatial Parts of 
Two-point Function 



In this appendix, we show that the temporal and spatial components of the two-point function 
jj(vac)/ii/ independent of the external momentum in the VM limit. 
From Eq. (|5.60j) . ~\j(™ c ^ h ^ u [ s rewritten as 



t 

r(vac)(b)/^ / — \ , T u 



Il^^M = N f jX£X»B(™)i*( Po ,p-M p ,0) 



= N f jB(™^(p o ,p-M p ,0), (E.l) 

where we define 

XZ = u»u-, + V*(g£-u»u-,). (E.2) 

In the above expression, we define the function £?( vac )/^ contributed to the diagram (b) in 
Fig. 15.11 as 



5 {vac) ^(Po,p;M p ,o) 

dk f d 3 k {2k-pY{2k-p) 1 



with 



i(2k)J (2tt)3 [kl - uflftko - po) 2 - (com 



"1 = v 2 l\ 



(E.3) 



K? = V a 2 \k-p\ 2 + M 2 p . (E.4) 
In terms of each components of _g( vac )^ ; the temporal and spatial parts of B^ u are given by 
B^\p ,p;M p ,0) 
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£ (wc)t (po,p; M p , 0) + (1 - V 2 )^B^ h (p ,p; M p , 0) 



B^ s {p ,p;M p ,0) 



B^ s (p Q ,p;M p ,0) 



V 

1 - v 2 - 2 



^5 (vac)L (p ,P;M p ,o) 



(E.5) 



V 2 p 2- 

By using the expressions in Appendix El the componets B^ vac ^, _E?( vac ) s and 5( vac ) L take the 
following forms: 

B< TOC > t ( P o,p;M p ,o) 



-1 (2^-p ) 2 -1 + 



(27r)H2^(^-p ) 2 -K) 2 2u p p (u p p+Po y-uj 2 
p-(2k-p) \ 1 2^-po | 1 + 

\ 2uv (a^ - p ) 2 - (u;p) 2 2a$ {u p p + p ) 2 - ^ 2 



Po 



£^ c > s (po,p;M p ,0) 



(2vr) 3 



(2& • p — p 

^2 



+ 



2^ - po) 2 - (^p 1 ) 2 + po) 2 - iol 



+ 



Pop ■ (2k -p) \ 1 2^y - po 



P 2 



B^ L (p o ,P;M p ,0) 

d 3 k p 2 p ■ (2k — p) f 1 



\ 2u n (u n - po) 2 - H) 2 ^ 2u p (cu p p + po) 2 - ul 
2uj w - po 1 2lj p p + po 



(27T) 



Pop' 



_2lv 7T (u n - po) 2 - (u p p ) 2 2u p p (u/p + po) 2 - ^ 2 



(E.6) 



Now we take the VM limit. Note that the VM condition for a* and a s implies that a 
velocity becomes equal to the pion velocity, V a — > V v for T — > T c as shown in Eq. (|5.21|) . The 
components £?( vac ) t , _B( vac ) s and £>( vac ) L are calculated as follows: 

limS^^pjM^O) 

d 3 k -1 l(k;po,p) + p- (2k - p)J(k;p ,p) 



VM 



\imB^ c >(po,p;M p ,0) 



(2vr) 3 ^ [fa -po) 2 - K) 2 ][K+Po) 2 - M?) 2 ] ' 
p ,p;M p ,( 

d 3 fc 1 1 1 
(2^^ [K - Po) 2 - M?) 2 ] [K + Po) 2 - M?) 2 ] 
p- (2k-p)) 2 lC(k;poip) +PoP- (2k - p)J(k;p ,p) 

p- (2k - p)J(k;p ,p) 



(E.7) 



hmfi( vac ) L (p ,p;M p ,0) 
c/ 3 A; 1 p 2 



(27r) 3 ^p 2 [K - po) 2 - K) 2 ][(^ + p ) 2 - (^) 2 ]' 
where we define the functions as 

J(fc;p ,p) = u; 2 [4u; 2 - 4(^) 2 - 3p 2 ] +p 2 [p 2 - (^) 2 ], 



(E.8) 



(E.9) 
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J(k ; p ,p) = -p [-3ul+p 2 -(u,%) 2 }, 

/C(fc;po,p)=^ 2 +p 2 -(^) 2 . (E.10) 
Using Eqs. (|E.7J) - (|E.9J) with these functions, we obtain ^O^M as 

d 3 k -1 I(kp ,P)+V 2 p-(2k-p)J(k;Po,P) 



(E.ll) 



(2vr)3 ^ [fa - Po y - {uiY} [k + Po y - fay] ' 

d 3 A; 1 1 1 

(27p^v;V -po) 2 - K?IlK + po) 2 - K) 2 ] 

(KV • (2k - p)) 2 /C(fc; p , P) + pIVZp- (2k -p)J(k; p , p)] . (E. 12) 

The integrand of the functions 5^ vac ^ and _E?( vac ) s in the VM limit are same as those of B^ 11 ^ 1 
and _E?( vac ) s with = 1 when we make the following replacement in ^( vac ) t ' s : 

Vjc^k, V w \k-p\ -> (E.13) 

The functions 5 (vac)t,s with M p -> are independent of the external momentum po and p [see 
Appendix |B] #1 . Thus we find that the functions B^^ and i?( vac ) s in the VM limit are 
obtained independently of the external momentum p and p: 

U m flM.( Po , p;M,,0) = 

Mm £<"">■( p„,fi; M„,0) = -V,^. (E.14) 

VM r 87T 



^ tl In Ref. where = 1 was taken, it was shown that the hadronic corrections n^ ;5 (jio,p; T) at the VM 
limit are independent of the external momentum p$ and p. The structure of the integrand in the vacuum part 
is the same as that in the hadronic part except for the absence of the Bosc-Einstein distribution function. Thus 
the vacuum part is also independent of po and p. 



Appendix F 



Current Correlator in Operator 
Product Expansion 



In this appendix, following Ref. [ZO] where the current correlator is discussed using the OPE 
in dense matter, we determine the Lorentz non-invariant effect in bare pion velocity Kr,bare 
through the Wilsonian matching. The tensor structure of the current correlator with Lorentz 
non-scalar operators in dense matter is same as that in hot matter. Thus in order to obtain 
the current correlator in terms of OPE variables, we simply make a replacement of the matrix 
element in dense matter with the one in hot matter. 

The axial-vector current correlator obtained from the OPE in dense matter is given 
by [SI 



GT(?o,?) = (gV-<rg 2 ' 



-1 



1 'i + ^b/ * 



2tt 2 



7T 



1 /a 



[i 2 J 6Q 4 \ 7T 



G- 



2,7i(y it — \ 2y 



Ana. 



u,d,: 



Q 

J^J^-^W 2 + g^q v q^ + + g^g^Q 2 } 

2 2k 



T=2 k = l 



Q 



_ A2k+T,T 



EE 

t=2 k=l 



q 2 



Y2k 



... Q^k . 



Q 4 



2k + T,T 



(F.l) 



with Q 2 = —q 2 . r = d — s denotes the twist, where s = 2k is the number of spin indices of the 
operator of dimension d. When we consider the operators of (r, s) = (2, 2) and (2, 4), G^ has 
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A al and A a0\a defined by #1 



A 
A 



4,2 _ 
a0 ~ 

2 

a0\o 



PaPf3 ~ -TQapP 



-P-A 4 ' 2 
2m 2 ' 



V 



PaP(3P\Pa ~ -^(PaPpgXa + PaP\9pa + PaPa9\p + PpP\9c 



+ PpPa9 a \ + P\Pa9afi] 



J^{g a p9\a + 9aX9/3a + SW#/3a) ^-^4' 2 , ( F - 2 ) 
4o J LTfl 



where p denotes the nuclear density and m the nucleon mass, and is the nucleon momentum 



and we take p h 



mu h 



m(l,0). From Eq. (|F.2|) . we have the following expressions of G 



L(0) 
A 



and G 



L(l) 



Hi 



#2. 



G (OPE)L(0) ( _ g2 , 



-1 » yr r(OPE)(0) 



7T 



+ 



2tt 2 
27ra s 

IF 

~w 

pm 



In 



6Q 4 

. 2 



a . 



7T 



u,d,s 



G 



(OPE)L(l), 



2pm d 



-q 

AU 



4 u + d _ 1 



5 pm 



_ A 



u+d 



(F.3) 



(F.4) 



Following the same procedure done in section l4~2*j we obtain the matching condition on the 
bare pion velocity: 



v: 



bare 



1 2 P m3 A u+d 

A^ 4 ' 



where we define Go as 

Gq 



8tt 2 



1 + ^ 

7T 



+ 



G 



2vr 2 e G 2 ) p 



(F.5) 



7T ' 



; 1408 a s (qq) 



™\2 
P 



3 A 4 27 A 6 

* 1 Eq. (8) in Ref. [75] does not have the factor In the following, we explain the reason why this factor ^ 
appears in Eq. (|F.2() : In the linear density approximation, the matrix element of an operator A is expressed as 
(G\A\G) = (0|A|0) + -^(p\A\p) where \G) denotes the ground state of nuclear matter jJSJ. The first term on 
the right-hand-side is the vacuum expectation value, which vanishes for operators with spin. In the second term 
\p) denotes a nucleon state with momentum p. In the present analysis, we consider the operators with spin 
2 and spin 4 which are dominant to the Lorentz non-invariant effect. The matrix elements of these operators 
with spin s are expressed as A s = (G\A S \G) — ^(p\A s \p) = ^-A s . ^2,4 in Eq. (8) in Ref. [Sj correspond to 



4,2 



(PaPp ~ \g a 0P 2 )A^ 2 = (p a pp - lgai3P 2 )^A^ 2 as showen in 



^2,4 ■ In terms of A s , A a p is expressed as Ajg 

Eq. era . 

# 2 Thc expression of Eq. (|F.4(1 is obtained from Eqs. (11) and (14) in Ref. |7H], where we neglect the 0(a s ) 
terms in A^ 2 and A 6 a 2 Xa . 
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+ 



2pm d 5pm 3 

'■ti-o 7T. — Aa 



(F.6) 



A 4 • 2 A6 "4 

In the above analysis, we consider the OPE in dense matter. The quantity in hot matter 
corresponding to pm 3 in Eq. (jF.5|) is determined as follows: From Eq. ()F.1|) . the term of 
(r, s) = (2, 4) is expressed as [ZO] 



1 r (i) * 

q 2 U A fj, 



T =2,s=4 



-G 



Mi) I 

A k=2,s=4 



-2G 



A |t=2,s=4 



20pm 3 d 



-,2^4 ^4 



(F.7) 



On the other hand, the same matrix element in hot matter is expressed as jZHj 



Ml 

Afi |t=2,s=4 



■8i ? "l" ( ( u^D u D x D a u + d lfl D„D x D a d 



Q 



10 



16-3 



16vr 4 T 6 „ /m,\ 7r 2 m 2 T 4 



x 



63 



q^q v q x q a 



10 



m„ 



Q 



10 



{u^UyUxU^ — tracesjv44 +d , 



where the functions Bi, B 2 and £>3 are defined by 

6 ^ , f~2 2 1 

"2/ vr - z — — 7, 



Bi(z) 



71- 



%f d yyW 



e« - 1 
1 



63 



dyy* x/y 2 - z 2 — - - 



8tt 6 

^ +d (p 2 = lGeV 2 ) = 0.255. 



r 

m\T 2 
48 



(F.8) 



(F.9) 



The momentum dependent part of Eq. (|F.8|) is expanded around q = as 

5 1 5 1 



q fJj q u q x q a 



Q 



10 



(u^UuUxUv — traces 



16 (-g 2 ) 3 4(-g 2 ) 4 



(F.10) 



From Eqs. (jF~Sjl and (jF~T0|l . G^/ with (r, s) = (2, 4) takes the form 



1 r (i) ^1 

^ |r=2,s=4 



12 



"9 



2^4 



16vr 4 T 6 5 /m, 
63 3 vT" 



tt m^T fw^ 
10 2 V T 



, <T 2 



77^ 
T 



:F.ir 



Comparing this expression with Eq. (IF. 71) . we find the following correspondence: 



pm 



16tt 4 T 6 
63 



T 



n 2 mlT A 
10~ 



B-) 



T 



m\T 2 
48 



B x 



T 



;F.12) 
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Thus the expression in hot matter corresponding to Eq. (jF.5|) is obtained as 



l-K 2 



7r,bare 

1 6 

Go 5 



16vr 4 T 6 



63 V T 



A 6 



^ +d . (F.13) 



Taking = 0, Eq. ()F.13|) is rewritten into the same form as Eq. (|4.31|) . 



Appendix G 



Hadronic Thermal Corrections to 
Susceptibilities 



In this appendix we summarize the hadronic thermal corrections to the two-point functions of 
Ap-Au, Vfj,-V v , V M -V„ and F„-V„. 

The four components of the hadronic thermal corrections to the two point function of 
A^- A v ^ are expressed as 

ITi(p ,p;T) = N f {a-l)A (0,T)-N f aM 2 B (p ,p;M p ,0;T) 

+ N f ^B t (p ,p;M p ,0;T) , (G.l) 

U s ± (p ,p;T) = N f (a-l)A (0,T)-N f aM 2 B (po,P]M p ,0;T) 

+ N f ~B s (p o ,p;M p ,0;T) , (G.2) 

ni( Po ,p;T) = N f ~B L (p o ,p;M p ,0;T) , (G.3) 

til(p ,P;T) = N f ~B T (p ,p;M p ,0;T) . (G.4) 

The two components IT' and IP of hadronic thermal corrections to the two-point functions 
of Vp-Vu, Vp-Vu and V ^-V u are written as 

Tl t v (p Q ,P;T) = Tl s v (p ,P;T) 
= U\(p ,p ] T)=Tll(p ,p ] T) 
= -n* V ||(po,p;T) = -n^n(p ,p;T) 
= -JV/i [A (M p ; T) + a 2 A (0; T)] - N f M 2 B (p ,p- M p , M p - T) . (G.5) 

Among the remaining components only flf is relevant to the present analysis. This is given 
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by 



nf(po,p;T) = N f ±B L (p ,p-M p ,M p ;T) + N f ( 2 ^ B L (p ,p; 0, 0; T) . (G.6) 

For obtaining the pion decay constants and velocity in subection 13.4.21 we need the limit 
of Po = p of n* ± and U s ± in Eqs. flUH) and (j(T2jl . With Eq. (fCTTfi^ . lT ± and fl s ± reduce to the 
following forms in the limit M p — > and a — > 1: 

nKpo=p + <c,p;T) —> -^J 1 2 (0;T) = --fT 2 , 

M p — >0, a— >1 Z Z4 



ni(p =p + ie,p;T) 



-^^?(0;T) 



24 



(G.7) 



In the static-low- momentum limits of the functions listed in Eq. ()G.17|) . the (IT^ — II ^ 
appearing in the axial-vector susceptibility becomes 



lim lim [n* ± (p ,p; T) - U^(p 0l p; T 

= -N f J*(P; T) + N f a Jf(M p ; T) - N f £- [l^M,, T) - 3^(0; T) 

p 



(G.8) 



For the functions appearing in the vector susceptibility relevant to the present analysis we have 

N f 



lim lim [lT y (p ,p; T)] 

p^0p ^0 



lim lim riTM / (po,.p; 

p^Opo-»0 1 



2J°_ 1 (M p ,T)-J 2 (M p ,T)+a 2 J 2 (0;T) , (G.9) 



4 L 



M p 2 J 1 °(M p ;T) + 2J 1 2 (M p ;T) 



(2 - af 



J 2 (0;T) . 



(G.10) 



In the following, we list the explicit forms of the functions that figure in the hadronic 
thermal corrections, Aq, Bo and B^ u in various limits relevant to the present analysis. 
The function A (M;T) is expressed as 



where J 2 (M;T) is defined by 

J™(M;T) 



A Q (M-T) = J 2 (M;T) 



d 3 k 



\n-2 



(2tt) 3 e -(fc;M)/T _ i 



(G.n; 



(G.12) 



with / and n being integers and u(k;M) = y M 2 + \k\ 2 . In the massless limit M — 0, the 
above integration can be performed analytically. Here we list the result relevant to the present 
analysis: 



J 2 (0;T) = .f^T) = ^T 2 



(G.13) 
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It is convenient to decompose B^ u into four components as done for in Eq. ([5.17)1 : 

= vTtfB 1 + {g" u - u»u u )B s + Pf ' B L + P£ U B T . (G.14) 
We note here that, by explicit computations, the following relations are satisfied: 

B\p^p-M,M-T) = B s {p,,p-M,M-T) = -2A (M;T) = -2J?(M;T) . (G.15) 



To obtain the pion decay constants and velocity in Section rj.4.2l we need the limit of po = p 
of the functions in Eqs. (|G.1J) and (|G.2jl . As for the functions M 2 Bo, B l and B s appearing in 
Eqs. (jHZQ) and (JUIJ), we find that, in the limit of M p going to zero, they reduce to 



M 2 p B (po=p + te J p;M p ,0;T) — > 



B t (p =p + ie,p;M p ,0-T) 



-2J 2 (0;T) = --T 



B s (p Q =P + ie,p;M p ,0;T) — ► -2^(0; T) = --T 



(G.16) 



The static-low-momentum limits of the functions appearing in the corrections to the axial- 
vector and vector susceptibility are summarized as 

lim lim [M 2 p B (p ,p; M p , 0; T)] = -J%(M P ; T) + J 2 (0; T) , 

p->0p ->0 H 

lim lim & M p ,0;T) - B L (p Q ,p;M p ,0;T 



-4 



-J 2 _ 1 (M p ;T) + ^(0; T) 
lim lim [M 2 p B (po,p;M p ,M p -T 



p— *Opo— >0 



^°i(M p ; T) — J\{M P \ T) 



lim lim [5 L (po,P; M p , M p ; T)] = -2M 2 J?(M p ; T) - 4J 1 2 (M p ; T) 

p^Opo— +U 

lim hm [5 L (po,p;0,0;T)] = -4J?(0;T) . 

p^Opo^O 



(G.17) 



Appendix H 

Explicit Calculation of Quantum 
Correction to AM 



In this appendix, we compute the quantum effects on the masses of 0~ (P) and + (Q*) heavy- 
light .M-mesons by calculating the one-loop corrections to the two-point functions of P and 
Q* denoted by Upp and Hq*q*. Here we adopt the following regularization method to identify 
the power divergences: We first perform the integration over the temporal component of the 
integration momentum, and then in the remaining integration over three-momentum we make 
the replacements given by 



A 



d 3 k 1 



(2tt)3 k 2 2V2n 2 



d 3 k 1 
(2tt) 3 I 



A 



d 3 k 



A 



12>/2- 



(H.1) 



7T- 



Here we use the t'Hooft-Feynman gauge for fixing the gauge of the HLS. 

The diagrams contributing to Upp are shown in Fig. IH.lt In the limit of zero external 



p 

(a) 



(d) 



P P • 



P P 



(b) 
n 

Q* 
(e) 



• P P 



P P 




(f) 



Figure H.l: Diagrams contributing to P-P two point function. 
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Figure H.2: Diagrams contributing to Q*-Q* two point function. 



momentum, the divergent parts of these contributions are given by 



n 



[a)[aP] 
PP 



llpp 



div 



div 



n 



(c)[pP) 

PP 



2tf 

n l 

2tf_ 

9 2 , 



(A 2 - 2M p lnA) + ^(^-M H \nA 
(4vr) 2 V p > 4vr 2 Vv^ 

A 3 4 9 Ml f A 

A 2 + ff 



24 v / 2tt 2 (4tt 



A 



4vr 2 Vv 7 ^ 



M^lnA 
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(d)[trO] 

PP 
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2tt 2 
2fc 2 



v/2 
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, A2 _, 2l a \ , 



„ (A 2 -2M 2 lnA) + 
F 2 [24V/2VT 2 (4vr) 2 ^ p J 4vr 2 



M 2 
p 



n( e )[ 7rC 
iipp 



n 



(/)[pQ] 
pp 
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A2 + M 2 



2P 

7T 

^!fc 2 I — M G In A 
2tt 2 V.V2 



(4tt 
A 



A 



4tt 2 



M G In A 



fH.2) 



The particles that figure in the loop are indicated by the suffix in square bracket; e.g., [7rP*] 
indicates that n and P* enter in the internal lines. Here and henceforth, we suppress, for 
notational simplification, the group factor C 2 (Nf) defined as (T a )ij(T a )ji = C 2 (Nf)Su. 

The relevant diagrams contributing to Hq*q* are shown in Fig. IH.21 The divergent parts 
of these contributions in the low-energy limit are expressed as 



n 



n 



{a)[*Q* 
Q*Q* 

(b)hQ] 
Q*Q* 

tt(c)[pQ*. 



div 



div 



div 
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2k 2 r 

2k 2 



M, 
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(4tt) 2 
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(A 2 -2M 2 lnA) + g(A_ MGlnA ^) 



M G a2 + M 2 



F 2 [24V27T 2 (4vr) 2 
J (l-k) 2 (-±=-M G lnA 
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47T 2 



Mg In A 



2tt 2 
2P 
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- A^(A 2 -2M 2 lnA) 



Ml 



M 2 ( A 
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Q*Q* 



n 



imp* 

Q*Q* 



div" 

= 3£ 

div 27T 2 



-A 2 + — | - M H In A 
(4tt) 2 4vr 2 Va/2 



A- 2 



A 

7! 



In A 



Now, let us compute the difference of Hq*q* — Upp- 



It is easy to show that llpp 6 " 1 



given in Eqs. ()H.2|) and (|H.3J) . we have 



dh 



exactly cancels with H^q*q* 



div 



(H.3) 

From the explicit forms 



n 



(6)[ttQ] 



n 



n 



n 



(e)KQ1 
PP 

(&)[ttP*] 

PP 



2k 2 A 3 
div ~ F 2 24^2 
2fc 2 A 3 



(H.4) 



div F 2 24^/2tt 2 ' 

Note that the logarithmic, linear and quadratic divergences in I1q*q* are exactly canceled by 
those in llpp. This cancellation simply reflects that the external particles are chiral partners. 
This immediately leads to 



u (b+e) 
LL Q*Q* 



n 



(b+e) 
PP 



. 



div 



(H.5) 



The cubic divergence in Hq*q* is exactly canceled by that in Upp, reflecting that the internal 
particles are chiral partners to each other. 

In a similar way, a partial cancellation takes place between T1q*q* anc ^ ^pp as wen as 
between HqIq* and Hp P : 



n 



n 



(a)[aQ*] 


llpp 




2fc 2 


A 3 


Q*Q* 


div 


*F 


24V2n 


(,d)[*P*] 


(a)[<xP] 
llpp 




2k 2 


- A 3 


Q*Q* 


div 


F 2 


_24y/2n 2 




jr{a+d) yt(a+d) 
1L Q*Q* 1L PP 


div 


k 2 



Ml 



+ 



Ml 



—= - Mq In A 



— - M H In A 



These lead to 

where we used M 2 = g 2 F 2 . The remaining contributions sum to 

-9 



rr(c+f) _ rj(c+f) 
Q*Q* L1 PP 



2k 2k2 -(M G -M H )\nA 



(H.6) 
(H.7) 

(H.8) 



div 27T 2 

By summing up the contributions in Eqs. (|H.5j) . ()H.7|) and (|H.8j) . we obtain the divergent 
part of the correction to the mass difference: 

,2 



n * 



Q*Q* 



n 



pp 



div 



-C 2 (N f ) ^(l-2k- A; 2 ) (M G - M H ) In A 



(H.9) 



where we reinstated the group factor C2{Nf). The logarithmic divergence in the above expres- 
sion is renormalized by the bare contribution given by 



IT 



n 



pp 



AM, 



bare 



bare 



(H.10) 
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Thus the RGE for the mass difference AM = Mq — M# has the following form: 

^=*™£( i -*-*'K < H11 > 

We should stress that this RGE is exactly the same as the one in Eq. ()6.32j) obtained by setting 
a = 1, i.e., F a = F n . 

Note that the complete cancellation of ir- and u-loop contributions occurs even when we 
take into account finite renormalization effects. 
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